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1 .1  Molecular  Been  Epitaxy* 

This  is  e  research  prograa  involving  aany  aspects  of  senieondnetor 
research  and  aolecnlar  bean  epitaxy.  The  research  effort  includes  materials 
growth  by  MBE.  studies  of  basic  electronic  and  optical  properties  of  layered 
structures  and  applications  to  high  speed  electronic  devices.  The  objective 
is  to  utilise  the  unique  features  of  MBE  in  the  preparation  of  novel 
structures  and  to  explore  new  opportunities  for  devices. 

In  this  period,  the  research  on  Al^Gaj.jAs/GaAs  single  and  aultiple 
interface  structures  continued.  In  addition,  we  began  studying  the 
incorporation  of  Bi  both  in  GaAs  and 

Multiple  iaterface  he teroj unctions  with  GaAs  layers  as  thin  as  25  &  were 
grown  in  a  wide  range  of  growth  teaperatures.  The  photolnainescence  intensity 
associated  with  the  quantua  wells  was  very  low  for  substrate  growth 
teaperatures  below  about  650°C  and  iaproved  steadily  as  the  substrate 
temperature  increased  to  700°C.  The  principal  electron  to  heavy  hole 
transition  appeared  to  have  shifted  downward  in  energy  unless  the  structure 
was  grown  at  abont  700°C.  Above  700°C,  the  emission  froa  the  quantua  wells 


*  This  work  was  supported  by  the  Joint  Services  Electronics  Prograa  (U.S. 
Army,  U.S.  Navy,  and  U.S.  Air  Force)  under  contract  N00014-79-C-0424. 
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diminished  entirely  while  that  from  the  Al^Gaj.jAe  improved  steadily  with 
increasing  substrate  temperature.  The  overall  conclusion  of  this  study  is 
that  the  heterointerfaoial  properties  are  very  important  and  dominate  the 
luminescence  spectra.  Investigations  of  the  single  interface  structures  with 
the  binary  on  top  of  ternary,  and  the  ternary  on  top  of  binary  shoved  a 
substantial  anisotropy.  While  it  ie  relatively  easy  to  obtain  good  interfaces 
when  the  Al^aj^^As  layer  is  grown  on  top  of  GaAs.  the  reverse  is  not  always 
reproducible  and  requires  a  precise  set  of  growth  conditions,  the  most 
dominant  growth  parsmeter  being  the  snbstrate  temperature.  At  a  substrate 
temperature  of  700°C,  the  atomic  surface  morphology  of  the  Al^a^^As  is  best, 
leading  to  a  specular  heterointerfaee.  This  interface  roughness  is 
responsible  for  the  low  Pi-  intensity  as  well  as  the  shift  in  energy. 

The  rough  surface  of  the  Al^Ga^^As  layer  can  be  improved  by  growing  a 
thin.  e.g..  50  X.  GaAs  layer.  As  a  resnlt  of  this  smoothing  effect*  the 
multiple  interface  structures  are  espeeted  to  be  a  little  less  sensitive  to 
the  growth  conditions  than  the  single  well  heterojunetion  structure.  The  use 
of  a  smoothing  layer  to  relax  the  growth  conditions  without  the  loss  of 
quality  was  demonstrated  using  GaAs  FBTs  on  ^Oaj.jAs  buffer  layers. 

Incorporation  of  donor  impurities  in  A1zGa1_zAs  was  also  investigated. 
Particular  emphasis  was  placed  on  Si  because  the  other  most  commonly  nsed 
donor.  Sn,  evaporates  quite  readly  from  the  surface  at  high  substrate 
temperatures.  Both  electrical  measurements  and  secondary  ion  mass 
speotroseopy  techniques  were  used  to  characterixe  the  films.  A  snbstrate 
temperature  range  of  600  -  720°C.  mole  fractions  between  0  and  1.0  and 
electron  concentrations  of  10**  -  5  x  10**  cm“*  were  investigated.  It  was 
discovered  that  caution  must  be.  exercised  in  determining  the  donor  activation 
energy  from  the  n  vs  T  measurements.  The  modulation  doping  effect  can  lead  to 
anomalously  high  activation  energies.  Unlike  Sn.  the  activation  energy  of  Si 
in  was  small.  ~  15  eV.  except  for  some  deepening  at  about  x  *.5. 

1.2  Modulation  fiflaiA  AIxG»i-xA*/GaAs  FETs* 

Selectively  doped  Al^a^gAs/GaAs  hetero structures  where  only  the  larger 


*  This  work  was  supported  by  the  U.  S.  Air  Force  Office  of  Scientific 
Researeh  under  contract  AFOSR  80-0084. 
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bandgap  material  ia  doped  with  a  suitable  donor  are  of  speoial  importance  from 
the  standpoint  of  studies  of  basie  phenomena  and  device  applications.  The 
objectives  of  this  program  are  to  study  in  detail  the  low  and  high  field 
transport  parallel  to  the  heterointerface  and  to  fabricate  short  ehannel  FETs 
for  switching  applications. 

During  the  last  period  our  investigation  of  single  and  multiple  interface 
undoped  and  modulation  doped  GsAs/Al^Gs^_xAs  heterostruotures  was  undertaken. 
This  program  differs  from  our  other  work  in  that  the  objective  is  geared 
toward  FETs.  Primarily *  the  electrical  properties  of  single  interface 
modulation  doped  structures  and  optical  properties  of  undoped  multiple 
interface  structures  were  studied.  In  addition*  by  measuring  transport 
properties  parallel  to  the  heterointerface  for  a  number  of  structures  grown 
under  different  conditions*  the  optimum  MBE  growth  conditions  were  determined 
for  both  binary-on-ternary  and  terns ry-on-binary  structures.  Structures  with 
the  ternary  on  top  of  binary  exhibited  much  better  performance  and  thus  were 
used  to  fabricate  FETs.  The  thic’ less  of  the  undoped  separation  layer  in 
these  structures  was  adjusted  to  provide  an  optimum  donor/electron  separation 
for  FET  operation. 

Optimized  normally-off  modulation  doped  FETs  with  3  pm  channel  lengths 
and  1  pm  gate  lengths  displayed  transconduetsnees  as  high  as  230  mS/nm  at  300 
K.  This  figure  is  expected  to  double  at  77  K.  It  is  predicted  that  such  a 
FET  will  exhibit  a  switching  time  of  8  psec  in  a  ring  oscillator  configuration 
at  300  K.  Needless  to  say*  operating  parameters  such  as  load  resistance  and 
power  supply  voltage  must  be  optimized  to  aohieve  this  performance. 
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2.1  Introduction 

Tka  prlaary  goal  of  tkla  raaaarek  ia  to  evaluate  tka  potential  of  higk- 
purity  IaP  and  InGaAaP  (iaeiadiag  tka  ternary  liait  of  tkia  lattica-aatckad 
quaternary  alloy  ayataa)  for  optical  and  microwave  device  applieatioaa.  Tka 
tkrnat  of  tka  JSEP  anpported  work  is  to  investigate  tka  hydride  growtk 
teckniqne  for  tka  growtk  of  IaP  aad  InGaAsP.  This  process  is  partiealarly 
iatarastiag  baeaasa  it  peraits  easy  control  of  tka  III-V  ratio,  ia  eoatraat  to 
tka  PClj  (cklorida)  growtk  technique,  aad  baeaasa  tka  growtk  of  P  containing 
eoapoaads  aaaas  to  be  aneh  aora  controllable  wkan  contrasted  to  tka  popular 
MOCVD  aad  MBB  growtk  taekniqnas.  With  otkar  anpport.  LPE  growtk  of  IaP  aad 
InGaAaP  ia  alao  being  studied,  and  tka  properties  of  tka  hydride-VPB  aad  LPE 
epitaxial  layers  are  cosipared.  Tka  techniques  used  ia  this  work  to 
characterise  the  material  include  tka  measurement  aad  analysis  of  Hall 
eoeffieieata  aad  reaistivity  measurements  over  the  4-350  K  temperature  range 
to  determine  tka  total  eoaeeatratioas  of  electrically  active  donors  aad 
acceptors.  Tka  influence  of  the  growth  conditions  on  the  identity  aad 
concentration  of  the  acceptor  species  iaeorporsted  in  the  epitaxial  layers  is 
studied  through  photolumiaescence  measurements  over  the  1.5-20  E  range.  Deep 
levels,  aad  the  influence  of  different  growth  parameters  on  the  concentrations 
and  energy  levels  of  these  centers  are  studied  through  DLTS  measurements.  The 
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crystalline  quality  and  degree  of  lattice  natch  are  studied  with  an  x-ray 
diffraetoneter  purchased  with  partial  support  using  JSEP  funds.  The  residual 
shallow  donor  inpurities  in  InP  are  studied  using  far  infrared  Fourier 
transfom  neasurenents  of  photothernal  ionization  photoconductivity  on  high 
purity  epitaxial  InP. 

2.2  Hydride  Vapor  Phase  Growth  &£.  High  Purity 

The  su»st  connon  figure  of  nerit  cited  for  high  purity  InP  is  the  najority 
carrier  srobility  at  77  K.  The  carrier  concentration  and  the  nobility  are 
generally  derived  fron  the  Hall  coefficient  and  resistivity.  There  is  a 
danger  in  using  the  77  K  nobility  as  the  sole  indication  of  sanple  purity  in 
that  the  neasurenent  technique  averages  the  Hall  constant  and  sasqple 
resistivity  over  the  sasqile.  The  electrical  data  obtained  fron  the  van  der 
Pauw  technique  are  of  sone  use  in  detecting  inhonogeneities.  In  a  honogeneous 
sanple,  the  ratios  of  resistanes  neasured  on  a  van  der  Pauw-Hall  effect  sas^le 
by  pemuting  contacts  should  be  detemined  by  sanple  and  contact  geonetry 
only.  They  should  be  the  sane  at  77  K  and  300  (.  Inhonogeneities  in  inpurity 
concentration  across  the  wafer  can  cause  differences  in  these  ratios  at  the 
different  tenperatures.  When  proper  caution  is  exercised,  the  77  E  nobility 
is  a  useful  neasure  of  the  saaple  quality.  The  77  E  carrier  concentration  and 
the  conpensation  curves  calculated  for  InP,  [HI]  provide  an  estinate  of  the 
concentration  of  eleotrieally  active  donors  and  acceptors  in  the  crystal. 
Figure  2.1  shows  the  77  E  nobilities  of  a  nunber  of  VPE  hydride  grown  saarpies 
for  a  wide  range  of  carrier  concentrations  as  well  as  the  corresponding 
theoretical  curves  of  lode  [113.  The  data  points  generally  follow  the  curve 
for  (N*  +  lO/(n  +  p)  ■  2  indicating  that  the  concentration  of  donors  is 
approxinately  three  tines  that  of  acceptors.  By  contrast,  high  purity  InP 
grown  by  LPE  appears  to  be  nearly  unconpensated. 

Table  2.1  is  a  eonpilation  of  the  77  E  and  300  E  Hall  data  for  the  ten 
best  InP  layers  grown  by  VPE  in  our  lab.  Also  included  for  conparison  are  the 
highest  nobility  VPE  sanples,  both  halide,  reported  in  the  literature.  The 
71,000  cn^/V-see  layer,  VPE  1-221  is  the  highest  nobility  reported  to  date  for 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Progran  (U.S. 
Any,  U.S.  Navy,  and  U.S.  Air  Force)  under  contracts  DAAG-29-7 8-C-0016  and 
N00014-79-C-0424 . 
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Fig.  2.1  (top  loft)  77  X  aobllltloo  os  o  function  of  corrlor  concentration 
for  VPE  Hydride  InP  plotted  with  the  Rode  curves  for  InP. 


Fig.  2.2  (top  right)  Hell  nobility  at  77  and  300  X  versus  Nn-NA  of  the  InP 
layers  grown  in  this  work.  The  solid  curves  show  Rode* a  calcu¬ 
lation  of  u77  for  two  compensation  ratios. 

Fig.  2.3  (bottom)  Bell  s»bility  at  77  and  300  X  versus  of  the  InGaAs 

layers  grown  in  this  work. 
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i  sample  grown  by  the  hydride  technique.  It  is  somewhat  more  compensated  than 
high  parity  material  grown  by  LPE  and  the  300  I  mobility  is  not  as  large  as 
one  would  expect  for  a  layer  with  such  high  77  K  mobility.  Some  samples,  for 
example  VPE  1-205,  with  lower  77  E  mobility  hawe  larger  300  E  mobility. 

One  interesting  aspect  of  the  data  is  the  existence  of  an  aging  effect 
first  reported  by  Eastman  in  InP  grown  by  LPE  [R2] .  The  aging  phenomenon  in 
which  the  77  E  mobility  in  n-type  samples  increases  and  77  E  carrier 
concentration  simultaneously  decreases  after  3  to  4  months  storage  at  300  E, 
was  attributed  by  Eastman  to  silicon  donors  changing  from  shallow  to  deep 
states.  The  Hall  mobilities  and  oarrier  concentrations  of  InP  grown  by  the 
hydride  technique  chsnged  similarly  when  stored  at  room  temperature  (aee  Table 
2.2).  Dramatic  increases  have  been  observed  in  as  short  a  time  as  one  week  in 
which  the  77  I  mobility  changed  from  56,800  cm2/V-sec  to  70,030  cm2/V-aec. 
The  greatest  overall  change  occurs  in  the  material  which  atarts  ont  with  the 
highest  mobility.  In  the  eonrse  of  3  months  the  77  E  mobility  of  VPE  1-221 
increased  from  71,260  cm2/V-sec  to  95,920  em2/V-see.  The  process  can  be 
accelerated  by  annealing  at  105°C.  The  original  77  I  mobility  and  carrier 
concentration  are  re-established  when  the  sample  ia  heated  to  temperatures  as 
low  as  400°C  (used  for  alloying  contacts) .  All  values  for  mobility  and 
carrier  concentration  quoted  above  were  taken  immediately  after  the  contacts 
were  alloyed  except  where  otherwise  noted. 

2.3  I£E  Guaxth  slIaEiM  IbQsAsP* 

Methods  for  obtaining  high-parity  InP  and  InGaAsP  alloy  epitaxial  layers 
on  ( 100) -oriented  InP  substrates  using  liquid  phase  epitaxy  have  been 
explored.  The  300  and  77  E  electron  mobilities  for  InP  and  InGaAs  are 
smumarixed  in  Figs.  2.2  and  2.3.  All  of  the  undoped  layers  were  n-type,  and 
the  InP  liquid  nitrogen  mobilities,  when  compared  with  the  theoretical  curves 
by  Rode  [Rl]  for  varying  compensation,  show  relatively  little  compensation. 
The  freexe-out  ratios,  (300  E)/Njj-Na  (77  E) ,  of  the  InP  and  InGaAs 

*  This  work  was  supported  by  the  National  Science  Foundation  under  con¬ 
tracts  NSF  BCS-79-17581  and  NSF  DMR-80-20250-73 ,  the  Naval  Research  Laboratory 
under  contracts  N00173-78-C-0129  and  N00172-79-C-0184,  and  the  Joint  Services 
Electronics  Program  (U.S.  Army,  U.S.  Navy  and  U.S.  Air  Force)  under  contracts 
DAA8-29-7 8-C-001 6  and  N00014-79-C-0424. 
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layers  are  generally  between  1.0  and  1.2.  The  nobilities  of  these  InGaAs 
layers  grown  on  (100)  InP  are  in  fair  agreement  with  the  nobilities  of  InGaAs 
layers  grown  on  (lll)B  InP  by  Oliver  [R3]  with  the  exception  that  for  equally 
high  liquid  nitrogen  nobilities  the  roon  temperature  mobilities  of  these 
layers  are  lower.  Oliver's  highest  roon  temperature  nobility  was  13,000 
cm2V"*s"*  and  in  this  work  the  highest  is  11,000  cm2V"*s"*.  The  InGaAs  layers 
with  Np-N^  for  (2-4)  x  1014  cm”3  and  liquid  nitrogen  mobilities  above  60,000 
cm^V^s"2  were  grown  at  687. 5°C,  and  the  rest  were  grown  at  627. 5°C t  In 
general,  the  growth  of  InP  and  InGaAs  above  650°C  has  resulted  in  higher 
purity,  probably  due  to  a  lower  distribution  coefficent  for  silicon  at  higher 
growth  temperatures.  The  InGaAs  layers  with  approximately  equal  to  1  x 
10*4  cm-3  were  grown  from  unbaked  solutions'  and  those  with  >  1  x  1014 
were  doped  with  tin.  Typically,  InP  layers  grown  from  unbaked  solutions 
yielded  doping  levels  closer  to  1  x  lO2^  cm-3  presumably  due  to  the  use  of 
lower  quality  InP  source  material. 

The  ability  to  obtain  high  purity  layers  has  been  highly  reproducible. 
At  the  present  time  11  of  the  last  13  InP  layers  grown  using  this  technique 
have  had  liquid  nitrogen  mobilities  between  50,000  and  70,000  cm2V~1s~1  with 
Njj-na  between  9.0  x  1014  2.5  x  10* 5  cm"3.  Also,  12  out  of  the  14  InGaAs 
layers  grown  at  687 .5°C  have  had  liquid  nitrogen  mobilities  between  55,000  and 
70,000  cm2V"1s"1  with  Njj-Na  between  1.5  and  4.5  x  1014  cm"3.  Fewer 

*B0. 83 ®*0.1 7^*0. 40*0. 60  ^  “  1.15  |im)  layers  were  prepared  for  van  der  Pauw 
measurements,  but  of  those  which  were  measured,  the  highest  purity  layer  was 

determined  to  have  N^-N^  *  5,6  x  10* 4  cm"3  with  room  temperature  and  liquid 
nitrogen  temperature  mobilities  of  4,600  and  24,800  cm2V"*s~2,  respectively. 
Most  of  these  layers  were  grown  between  runs  where  device  structures  were 
grown  which  included  layers  doped  with  either  xinc  or  tin  on  InP  substrates 
doped  with  either  zinc  or  tin.  Thus,  these  results  are  highly  reproducible 
and  the  use  of  dopants  in  the  boat  has  not  had  an  adverse  effect  on  the  purity 
of  layers  grown  in  subsequent  runs. 
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Measurements  si  Shallow  Acceptor** 


Variable  temperature  (1.4  -  20  I).  high  resolution  (0.2  %) 
photoluaineaeence  aessureaents  were  performed  on  high-purity  undoped  InP 
layera  grown  at  550°C  uaing  a  III:V  gaa  ratio  of  one  and  an  HC1  partial 
preaanre  of  0.02.  The  ezeiton  region  of  a  highly  eoapensated  sample  is  shown 
in  Fig.  2.4a.  The  sharp  line  strneture  is  typical  of  high-quality  InP  and 
includes  lines  produced  by  the  reeoabination  of  free  ezeitons  (FE) .  neutral 
donor-bound  ezeitons  (D°,X) ,  ionized  donor-bound  ezeitons  (D+,X) ,  and  neutral 
acceptor-bound  ezeitons  (A°.X).  At  high  ezeitation  intensity,  the  band-to- 
acceptor  region  of  this  saaple  (Fig.  2.4b)  shows  the  "two-hole"  (HI)  replica 
of  the  (A°.X)  doublet  for  Zn  acceptors  superiaposed  on  the  broader  emission 
bands  due  to  Zn.  A  donor-to-aeceptor  peak  involving  an  unidentified  (21  aeV) 
shallow  acceptor  and  a  sharp  line  due  to  an  ezeiton  bound  to  a  deep  center  are 
also  observed. 

At  low  ezeitation  intensity,  aeaaureaenta  of  the  band-to-acceptor  region 
of  the  saaple  1-221  (see  Table  1.1)  as  a  function  of  teaperature  permit  the 
identification  of  the  donor- to-aceeptor  (D°-A°)  and  conduction  band-to- 
acoeptor  (e-A°)  peaks  involving  two  different  acceptors  of  different 
ionisation  energies.  (See  Fig.  2.4c).  The  (D°-A°)  peaks  diainish  in 
intensity  relative  to  the  (e-A°)  peaks  as  the  teaperature  is  increased,  due  to 
theraal  ionization  of  the  donors.  This  trend  peraits  proper  identification  of 
the  various  peaks  when  they  overlap,  as  in  the  present  case.  The  ionization 
energies  of  the  two  acceptors  are  45.4  and  40.5  aeV.  which  are  those  for  Zn 
and  the  so-called  "A^*  acceptor  (which  aay  be  C  or  Mg),  respectively. 

The  deep  region  of  the  spectrua  (Fig.  2.4d)  shows  a  aeries  of  sharp  lines 
which  are  various  LO  and  TO-  phonon  replicas  of  the  FE  eaission,  superiaposed 
on  a  broad,  structured  eaission  of  unknown  origin  which  has  not  peviously  been 
reported.  Future  work  will  be  directed  toward  the  identification  of  sons  of 
the  unknown  spectral  features  discussed  here. 


*  This  work  was  supported  by  the  Aray  Research  Office  under  contract  DAAfi- 
29-82-E-0059.  the  National  Soienoe  Foundation  under  contract  NSF  DMR-80- 
20250-73.  and  the  Joint  Services  Electronics  Prograa  (U.S.  Aray.  U.S.  Navy, 
and  U.S.  Air  Force)  under  contracts  DAAG-29-78-C-0016  and  N00014-79-C-0424. 
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AND  DEVICES 


Fig.  2.4  (a)  Exclton  region  of  the  photolumineacence  spectrum  of  a 

highly  compensated  InP  sample  (N^qq  -  6.2  •  10lftcm  ,  frozen 
out  at  77  K.  (b)  Band- to-accep tor  region  of  this  sample  at 
high  excitation  intensity,  (c)  Band -to-accep tor  region  of 
sample  1-221  at  three  different  temperatures,  (d)  Deep 
regions  of  a  high-purity  n-type  sample. 
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2.5  jjmiss.  Transform  Spectroscopy  si  gfelllgV  Poaors* 

The  aechanisa  for  the  process  of  phototheraal  ionization*  represented 
scheaatically  in  Fig,  2.5,  oan  ba  understood  by  considering  an  electron 
initially  occupying  a  donor  gronnd  state.  The  temperature  of  the  sample  is 
4.2  K  so  that  few  electrons  are  in  the  conduction  band.  Absorption  of  a  far 
infrared  photon  of  appropriate  energy  oan  induce  a  transition  to  one  of  the 
donor  eneited  states,  from  which  the  still  bound  electron  nay  be  excited  to 
the  conduction  band  by  absorption  of  phonon(s),  contributing  to  the  electrical 
conductivity  of  the  sample.  A  photothermsl  ionization  spectrum  then  consists 
of  a  plot  of  photoeonductive  response  versus  photon  energy  or  wavenuaber,  and 
contains  peaks  at  the  various  donor  transition  energies. 

The  ground  state  energies  of  shallow  donors  in  InP  are  very  nearly 
hydrogenie  but  the  energy  of  the  donor  ground  state  energy  differs  slightly 
from  the  value  predicted  froa  hydrogenie  theory,  and  since  this  difference  is 
donor  speoies  dependent,  each  hydrogenic  transition  in  a  phototheraal 
ionization  speetrna  contain  several  closely  spaced  peaks.  Bach  peak 
corresponds  to  a  particular  donor  species  present,  and  has  an  amplitude  which 
is  a  aeasure  of  the  relative  concentration  of  the  donor  species.  The 
nultiplet  structure  is  aost  easily  resolved  at  high  magnetic  fields  (about  3 
T)  in  the  hydrogenic  ls-2p  (a  -  -1)  transitions.  For  this  reason,  the  spectra 
which  appear  in  this  paper  show  only  these  transitions.  The  phototheraal 
ionization  spectra  of  Fig.  2.6  show  the  donor  peaks  present  in  high  purity  VPE 
and  LPE  InP  grown  in  our  laboratory.  The  LPE  samples  typically  show  three 
donor  peaks  labeled  A,  B,  and  C.  The  amplitudes  of  these  peaks  are  related  to 
the  relative  concentrations  of  the  associated  donors.  For  exaaple  the 
photoresponse  of  the  sas^le  at  the  energy  of  peak  B  is  largely  due  to  the 
Stark  broadened  low  energy  tail  of  peak  C.  By  nuaerically  fitting  a  sum  of 
individual  peaks  of  the  characteristic  shape  with  adjustable  amplitudes  to  the 
spectra,  it  was  possible  to  deteraine  the  typical  relative  concentrations  of 
donors  A,  B,  and  C  in  LPE  samples  as  approxiaately  0.2  to  0.1  to  0.7.  The 
spectra  for  the  VPE  staples  show  only  the  peaks  A  and  C  in  the  very  different 


*  This  work  was  supported  by  the  Office  of  Naval  Research  under  contract 
N00173-C-0184  (Navy  Subc.  81-15  froa  Vashington  University)  and  the  Joint  Ser¬ 
vices  Electronics  Program  (U.S.  Any,  U.S.  Navy,  and  U.S.  Air  Force)  under 
contracts  DAAG-29-78-C-0016  and  N00014-79-C-0424. 
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•5  Two  representations  of  the  photothermal  ionization  process: 

(a)  illustrating  the  multlstep  photon-phonon (s)  excitation*  H 

(b)  illustrating  the  difference  in  ground  state  (Is)  energies 
for  3  donor  species. 


6  Photothermal  ionization  spectra  for  InP  grown  by  LPE  and 
hydride  VPE  (PH.-VPE). 
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relative  concentrations  of  0.7  to  0.3. 

Carefully  controlled  doping  experiments  have  been  performed  by  several 
research  groups  to  identify  various  donor  species  with  the  associated 
photothermal  ionization  peaks  for  GaAs.  Such  experiments  are  difficult 
becanse  of  the  small  range  of  donor  concentrations  between  that  of  the  purest 
material  which  can  be  grown  by  a  given  technique  and  that  for  which  impurity 
interactions  degrade  the  spectra  by  broadening  and  distorting  the  spectral 
peaks.  Until  recently,  nearly  all  InP  grown  was  of  insufficient  purity  to 
resolve  these  peaks  and  so  very  few  donor  identification  experiments  have  been 
performed  for  InP.  Using  intentionally  doped  LPE  grown  InP.  Stradling  [R4] 
has  identified  photothermal  ionization  peaks  for  Sn.  Si,  and  S..  This  work 
indicates  that  the  peaks  A  and  C  are  Si  and  S,  respectively,  but  these 
identifications  should  be  regarded  as  tentative. 
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3.1  Idxit  Daaaity  Tara  Traaaport  lgfl.  Pram  let  ^ 

Quasi-Ona-Diaaational  Conductor!* 

Aa  idaal  oaa  diaaaaioaal  eoadaetor  is  a  chain  of  atoais  ia  a  liaaar 
periodic  array.  Mary  liaaar  crystal a.  such  aa  NbSe^  and  TaS3,  approziaata 
this  idaal  quite  oloaaly  baeaaaa  tba  aetallic  boada  along  tba  axis  ara  strong 
vbila  tba  Van  dar  Vaala  forcaa  batvaaa  ebaiaa  ara  relatively  vaak.  Sack  one- 
diaaaaioaal  aatals  ara  aastabla  to  tba  foraatioa  of  aa  aaargy  gap  at  tba  Farai 

w*'r*T*otor»  *p#  bale*  a  eritieal  teaperature,  Tc  (~  215  K  for  TaSj) .  This 
eaasaa  a  aaparlattioa  aad  aa  assoeiatad  alaotroa  ekarga  daaaity  tava  (CDV)  to 
fora  with  a  wavelength  of  x/kp, 

Ezpariaaats  [R1-R5]  ka ve  shown  that  ekarga  daaaity  wares  oaa  boys 
cokaraatly  through  the  lattlea  ia  NbSa3  sad  TaSj  wkaa  tka  applied  alaotrie 
field  azeaads  a  eartaia  threshold  value  Rj.  Joke  Bardaaa  developed  a  aodal  of 
CDV  dapiaaiag  ia  which  tka  CDV  tuaaals  cokaraatly  through  a  piaaiag  gap  ia  a 
aaaaar  aaalogous  to  Zaaar  tuaaeliag  [R6.R7].  Ia  additioa  ka  applied  pkotoa 
assisted  tuaaeliag  theory  [R7]  to  predict  tka  frequency  dapeadaat  eoaductivity 
of  NbSa^ .  Ezpariaaats  [R8]  parforaad  at  UCLA  have  skowa  azoellaat 
quaatitativa  agraaaaat  with  tka  tuaaeliag  aodal  pradietioas  of  tka  field  aad 
frequeaey  dapeadaat  eoaduetivities  of  NbSe^  and  TaSj,  azeapt  that  there  saaas 

This  work  was  supported  by  tka  Joiat  Services  Elaetroaics  Prograa  (U.S. 
Any,  U.S.  Navy,  U.S.  Air  Force)  uadar  eoatract  N00014-79-C-0424  aad  by  the 
Natioaal  Sciaace  Fouadatioa  uadar  great  ECS  81-20038. 
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to  bo  oa  additional  "classical*  coaponaat  to  tbo  a.e.  conductivity.  However, 
tho  raaalts  of  tbosa  axporiaaats  have  alao  boaa  iatarpratad  aaiag  a  pataljr 
elaaaieal  overdamped  harmonic  oscillator  modal. 

Tho  priaoipal  aim  of  tbit  work  ia  to  parform  aa  azparimaat  for  wkiek  tka 
tunneling  aad  olassieal  modal a  giro  qualitatively  diatiaot  pradietioaa.  Ia 
tkia  way  wa  kopa  to  aitkar  prove  or  disprove  tka  taaaaliag  hypothesis.  If  tka 
taaaaliag  modal  proves  to  bo  baaioally  correct,  tkaa  matarials  amok  aa  TaSj 
may  ba  asafal  ia  radio  frequency  aad  mierowawa  applieatioaa  (e.g.,  aoaaitiwa 
dataetora  aad  mizars). 

Tka  azparimaat  parformad  at  UCLA  [U]  maaaariag  tka  d.e.  conductivity  of 
NbSa^  In  tka  praaaaea  of  largo  a.o.  aigaala  waa  iatarpratad  by  oar  group  [3.4] 
ia  tarma  of  tka  taaaaliag  modal.  Tka  azparimaat,  wkiek  kad  previously  boaa 
iatarpratad  aaiag  a  elaaaieal  modal*  akowed  azeallaat  qualitative  aad 
raaaoaabla  quantitative  agreameat  witk  Bsrdeea’s  taaaaliag  modal  eoaplad  witk 
pkotoa  aaaiatad  taaaaliag  tkaory.  Vary  raoaatly  field  aad  freqaeaey  dapaadaat 
raetifieatioa  measurement*  oa  TaSj  performed  by  our  group  kave  yielded  re salt a 
wkiek  qualitatively  eaa  easily  ba  azplaiaad  aaiag  tka  taaaaliag  kypotkaaia  bat 
are  difficult,  if  uot  impossible,  to  iaterpret  classically.  Ia  tka  aaar 
future  wa  plaa  to  make  a  qaaatitativo  oomparisoa  of  raetifieatioa  measurameats 
<*itk  taaaaliag  tkaory  pradietioaa  for  a  wide  raaga  of  bias  voltages  aad 
frequencies. 

3 .2  Pi aeharae  Annealint* 

Stadias  have  eoatiaaad  oa  tho  use  of  a  glow  disekarga  ia  Helium  to 
gaaarata  aa  alaetroa  beam  for  aaaealiag  ioa-implaated  soaieoadactor  samples 
[3.3].  A  aaw  diagaostie  tool  -  laser  roflaetemetry  -  has  boaa  developed  to 
allow  aeoarate  determ iaatioa  of  silieoa  sample  temperature  aad  ragrowtk 
velocity  [3.5]. 

Laser  refleetometry  eoasists  of  diraetiag  light  from  a  Helium-Neon  laser 
at  432. S  am  oato  a  semiconductor  suable  dariag  aaaealiag.  Tka  raflaetad  beam 
is  syaehroaously  dataetad  aaiag  a  vacuum  photodiode  aad  loek-ia  amplifier  to 
remove  baekgrouad  sigaal  due  to  blaek  body  radiatioa  from  tka  sample.  As  tka 


*  This  work  was  supported  by  tka  Joiat  Services  Eleetroaies  Program  (U.S. 
Army,  U.S.  Navy,  U.8.  Air  Force)  aadar  eoatract  N00014-79-C-0424. 
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silicon  is  Rested  by  tbs  eleetren  besa,  its  index  of  rsfesotion  inereases, 
increasing  tbs  reflectivity  of  tbs  sample.  The  relation  between  index  of 
refrsetion  and  temperature  for  silicon  is  known  [19],  so  a  relation  between 
temperature  and  reflectivity  can  be  calculated.  Thus,  tbe  temperature  of  a 
silicon  ssaple  esn  be  monitored  during  annealing. 

Refleetometry  also  provides  an  acourate  measure  of  tbe  regrowtb  time 
during  annealing  [RIO].  Tbe  amorpbixed  surface  layer  bas  an  index  of 
refraction  wbieb  differs  from  tbst  of  tbe  single-crystal  substrate.  Tbe 
discontinuity  in  refractive  index  censes  a  reflection  of  tbe  laser  beam  from 
tbe  amorpbous-single  crystal  interface,  wbieb  interferes  with  tbe  light 
reflected  from  tbe  sample  surface.  During  regrowtb.  tbe  interface  proceeds  to 
tbe  snrface  by  solid  phase  epitaxy,  causing  tbe  path  length  of  tbe  beam 
reflected  from  tbe  interface  to  decrease.  Thus,  tbe  phase  of  tbe  beam 
reflected  from  tbe  interface  eontinnally  changes  with  respect  to  tbe  phase  of 
tbe  besm  reflected  from  tbe  surface,  causing  maxima  and  minima  or 
"interference  fringes"  in  tbe  reflectivity  as  tbe  sample  regrows.  Ultimately, 
tbe  reflectivity  reaches  tbe  value  for  aingle  crystal  silicon.  Using  this 
technique,  tbe  regrowtb  rate  was  measured  as  a  function  of  electron  beam  power 
and  sample  temperature.  Our  results  for  annealing  large  area  samples  show 
excellent  sgreement  with  laser  annealed  samples  [Rll] ,  indicating  that  tbe 
discharge  anneal  system  shows  promise  an  an  alternative  to  laser  annealing 
ay stems. 

Studies  were  also  performed  on  tbe  .annealing  of  single  crystal  silicon 
with  a  polycrystalline  silioon  layer  grown  on  top.  As  deposited,  tbe  sables 
bad  sheet  resistivities  grester  than  200  0  per  square.  Our  study  shows  tbst 
with  tbe  discharge  system,  these  samples  can  be  annealed  yielding  sheet 
resistivities  of  15  0  per  square  and  less.  These  results  are  obtained  over 
large  area  samples  with  anneal  times  less  tbsn  five  seconds. 

3.3  ElUll-iiliilii  Etching* 

Recent  work  in  tbe  srea  of  plasma-assisted  etching  bas  been  directed 
toward  tbe  design  and  oonstruotion  of  an  RF  plasma  discharge  system.  Previous 


*  This  work  was  supported  by  the  Joint  Services  Electronics  Program  (U.S. 
Army,  U.S.  Navy,  U.S.  Air  Force)  under  contract  N00014-79-C-0424. 
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MiMtek  utilized  a  DC  hollow-cathode  discharge  to  g#  crate  etehiog  species 
vhila  the  new  design  employs  aa  inductively  coupled  Si'  is  charge.  This  offers 
several  advantages.  High  frequency  discharges  generally  tend  to  provide:  (1) 
a  higher  degree  of  ionization  and  dissociation,  (2)  a  high  ratio  of  electron 
teaperature  to  ion  and  gas  temperatures,  and  (3)  a  wider  range  .of  operating 
gas  pressures.  Furthermore,  sputtering  effeots  due  to  the  electrode  asterial 
are  eliminated  hy  the  inductively-coupled  configuration. 

A  major  improvement  in  diagnostics  is  provided  hy  the  addition  of  a 
quadrupole  mass  spectrometer.  This  system  features  the  ability  to  detect  both 
positive  and  negative  ions.  The  quadrupole  contains  its  own  source  of 
ionization  which  can  be  switched  off.  The  system  can  then  be  utilized  to 
directly  detect  species  in  the  plssaa  along  with  providing  means  to  analyze 
downstream  products.  Currently,  the  QMS  is  connected  downstream  of  the 
discharge-sample  region  via  a  variable-leak  valve.  A  further  improvement 
would  connect  the  quadrupole,  via  a  pinhole  leak,  to  the  discharge-sample 
region  and  would  enable  the  detection  of  many  of  the  ionic  species  involved  in 
the  etching  process. 

Much  of  the  work  currently  being  done  in  this  field  concerns  the 
"Reactive  Ion  Etching"  scheme.  This  differs  from  conventional  "plasma 
etching"  in  that  physical  processes  as  well  as  chemical  are  used  to  remove 
materiel.  Concern  has  been  raised  recently  that  damage  to  the  semiconductor 
occurs  due  to  the  bombarding  ions.  Surface  states  at  the  SiOj-Si  boundary, 
for  example,  are  thought  to  be  caused  by  this.  Machines  designed  to  exploit 
the  RIE  scheme  employ  a  parallel  plate  IF  discharge  whereby  the  electrode  upon 
which  the  samples  are  placed  assumes  a  negative  bias  proportional  to  the  ratio 
of  the  surface  area  of  the  electrodes.  The  system  configuration  we  have 
constructed  avoids  this  and  allows  a  DC  bias  to  be  superimposed  on  the  sample. 
Therefore,  we  hope  to  study  the  influence  of  ion  bombarding  effects  on  the 
anisotropy  of  etch  as  well  as  damage  to  the  semiconductor. 

Preliminary  etching  runs  have  been  conducted  and  etching  has  been 
achieved.  Further-  work  to  determine  etch  rate  versus  various  parameters  of 
the  system  are  being  conducted  in  an  effort  to  "calibrate"  the  system  with 
respect  to  known  results  in  the  literature.  In  particular,  attention  is 
focused  on  observing  etch  rate  versus  sample  bias  in  the  RF  system  as  compared 
to  the  results  already  obtained  in  the  DC  hollow-cathode  configuration. 
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3 .4  Laser-Induced  Growth  si  Semiconductor  Films* 

The  objective  of  this  resasroh  program  is  to  investigate  tho  ehasiitzy 
sad  optiesl  propartias  of  axoitsd  atoms  sad  aolaealsr  radicals  producad  by 
ultraviolet  photodissooiatioa  aad  tha  iataraetioa  of  thasa  axeitad  species 
with  a  surface  to  goaarato  a  thia  film.  Specifically,  aa  altraviolat  axeimar 
lasar  is  baiag  asad  to  photodissoeiata  iaorgaaie  molaculas  that  eoataia  a 
saaiooaductor  atom  (such  as  SiH^,  GaHj  or  SiClj) .  Tho  apeetrosoopy  of  those 
lasar-gaaaratad  plasmas  as  wall  as  tha  olaetrieal  aad  ohamieal  propartiaa  of 
tha  resultiag  Ga  or  Si  filaa  ara  wader  iaveatigatioa. 

Iaitial  axpariaaats  iavolved  tho  photodisaoeiatioa  of  SiH^,  (CHjJ^Si  aad 
in  the  altraviolat  us lag  oxoiaar  lasars  aad  aaoeaadad  ia  growiag  Si  aad 
Ga  films  oa  amorphous  SiOj  substrates.  Thasa  experimeats  ware  coaductsd  at 
room  tamperatura  aad  tha  resultiag  Ga  or  8i  films  wars  fouad  to  be 
polyerystalliae  with  graia  sixes  of  ip  to  0.5  |im.  However,  tha  films  ware 
highly  resistive  due  to  iatorgraaular  voids  formed  as  a  result  of  tha  large 
dapositioa  rates  (-  1  am  -  s"1) . 

Ia  more  raeaat  experimeats,  a  Cf  <X>2  laser  has  baaa  used  ia  conjunction 
with  tha  axeimar  lasar  to  produee  films  of  much  higher  quality.  Tha  C02  laser 
maiataias  tha  spatial  selectivity  of  tha  teohniqae  aad  serves  to  awderstely 
hdat  the  substrate,  thas  eahsaeiag  tha  surfaea  mobility  of  sa  absorbed 
semieoaductor  atom.  However,  siaee  it  is  the  exoimer  laser  that,  by  a  noa- 
tharmal  process  frees  the  Ga  or  Si  atom  from  the  pareat  molecule,  the  films 
caa  be  deposited  at  temperatures  at  least  100-150°C  lower  thaa  those  required 
for  pyrolytic  dapositioa.  For  example,  Ge  films  growa  at  a  substrata 
temperature  of  300°C  for  10-20  minutes  are  -  1  pm  thick  aad  exhibit  room 
temperature  mobilities  of  **  100  cm^-V^s”* .  Hall  measurements  also  indicate 
carrier  eoaeaatratioas  of  ~  10*®  cm"3  aad  x-ray  studies  using  a  powder 
diffractometer  show  the  films  to  be  polyerystalliae  with  no  preferred 
orieatatioa.  Substrate  temperature  measurements  show  that  the  axeimar  laser 
heats  the  substrate  by  10°C  or  less.  Also,  we  have  demonstrated  that  (111)  Ge 
caa  be  epitaxially  growa  oa  ll02  oriented  sapphire  substrates  by  pyrolyxiag 
GeH^,  Current  experiments  are  foousiag  oa  accomplishing  the  same  goal  at 

*  This  work  was  supported  by  the  Joiat  Services  Electronics  Progrsm  (U.S. 
Army,  U.S.  Navy,  U.S.  Air  Force)  under  contract  N00014-79-C-0424  aad  by  the 
Air  Force  Office  of  Scientific  Researeh  under  coatraot  AFOSR  79-0138. 
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lower  substrate  teapera tares  bp  photodissociating  GeH^ . 

Eaission  froa  atoaie  Ge  and  the  geraylene  (GeH) ,  bas  also  been  observed 
and  profiles  of  fluorescence  intensity  versns  distance  froa  tbe  substrate  are 
being  aeasured.  These  results  will  later  be  incorporated  into  a  ooaputer 
aodel  of  the  gas  phase/surfaee  kinetics  in' an  effort  to  better  understand  the 
processes  involved  in  the  nueleation  of  Ge  filas. 
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4.1  Introduction 

This  research  involves  study  of  basic  properties  of  semiconductors, 
methods  of  devioe  processing,  and  new  device  concepts.  Both  theoretical  and 
experimental  methods  are  employed  in  each  of  these  categories.  We  are 
examining  a  variety  of  hot  electron  phenomena  and  their  effects  on  present  and 
future  devioe  performance,  espeoially  in  connection  with  modulation  doping. 
Ve  are  studying  ion  implantation  and  annealing  of  Si  and  III-V  compounds, 
ineluding  laser  and  electron  beam  processing.  These  experimental  studies 
include  examination  of  deep-level  impurities  and  defects  arising  from 
implantation  and  annealing.  Several  aspects  of  this  work  are  done  in 
collaboration  with  other  units,  particularly  the  studies  of  materials  grovn  by 
molecular  beam  epitaxy  (MBS,  Unit  1). 
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4.2  Laser  £Bd  Electron  Been  Anneal in«  si  Seniconductors* 

In  our  studies  of  the  physios  and  Materials  problens  underlying  new 
device  processing  techniques*  we  have  exanined  several  effects  of  transient 
annealing  in  ssnieonduetors.  While  these  Methods  .exhibit  great  pronise  for 
More  flexible  processing  of  devices,  a  better  understanding  of  the  annealing 
nechaniSMS  is  required  for  reliable  application  to  fabrication.  In  the 
previous  reporting  period  we  described  extensive  experinents  in  pulsed  laser 
annealing  of  silicon.  During  this  period  we  have  published  results  of  that 
work  applied  to  BF^*  anorphized  Si.  We  have  also  exanined  electrical 
activation  and  residual  defeets  (using  DLTS)  in  anorphized  Si  recrystallized 
by  swept  line  electron  bean  annealing.  We  have  also  presented  results  of  a 
new  electron  bean  annealing  approach  in  which  electrons  are  derived  fro*  a 
plasna  discharge. 

Considerable  work  haa  been  done  on  recrystallization  of  deposited  silicon 
on  Si<>2  using  lateral  seeding  froa  windows  in  the  oxide  to  underlying 
crystalline  Si.  The  approach  used  ie  heating  by  a  narrow  carbon  strip  which 
traverses  the  sas^le  surface  in  a  swept  line  pattern.  Large  single  crystals 
have  been  obtained  by  this  Method*  but  these  crystals  contain  low-angle  grain 
boundaries.  Studies  are  underway  on  the  crystal lographio  and  electrical 
properties  of  these  filns. 

4.3  iMPlantation  and  Q*  CoMUQUud  Seniconductors** 

In  the  last  reporting  period  we  described  studies  of  encapsulation  and 
annealing  in  InP.  Based  upon  those  results*  we  have  carried  out  extensive 
neasurenents  of  inpurity  Migration  in  inplanted  InP.  We  find  that  iuplanted 
inpurities  (e.g.*  Be)  undergo  anoualous  Migration  during  annealing  of  the 
danage*  and  also  that  eonpensating  inpurities  (Fe*  Cr)  used  to  nake  seni- 
insulating  InP  Migrate  considerably.  Correlations  were  observed  between  the 


*-This  work  was  supported  by  the  Joint  Services  Electronics  Progran  (U.S. 
Amy,  U.S.  Navy,  and  U.S.  Air  Force)  under  Contract  N00014-79-C-0424*  by  the 
Army  Research  Office  under  Contract  DAAG  29-80-C-0011 ,  and  by  the  Office  of 
Naval  Eesearoh  under  Contract  N0001 4-81 -1-0431. 

•*  This  work  was  supported  by  the  Joint  Services  Electronics  Progran  (U.S. 
Amy,  U.S.  Navy,  and  U.S.  Air  Force)  under  Contract  N00014-79-C-0424,  and  by 
the  Office  of  Naval  Research  under  Contract  N00014-76-C-0806. 
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migration  of  implanted  Be  and  either  Fe  or  Cr.  After  annealing,  flat  taila  of 
Be  extending  over  several  microns  have  been  observed,  with  an  accompanying 
decrease  in  Fe  or  Cr  over  the  same  regions.  Similar  flat  tails  were  observed 
in  the  case  of  implanted  Mg.  The  migration  of  implanted  Si  in  InP  is  mnch 
less  pronounced;  however,  Fe  redistribution  is  observed  which  can  seriously 
affect  the  electrical  properties  of  Si  implanted  layers. 

Studies  of  anomalous  impur i ty  migration  have  also  been  initiated  in  GaAs. 
Preliminary  results  indicate  that  iaplanted  S  in  GaAs  diffuses  much  less 
during  annealing  when  Si  is  present.  However,  this  effect  is  pronounced  only 
at  relatively  low  S  doses.  At  higher  S  concentrations,  the  deep  migration 
coansonly  observed  is  not  ia^eded  by  Si  co-implantation.  Work  is  currently 
underway  with  GaAs  layers  doped  with  Si  during  growth,  to  separate  out  the 
effects  of  defects  and  chemical  effects  between  the  impurities. 

4.4  Electronic  Transport  H  Hiah  Energies* 

The  goal  of  this  research  has  been  to  develop  a  complete  theory  of 
transport  phenomena  at  high  energies  (>  1  eV  above  the  conduction  band  edge) 
in  bulk  semiconductors.  In  previous  years  we  have  developed  a  Monte  Carlo 
simulation  which  included  a  realistic  band  structure  as  calculated  by  the 
empirical  pseudopotential  method.  The  scattering  rates  were  essentially 
computed  by  the  Golden  Rule  with  the  deformation  potential  constants  regarded 
as  adjustable  parameters.  From  the  fit  to  a  large  number  of  experimental 
results  (impact  ionization  rate,  saturation  velocity,  effects  such  as  electron 
collection  in  double-heterojunctions  and  real  space  transfer)  we  have  obtained 
a  reliable  set  of  coupling  constants.  The  scattering  rates  are  proportional 
to  the  square  of  the  coupling  constants  and  the  final  density  of  states. 
Because  of  the  peaks  in  the  density  of  states  there  are  regions  in  the  bands 
exhibiting  extremely  high  scattering  rates.  These  high  rates  make  it 
desirable  to  use  a  more  rigorous  approach  than  the  use  of  the  Golden  Rule  to 
include,  e.g.,  collision  broadening  effects.  Therefore  we  calculated  the  self 
energy  I  from  the  relation  (see  e.g.  J.  M.  Ziman,  Elements  Advanced  Quantum 
Theory .  Cambridge  Univ.  Press  1969) . 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program  (U.S. 
Army,  U.S.  Navy  and  U.S.  Air  Force)  under  contraot  N00014-79-C-0424  and  by  the 
Office  of  Naval  Research  under  contraot  N00014-76-C-0806 . 
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£  -  li  +  i8  B(H  -  .  TS=H  +  r  444* 

from  E  the  scattering  orossection  is  obtained  via  the  optioal  theorem.  Fig.  1 
ehows  2  In  as  well  as  the  scattering  rate  as  calculated  by  the  Golden  Rule 
(G) .  To  lowest  order  the  two  curves  should  be  identical.  Notice  the 
deviation  between  the  curves  at  high  energies.  Te  have  aieanwhile  performed 
calculations  of  the  impact  ionisation  rate  including  these  effeots.  The 
qualitative  features  of  the  theory  of  impact  ionisation  remain  the  aame;  there 
are*  however*  quantitative  differences.  These  become  especially  important  at 
large  energies  and  will  influence  significantly  processes  such  as  the  emission 
of  electrons  from  silicon  into  silicon  dioxide  which  is  so  important  for 
stability  considerations  of  micron  and  submicron  devices.  Detailed 
assentations  are  in  progress. 

4.5  Transport  in  Layered  Hotoroiunctidn  Structures* 

Ve  have  investigated  the  possibility  of  storage  and  switching  between  the 
layers  of  GaAs  separated  by  Al^Ga^^As  layers.  The  test  structure  shown  in 
Fig.  2  was  fabricated  using  standard  photolithographic  and  ehemical  etching 
techniques.  Four  ohmic  contacts  were  formed  by  vacuum  evaporation  of  Au-Sn 
and  subsequent  alloying  at  450°  for  30  sec  in  flowing  R 2,  The  resulting 
structure  had  a  dumbbell-shaped  oenter  region  containing  all  nine 
heterostructure  periods*  each  of  which  was  contacted  at  each  end  by  an  ohmic 
contact  (labelled  side  contact  in  Fig.  2)  that  was  diffused  down  through  most 
ef  the  layers.  The  areas  containing  the  main  contacts  (see  Fig.  2)  were 
etched  down  so  that  the  top  three  GaAs-il^Qaj^As  layer  pairs  were  removed* 
leaving  only  the  lower  sir  periods  contacted  by  the  main  oontaots.  If 
voltages  are  applied  only  to  the  main  contacts  (no  applied  voltage  on  side 
contacts)  then  only  the  GaAs  layers  in  the  bottom  six  periods  will  have  high 
eleetrio  fields.  Some  of  the  electrons  that  escape  from  the  GaAs  layers  in 
these  bottom  six  periods  can  propagate  into  the  upper  electrically  isolated 
layers  nearer  the  surfaoe  in  the  dumbbell  bridge.  If  the  sample  is  cooled  to 
cryogenic  temperatures,  the  electrons  that  lose  energy  in  these  isolated 


*  This  work  was  supported  by  the  Office  of  Naval  Research  under  contract 
N00014-76-C-0806  and  the  Army  Research  Office  under  contraot  DAAG29-80-C-0011 . 
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Figure  4.1  Imaginary  and  real  (inset)  part  of  the  self  energy  as 

calculated  in  the  Fock  approximation  (full  line)  and  the 
Golden-rule  result  (dashed  line)  fcr  the  conduction  band 
of  GaAs.  For  weak  coupling  the  two  curves  should  be  identical 
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layers  will  become  trapped  in  the  GaAs  regions  and  will  not  be  able  to  gain 
enough  energy  from  the  crystal  lattice  or  fringing  low  electric  fields  to 
escape  by  thermionic  emission.  Thus,  as  soon  as  the  electric  field  between 
the  main  contacts  becomes  high  enough  to  cause  electrons  to  be  emitted  from 
the  lower  GaAs  layers*  a  fraction  of  these  electrons  can  propagate  to  the 
upper  isolated  GaAs  layers  where  they  are  trapped,  provided  that  the  side 
contacts  do  not  connect  through  to  the  lower  layers  interconnected  by  the  main 
contacts.  The  resulting  current  between  the  Min  contaots  will  be  reduced  for 
all  following  measurements  of  the  current  at  lower  fields.  Suppose  a  voltage 
is  then  applied  only  to  the  side  contacts  so  that  the  resulting  electric  field 
in  the  dumbbell  bridge  region  is  high  enough  to  allow  thermionic  emission  of 
the  electrons  in  all  the  upper  GaAs  layers  as  well  as  those  between  the  main 
contacts.  The  initially  higher  concentration  of  electrons  stored  in  the  top 
three  periods  will  then  be  redistributed  to  the  other  periods,  thus  restoring 
the  original  carrier  distribution  of  the  sample  before  any  voltages  were 
applied.  The  current-field  characteristics  between  the  main  oontacts  at  low 
fields  will  then  also  return  to  their  original  values. 

The  results  of  current-field  measurements  between  the  main  contacts  on  a 
sample  cooled  to  10  K  is  shown  in  Fig.  3.  Measurements  were  performed  using 
700  ns  voltage  pulses  at  low  repetition  rates  to  reduoe  sample  heating. 
Measurements  of  the  current  vs.  field  shown  in  the  graphs  were  taken  at  600  ns 
after  the  beginning  of  each  pulse,  although  no  time  dependence  was  observed  in 
the  pulses  between  1  and  700  ns.  The  sample  was  mounted  on  a  temperature 
controlled  cold  finger  in  an  evacuated  sample  chamber  and  standard  50-0 
sampling  oscilloscope  and  x-y  recorder  techniques  were  used.  All  measurements 
were  performed  with  the  sample  in  the  dark.  The  upper  solid  line  in  Fig.  3 
shows  the  virgin  enrrent-field  characteristics  between  the  main  oontacts  when 
the  first  voltage  was  initially  applied  to  the  main  contacts  (no  voltage 
between  side  contacts) .  After  the  eleotric  field  exceeded  a  threshold  value 
necessary  to  cause  a  real-space  transfer  of  electrons,  the  voltage  was 
returned  smoothly  to  zero,  resulting  in  the  lower  solid  curve. 

The  significantly  lower  current  values  measured  as  the  voltage  was 
returned  to  zero  indicates  that  a  fraction  of  the  electrons  have  been 
transferred  and  stored  and  can  no  longer  contribute  to  the  ourrent.  Repeating 
the  same  voltage  scan  several  seconds  later  (dashed  curve  in  Fig.  3)  revealed 
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only  a  slight  change  in  the  reduced  current.  Another  scan  of  the  voltage 
several  ainutes  later  (dot-dashed  curve)  exhibited  only  a  saall  further 
increase  in  the  reduced  current,  indicating  the  number  of  stored  electrons  had 
changed  very  little.  Immediately  thereafter,  a  high  field  vas  applied  between 
the  side  contacts  (no  applied  voltage  on  main  contacts)  to  redistribute  the 
stored  electrons.  That  following  voltage  scan  between  the  main  contacts 
resulted  in  a  retracing  of  the  initial  virgin  solid  line  current-field  curve. 
The  above  pattern  of  voltage  applications  was  repeated  many  times  with  the 
same  results.  The  experiments  are  still-  inconclusive  as'  to  whether  the 
electrons  are  stored  in  the  upper  GaAs  layers  or  in  deep  traps  in  the  AlxGsj_ 
jAs  or  in  surface  states,  since  a  complete  electrical  separation  of  the  layers 
(side  contacts)  is  difficult  to  achieve.  Also,  the  density  of  deep  traps  is 
not  well-known  in  these  samples  and  more  research  needs  to  be  done.  However, 
the  effect  not  only  proves  the  concepts  of  real  space  transfer  as  discussed 
before,  but  also  shows  the  device  potential.  In  addition,  intereating 
physical  parameters  may  be  determined  by  time  resolved  measurements  of  the 
kind  described  above,  sinoe  all  these  effects  are  extremely  sensitive  to 
electron-phonon  interaction,  quantum  mechanical  transmission  coefficients, 
etc. 

4.6  Xmwms.*  lid  TUHtift  Safari  SA9n  fisxiaai* 

Ve  have  used  our  Monte  Carlo  simulation  program  to  investigate  the 
possibility  of  ultrafast  collision  free  transport  as  proposed  by  Eastman  at 
Cornell.  The  result  of  our  simulation  is  that  in  a  narrow  range  of  parameters 
ultrahigh  speeds  are  possible.  Collectors  and  drain  regions  of  current 
semiconductor  devices,  however,  are  outside  this  parameter  range  (because  of 
the  high  voltage  drops)  and  therefore  represent  the  bottleneck  for  speed. 
This  problem  is  fundamental  in  nature  and  can  only  be  circumvented  by 
constructing  devices  with  extremely  short  collector  or  drain  regions.  Fig.  4 
shows  typical  results  of  transit  times  vs  device  length  for  GaAs  with  the 
electric  field  as  parameter.  The  electrons  are  assumed  to  be  injected  at  300 
meV  energy  and  transverse  the  distance  indicated  at  the  x-axis.  Notice  the 

paradoxical  increase  in  transit  time  with  increasing  field,  which  is  csused  by 
the  population  of  higher  conduction  band  minima. 

•  This  work  was  supported  by  the  Joint  Services  Electronics  Program  (U.S, 
Army,  U.S.  Navy,  and  U.S.  Air  Force)  under  contract  N00014-79-C-0424  and  the 
Army  Research  Office  under  contract  DAAG29-80-E-0069 . 
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AM  Later 


Ve  are  studying  energetic  particle-surface  interactions  and  stiaulated 
gas  phase  reactions  which  control  the  nneleation  and  growth  kinetics, 
eheaistry,  and  physical  properties  of  ooaponnd  and  alloy  seaieondnotors  grown 
froa  the  vapor  phase  by  UHV  ion  beaa  sputtering,  plasaa-assisted  cheaical 
vapor  deposition,  accelerated  beaa  MBE,  and  laser-induced  cheaical  vapor 
deposition.  The  eoaaon  feature  of  these  techniques  is  that  crystal  growth 
proceeds  under  non-equilibriua  theraodynaaic  conditions  through  the  production 
of  highly  reactive  gas  phase  speeies:  excited  atoas,  aetastables,  radicals, 
and  ions.  Such  species  transfer  energy  to  the  growth  surface  upon 
condensation  thereby  altering  the  surface  reactivity  as  well  as  adsorption  and 
adatom  diffusion  kinetics  allowing  fila  growth  at  lower  teaperatures,  a  wider 
range  in  controlling  doping  concentrations  and  tailoring  fila  properties,  and 
the  growth  of  unique  aetastable  aaterials.  This  work  is  being  pursued  from 
both  an  analytical  and  an  experimental  point  of  view  to  establish  a  detailed 
understanding  of  fuudaaental  fila  growth  aeohanisas.  Ve  have  recently 
published  two  invited  review  papers  in  this  area  [5.3,  5.4].  Besults  froa 


*  This  work  was  supported  by  the  Joint  Servioes  Electronics  Prograa  (U.S. 
Any,  U.S.  Navy  and  U.S.  Air  Force)  under  contraot  N00014-79-C-0424 . 
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this  research  have  s  wide  range  of  applications  in  addition  to  crystal  growth 
inelnding  tha  active  rasaareh  fialda  of  rase t Its  ion  etching,  ion  beam 
lithography,  nicroehaaical  analyais,  plasma  chemistry,  and  laaar-mataxials 
intaraetions. 

5 .2  Ion/Snrfaca  Intaraetions 

Va  have  developed  general  modala  for  predicting  alaaantal  incorporation 
probabilities  of  aeealaratad  spaoias  incident  upon  growing  film  surfaces 
during  vapor  phase  deposition  [5.4].  Trapping  (low  energy  implantation), 
scattering,  re-evaporation,  and  preferential  sputtering  are  accounted  for. 
More  recently  we  have  added  terms  for  eollisional  mixing  and  segregation  of 
dopants.  As  discussed  briefly  below,  the  use  of  these  models  in  analyzing 
experimental  data  not  only  leads  to  a  more  detailed  understanding  of  film 
growth  processes  but  also  gives  predictive  capability  in  designing  new 
experiments  or  considering  the  growth  of  new  materials. 

5.2.1  BttflJLMtlgS  2S&  Growth  KUtllftl 

Initial  experiments  using  both  in-situ  electron  diffraction  and 
replication  transmission  electron  microscopy  have  been  carried  out  to 
investigate  the  effect  of  low  energy  (<  200  eV)  ion  bombardment  on  nucleation 
kinetics.  Preliminary  results  for  Go  and  GaAs  on  amorphous  SiO^  substrates 
show  more  than  an  order  of  magnitude  increase  in  grain  size  as  well  as  a  large 
inerease  in  the  degree  of  preferred  orientation,  (100)  and  (110) , 
respectively. 

Ion/surfaee  interaction  effects  have  been  found  to  greatly  affect  the 
growth  kinetiea  of  sputter  deposited  single  crystal  GaAs.  Experiments  were 
carried  out  on  (100)  GaAs  substrates  as  a  function  of  deposition  rate  R, 
growth  temperature  T#,  incident  As/Gc  flux  ratio  ^As^Ga*  *ad  *on  *eo*l***tion 
voltage  V  .  Atomic  Ga  and  As  beams  were  obtained  by  sputtering  from  a  single 
crystal  undoped  GaAs  wafer  while  the  simultaneous  evaporstion  of  arsenic  from 
an  effuaion  cell  provided  an  As^  overpressure.  The  results  show  that  in 
direct  contradistinction  to  MBS  growth  using  thermal  beams,  R  depends  not  just 
ob  Jqs  but  also  on  the  surface  composition  of  the  growing  film  and,  hence,  on 
Tg  and  the  relative  Ae4  overpressure.  This  is  primarily  due  to  secondary 
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spattering  iffteti.  An  analysis  of  tho  data  baaad  upon  earlier  stadias 
involving  the  rsaotion  af  thermal  aolsoalar  As^  Aad  0a  heated  OaAs 
surfaces  allows  a  determination  of  the  snrfaos  binding  snsrgy  of  6a  as  a 
f ana t ion  of  tbs  As  snrfaos  eororago. 

3.2.2  gitutil  Xiagafliiiifla  PirtrtlUtln 

Aebiowablo  doping  concentrations  of  naay  important  dopants  in  MBE  f ilas 
grown  mador  praotieal  conditions  has  boon  1  ini  tod  by  low  tbomodynaaio 
sticking  probabilities.  However,  accoloration  of  these  dopants  to  relatiwoly 
low  energies,  a  fow  hundred  #?,  has  boon  shown  by  soworal  investigators  to 
result  in  eohameeawats  in  not  incorporation  probabilities  in  6aAs  and  Si  of  up 
to  eight  orders  of  nagnitudo.  Vo  reported  sinilar  effects  for  S  doping  of 
sputter  deposited  OaAs  last  year. 

Using  the  aodel  deseribed  briefly  in  section  5.2  above,  we  have  analysed 
our  experimental  results  for  aeoelerated  8  doping  of  sputter  deposited  GaAs 
15.31  and  As  doping  of  MBS  Si  (grows  by  Dr.  Ota.  Bell  Laboratory,  Boladale). 
In  both  oasos,  several  orders  of  nagnitudo  increases  were  observed  in  the 
overall  incorporation  probability,  however  the  detailed  nechanisas  involved 
were  quite  different,  in  the  former  ease.  S  concentrations  (S]  of  up  to  1.2  x 
1021  cm3,  as  aeasored  by  secondary  ion  mass  spectrometry,  were  obtained  using 

•*  HjS/Ar  glow  discharge  source.  The  primary  mechanisms  of  incorporation 
involved  a  competition  between  dissociative  chemisorption  coupled  with 
trapping  and  preferential  secondary  sputtering  giving  rise  to  a  minimum  in  [S] 
at  a  critical  substrate  potential.  [S]  increased  linearly  with  increasing' 
partial  pressure  but  was  not  affected  by  varying  the  As/Ga  flax  ratio  between 
12  and  40.  Very  sharp  doping  profiles  (see  also  section  5.2.3)  with  no 
indication  of  surface  segregation  were  obtained  for  progremmed  doping  steps. 
In  the  Si:As  ease,  the  total  ion  current  was  low  enough  to  neglect  sputtering, 
but  enhanced  near-surface  diffusion  and  segregation  as  well  as  trapping  had  to 
be  considered. 

Ihile  the  growth  parameters  vers  very  different  in  the  above  experiments, 
the  calculated  results  agreed  very  well  with  measured  results  in  both  eases. 
The  model  is  now  being  used  to  make  predictions  in  other  systems  and  to 
provide  a  general  basis  for  understanding  accelerated  doping. 
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s.2.3  SmmUga  iM  fissiii 

A  segregation  aodel  was  developed  for  predicting  surface  accnanlation  of 
dopants  and  the  corresponding  dopant  depletion  region  near  the  film/ substrate 
interface  daring  MBE  and  sputter  deposition.  The  aodel  does  not  require  a 
rate— kinetic  liaitation  to  lattiee  incorporation  but  inatead  is  based  on  a 
surface  segregation  potential  controlled  by  the  strain  free  energy  resulting 
froa  aisfit  solute  atons,  "bond  breaking"  terns  due  to  a  decrease  in  the 
surface  free  energy  when  bulk  species  with  lower  enthalpies  of  snbliaation  are 
exchanged  with  surface  a ton  species,  and  the  excess  entropy  of  aixing.  The 
depth  distribution  of  dopant  throughout  the  entire  fila  can  be  calculated  by 
solving  the  full  position  and  tiae  dependent  second  law  of  diffusion 
accounting  for  both  diffnsional  and  drift  flnxea.  The  segregation  potential 
enters  into  the  latter  tera  through  the  drift  velocity.  The  equations  are 
solved  in  differential  difference  fora  subject  to  boundary  conditions 
depending  on  the  nature  of  the  experiaent. 

As  a  firat  application  of  the  aodel  [5.10],  we  applied  it  to  the  well 
known  problea  of  Sn  aurface  accnanlation  in  MBE  grown  GaAs  where  the  anrface 
concentration  0gn  i*  typically  3  to  4  orders  of  aagnitude  greater  than  that  of 
the  bnlk  and  a  Sn  deficient  region  up  to  100  na  wide  is  left  at  the 
fila/ substrate  interface  under  noraal  growth  conditions.  Calculated  Sn 
profiles  were  found  to  agree  extreaely  well  with  aeasured  profiles,  obtained 
uaing  secondary  ion  aass  speetroaetry  (SIMS)  in  several  series  of  GaAs:Sn 
saaples  (grown  in  collaboration  with  Prof.  Morkoc)  where  the  growth 
teaperature,  the  As^/Ga  flux  ratio,  and  the  Sn  flux  was  varied.  By  adding 
teras  accounting  for  trapping  and  collisional  aixing  (recoil  ^plantation)  we 
can  also  aodel  accelerated  beaa  doping  experiaents  and  are  now  applying  the 
aodel  to  interpret  and  explain  published  results  for  accelerated  doping  of 
MBE-grown  Si.  Finally,  the  aodel  can  also  b«  used  to  predict  experiaental 
conditions  leading  to  better  control  over  dopant  distributions.  Using  these 
ideas,  we  have  deaonstrated  a  reduction  of  9-  in  GaAs  by  -  two  orders  of 
aagnitude  with  a  correspondingly  sharper  profile  at  the  fila/ substrate 
interface  for  the  sane  Cg&  used  in  the  above  experiaents. 
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5.2.4  Growth  si  Simla  Crystal  M.ejgitafrlp  g9li99f4B9*9» 

As  reported  lest  year,  we  ere  carrying  out  the  first  detailed  study  of 
the  growth  and  physical  properties  of  new  single  crystal  awtastable 
semiconductors.  The  hey  feature  in  stabilising  the  growth  of  these  materials 
is  the  controlled  use  of  low  energy  ion  bombardment  during  deposition  to 
modify  elemental  sticking  probabilities  and  adatom  diffusivities  and  to 
promote  collisional  mixing.  The  research  described  in  sections  5.1  through 
5.2.3  above  has  provided  a  much  better  understanding  of  the  growth  of  these 
unique  materials.  During  the  past  year  we  have  concentrated  our  efforts  on 
the  study  of  (GsAs)^_z6ex  because  of  the  importance  of  the  end-members,  the 
interest  in  Ge/GaAs  heterostructures,  and  the  fact  that  it  is  representative 
of  a  new  subclass  of  potentislly  important  alloys,  (III-V) j_zIVz. 

Epitaxial  metastable  (GaAs)^_zGex  alloys  with  compositions  ranging  from  x 
*  0  to  x  *  1  have  been  grown  on  (100)  GaAs  substrates  by  ion  beam  sputtering 
in  an  ultrahigh  vacuum  system.  Electron  channeling,  double  crystal  X-ray 
diffrae  tome  try,  and  X-ray  t^graphy  analyses  indicate  that  the  films  are  of 
very  high  crystalline  perfection.  Either  n-type  or  p-type  conduction,  with  n 
and  p  varying  over  several  orders  of  magnitude,  could  be  obtained  by  varying 
the  film  composition,  the  growth  temperature,  and  the  As  overpressure  during 
deposition.  The  equilibrium  GsAs-Ge  pseudobinary  phase  diagrsm  has  been 
determined  by  differential  thermal  analysis  to  be  a  simple  eutectic  with  an 
invariant  temperature  and  composition  of  880°C  and  ~  18  mole  %  GaAs.  A  great 
deal  of  insight  into  the  bonding  structure  of  these  alloys  hss  been  achieved 
by  fitting  calculated  liquidus,  solidus,  snd  soluns  boundaries,  based  on 
pairwise  potential  quasi-chemical  models,  to  the  experimental  results.  Such 
models  inherently  contain  not  only  free  energy  (i.e.  activity)  vs  x  data  but 
also  require  knowledge  of  nesrest-neighbor  bonding  configurations. 

Collaborations  established  with  Prof.  Dow  (theory).  Prof.  Raceah  at  UICC 
(Raman),  and  Dr.  Aspnes  at  Bell  Laboratory,  Murray  Hill  (spectroscopic 
ellipsometry)  have  also  provided  a  more  quantitative  understanding  of  these 
alloys.  Raman  data  show  single  mode  behavior  across  the  alloy  system  with  the 
LO  mode  obeying  a  Yegard's  law  type  of  behavior  while  the  TO  mode  frequency 
increased  approximately  linearly  up  to  x  ■  0.75  and  then  much  more  rapidly 
thereafter.  The  elemental  mass  differences  were  too  smell  to  give  rise  to 
localised  modes  but  strong  alloy  broadening  was  observed.  The  broadening 
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effects  correlate  with  disorder  clone  the  (111)  bonding  direction.  Optical 
absorption  aeasurements  indicated  a  large  negative  bowing  in  the  direct  EQ 
band  gap  as  a  function  of  alloy  composition.  However  the  shape  of  the  bowing 
can  not  be  fitted  by  staple  dielectric  parabolic  aodels.  Ellipsoaetric  data 
also  indicate  large  negative  bowing  in  the  direct  and  Ej  transitions. 

Very  recent  aeasnreaents  nsing  analytical  scanning  transaission  electron 
aicroscopy  coabined  with  energy  dispersive  X-ray  spectroscopy  indicate  that 
the  alloys  are  hoaogeneous  to  within  the  3  to  5  na  resolution  of  the 
techniques  while  phase  separation  was  observed  after  annealing  near  800°C. 
This  technique  offers  potential  for  in-situ  phase  transforaation  studies  as 
well  as  for  resolving  current  arguments  concerning  clustering  in 
superlattices. 

5.2.5  Reactive  Ion  Etching  of  GaAs 

The  first  detailed  study  of  the  reactive  ion  etching  (RIE)  of  GaAs  has 
been  carried  out  by  our  group  [5.8].  Ve  investigated  the  RIE  of  (100)  GaAs  in 
pure  CCljFj  and  CF^  discharges  as  well  as  in  mixtures  of  Ar  and  CC14,  CCljFj. 
or  CFj.  Anisotropic  etching  with  removal  rates.  R,  of  up  to  800  nm/ain  have 
been  obtained  in  reactive  discharges  operated  at  a  pressure  of  40  aTorr  and  a 
target  voltage  of  -3kV.  The  physical  sputtering  rate  in  pure  Ar  discharges 
operated  under  the  saae  conditions  was  only  40  nm/ain.  A  combination  of  in- 
situ  optical  eaission  and  absorption  spectroscopies  have  been  used  to  show 
that  in  both  pure  and  dilute  (up  to  90  aole  %  Ar)  halocarbon  discharges, 
physical  sputtering  of  atomic  Ga  and  As*  is  not  a  primary  etching  mechanisa  for 
GaAs,  although  ion  bonbardaent  does  plsy  a  critical  role  in  the  overall 
process.  Transient  glow  disehsrge  optical  spectroscopy  aeasureaents 
deaonstrated  that  while  F.  increased  with  increasing  Cl/F  ratios  in  the  etch 
gas,  the  steady  state  carbon  concentration  at  the  GaAs  surface  also  increased 
indicating  that  carbon  accumulation  is  not  the  rate  liaiting  step  to  etching. 
Rather,  the  rate  limitation  is  provided  by  the  desorption  kinetics  of  gallium 
halides  which  we  believe  are  ejected  primarily  (except  in  pure  CCl^)  through 
ion-assisted  processes  as  the  reduced  radicals  GaF^  and/or  GsClz  (x-1  or  2) . 
We  have  proposed  a  phenomenological  model  to  provide  a  qualitative  description 
of  the  etching  behavior  of  GaAs  in  mixed  hslocarbon/ inert  gas  discharges. 
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5.3  klfl  Stimulated  Reactions 

In  collaboration  with  the  Solid  State  Phyeies  Division  of  Eastman  Kodak 
Research  Laboratories,  epitaxial  regrowth  of  Ge/GaAs  heterostrnetnres  by 
scanned  laser  annealing  of  amorphous  Ge  films  on  GaAs  substrates  has*  been 
studied  as  a  function  of  laser  power  and  scan  rate  [5.6].  At  least  eight 
regimes  representing  different  film  regrowth  characteristics  were  observed. 
Of  these,  two  were  of  primary  interest.  At  low  powers  (between  ~  1.6  and  3.2 
W  for  a  beam  diameter  of  ~  40  pm)  and  scan  rates  between  1  and  400  cm/sec, 
polycrystalline  Ge  with  a  (100)  preferred  orientation  was  formed  by  an 
"explosive"  crystallization  mechanism.  At  higher  powers,  and  over  a  scan  rate 
range  of  25  to  400  cm/ sec •  single  crystal  metastable  (GaAs)j_xGex  alloys  were 
obtained  by  liquid  phase  regrowth,  typical  fil*  resistivities,  p,  were  as 
follows:  as-deposited  p  ■  1800-cm;  polycrystalline  films,  p  ■  3  x  10~^  O-cm; 
single  crystal  films,  p  »  9  x  10“*  0-cm. 

5.4  Lu££  g.UWLUted  Chemical  Vapor  Deposition* 

In  collaboration  with  Prof.  Eden  of  the  Electrical  Engineering  Dept,  at 
UIUC,  we  have  grown  the  first  polycrystalline  Si  and  Ge  films  [5.1]  on 
amorphous  SiO^  substrates  (average  substrate  temperatures  <  120°c)  by  the 
photodissociation  of  Si^/ty  or  GeHq/He  mixtures,  respectively,  using  pulsed 
ArF  (193  nm)  and  KrF  (248  nm)  exeimer  lasers.  For  both  Si  and  Ge,  the  film 
growth  rate  exhibited  a  strong  dependence  on  laser  wavelength  and  beam 
intensity  I  where  1  <  I  £  10  MW.  -  cm~2«  As-deposited  films  had  average  grain 
sizes  of  up  to  0.5  pm  and  the  grains  were  equiaxed  with  a  random  orientation. 
Ge  films  grown  at  slightly  elevated  temperatures,  ~  300°C,  were  p-type  with  a 
carrier  concentration  of  2  x  10*^  cm“^  and  a  corresponding  mobility  of  120 
cm^/V-sec.  Ge  films  doped  with  -  10^®  cm“^  A1  w:Jte  obtained  by  the 
simultaneous  photodissooiation  of  A1(CH3)3  aa<j  OeH^.  Optical  absorption  and 
emission  experiments  are  being  carried  out  to  investigate  photon-molecule  gas 
phase  reactions. 


*  This  work  was  supported  by  the  Office  of  Naval  Research  under  contract 
N00014-81-K-0 56 8 . 
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6.  taatOVAVE  ACOUSTICS* 
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6.1  Objective 

OaAs  is  •  unique  material  in  that  it  kis  semiconductor  propsrtiss  tilt 
asks  it  ideal  for  high-speed  sleetronio  deviees,  and  it  has  piaaoalaetrie 
propartias  tkat  make  it  nsakla  for  aeonstie  deviees.  Vary  powerful  monolithic 
signal  processors  will  be  feasible  on  OaAs  because  tko  electronic  deviees 
provide  tke  gain  and  control  while  the  aeoustie  devices  provide  the  very 
eoapact  low  power  fast  nesMry  required  for  signal  processing. 

When  a  surface  scoustic  wave  propagates  along  tke  aetallized  surface  of  a 
OaAs  substrate  it  produces  potential  wells  buried  a  few  aierons  or  so  under 
tke  surface  tkat  nova  along  with  tke  wave  at  the  velocity  of  sound.  These 
traveling  wave  potential  wells,  assoeiatcd  with  tke  aeoustie  wave,  are  used  to 
capture  tke  eleotrons  and  transport  thca  in  synekronisa  with  tke  wave.  Ike 
bbjeetive  of  this  research  is  to  investigate  electronic  transport  by  tke 
traveling  wave  potential  walls  for  various  configurations  of* multi- layer  OaAs 
substrates. 


6.2  ttiiMip tin  ia4 

Ike  principle  of  electronic  transport  by  surface  acoustic  waves  on  OaAs 
substrates  kas  bean  daaonstratad.  Tke  first  device  deaoastrated  this 
characteristic  as  a  delay  line  which  is  in  away  ways  siaiilar  to  a  conventional 
CCD.  Tke  ekarga  is  injected  at  one  point,  propagated  along  by  potential 
wells,  and  then  eolleeted  at  a  later  point.  Tke  tine  during  which  it  reaains 

•  This  work  was  supported  by  tke  Joint  Services  Electronics  Program  (U.S. 
Army,  U.S.  Navy,  and  U.S.  Air  Force)  under  contract  N00014-79-C-0424 . 
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is  tbs  substrate  represents  signal  delay.  The  difference  between  the  aeonstic 
char go  transport  (ACT)  devices  and  the  conventional  CCD  is  that  the  ACT 
devices  have  no  gate  strnctnres,  and  the  charge  is  carried  strictly  by  the 
potential  wells  of  the  acoustic  wave.  The  following  paragraphs  briefly 
describe  the  first  experiaents  and  the  results  of  those  experiaents. 

A  photograph  of  the  device  shows  the  overall  layout  of  the  ACT  and  the 
physical  geoaetry  of  the  aultilayer  buried  channel  structure  is  described  as 
follows:  The  layers*  grown  by  aolecular  beaa  epitaxy  on  a  [100]  cut  Cr  doped 
GaAs  substrate,  consist  of  a  1.7  pa  thick  P  layer  (Na  -  10**  ea~3)  and  a  4.7 
|ia  N  layer  (Nd  ~  10**  ca-*).  The  device  transfer  channel  located  under  the 
Sehottky  Barrier  region  is  oriented  in  the  <110>  direction,  and  is  delineated 
by  a  preferential  aesa  etch  to  a  depth  of  5  pa.  Obaic  oontaets  to  the  input 
and  output  diodes  and  to  the  exposed  P  layer  are  foraed  using  the  liftoff 
process.  After  the  contacts  are  alloyed,  the  A1  Sehottky  barrier,  acoustic 
interdigital  transducer,  and  interconnect  pads  are  foraed  using  standard 
photolithographic  techniques.  The  interdigital  transducer,  which  has 
approxiaately  100  electrode  pairs,  generates  7.9  pa  wavelength  SAW  over  a  beaa 
width  of  ISO  wavelengths  at  a  frequency  of  367.1  MHz.  The  transport  channel 
is  1  aa  wide  and  1.35  aa  long  corresponding  to  an  acoustic  delay  between  the 
input  and  output  diodes  of  0.S  ps. 

The  Sehottky  barrier  region  > j  depleted  with  DC  bias  and  the  SAT 
generated  by  the  transducer  propagates  froa  the  input  diode  to  the  output 
diode.  Charge  carriers  are  injected  into  the  wave  by  applying  a  short 
negative  pulse  to  the  input  diode.  The  electrons  are  bunched  and  transported 
along  the  ohannel  by  the  potential  wells  of  the  SAT.  Then  the  wave  aaplitude 
is  large,  the  injected  electrons  are  quickly  bunched  and  aoved  synchronously 
with  the  wave,  at  precisely  the  sound  velocity.  Upon  reaching  the  output 
diode  the  delayed  charge  packets  are  swept  out  of  the  wave  by  the  applied 
reverse  bias  resulting  in  a  current  spike  in  the  output  detection  circuit. 

In  the  absence  of  loading  and  diffusion  effects,  the  aininua  traveling 
wave  elaetrie  field  necessary  for  the  synchronous  transport  is  about  40V/ca. 
In  our  experiaents  the  aeonstic  power  is  10  aT  resulting  in  a  peak  electric 
field  of  about  600  T/aa.  The  fields  will  be  aa  order  of  aagnitude  greater  in 
future  experiaents.  The  output  diode  current  obtained  when  the  SAT  transducer 
is  driven  by  a  967  Ms  source  is  found  to  be  a  nearly  perfect  delayed  replica 
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of  tk«  input  vmfon. 

There  arc  several  fundamental  features  of  the  ACT  device: 

(1)  Synchronous  electron  transport  is  induced  using  the  SAV  as  a  natural 
“clock**  eliainatinf  the  need  for  transport  gates  and  the  multiphase*  clock 
drivers.' 

(2)  The  interfering  "clock"  noise  is  eliminated. 

(3)  The*  ACT  is  expected  to  be  monolithically  integratable  with  high  speed 
electronic  components. 

(4)  High  clock  frequencies  and  large  time  bandwidth  products  appear  to  he 
feasible  because  the  charge  transport  efficiency  is  high. 

The  experimental  results  show  that  the  ACT  has  the  potential  for  the 
implementation  of  high  speed,  large  '  time-bandwidth  product,  large  dynamic 
range,  monolithically  integrated  signal  processors. 

6 .3  lecent  Developments  ip  LAW  Techno 1 oar* 

6.3.1  Kfiftlfftt  IM  qglT9lYf.il 

Line  Acoustic  Wave  (LAW)  convolvers  with  a  measured  internal  efficiency 
factor  of  -40dbm  have  been  demonstrated.  At  least  lOdb  of  the  losses  are  due 
to  bulk  wave  generation,  and  resistance  in  the  conductors.  If  these  loss 
mechanisms  arc  reduced,  a  ~30dbm  Internal  efficiency  factor  should  be 
possible.  This  is  about  two  orders  of  suignitude  more  efficient  than  the 
Surface  Acoustic  Wave  convolvers  presently  in  use. 

The  dynamic  range  of  LAW  devices  surprisingly  high.  Devices 

teplemented  with  narrow  bandwidth  IDT  trans  vlaeed  along  the  edge  have 

been  operated  with  input  power  levels  in  excess  of  +  27 dim.  The  LAI 
convolvers  are  naturally  resistant  to  second  harmonic  generation  and  no 
saturation  affects  could  be  measured  st  these  power  levels.  Bowever,  the  wide 
band  LAW  IDTs  are  destroyed  while  operating  at  such  high  powers.  The  failure 
mechanisms  invariably  involve  the  high  electric  fields  sssoeiated  with  the 
wideband  IDT,  and  these  failures  will  be  reduced  with  the  Multi-Strip  Mode 

•  This  work  was  supported  by  the  U.S.  Air  Force  under  contract  F19628-81- 
1-003 1 . 
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Cetmtir  transducer.  (MSMC) 


<•3*2  1 IdjfcUd  K8HC  m  T»»4i9fU 

Broadband  MSMC  LAV  transducers  ara  iapleaented  by  generating  surface 
aeoaatie  waves  aad  then  convertiag  the*  to  LAVs.  The  efficiency  of  the 
surface  aeoustie  wave  to  LAV  node  conversion  proeess  is  expected  to  reach  the 
1  to  2  db  level.  It  should  then  be  possible  to  build  a  MSMC  LAV  transducer 
with  an  insertion  loss  of  about  Sdb.  This  will  be  further  reduced  by  two  to 
three  db  if  unidirectional  SAV  transducers  are  used  as  part  of  the  MSMC.  The 
bandwidth  of  the  MSMC  LAV  transducer  has  been  increased  froa  2%  last  year  to 
the  present  value  of  7%  and  devices  with  20%  to  30%  bandwidth  are  presently 
being  processed. 
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7 .1  Dlffnaloa  sad  Reactivity  ea  Metals* 

Part  of  oar  affort  baa  baaa  ooaeaatratad  oa  aaderstaadlag  how  polys toale 
aolaaalas  deeosppose  oa  a  aolid;  ia  thia  work  we  haws  baaa  fooaslag  npon  tka 
iataraatioaa  of  aethaae  with  taagstoa  sarfaeaa.  Tbaaa  iatsraotioas  aro 
partionlarly  iatsrsstiag.  Deeoaposltion  of  aethaae  oa  as tala  rsqalrsa  a 
siaoabla  aetlwatioa  energy  aad  thia  systaa  thorafora  afforda  a  good 
opportaaity  for  exploring  ths  aaohaaisa  of  saoh  raaetioaa.  Ia  tha  past,  oar 
stadias  haws  baaa  dons  oa  awaporatad  filas;  for  thass  ws  haws  bailt  up  a 
eoharaat  wiaw  of  tha  aolsoalar  awaats  oa  tha  aarfaea  [U.7.8].  Kaeaat  work 
haa  straasad  tha  rola  of  sarfaos  straetara  ia  tha  deooaposltiom  proaass.  For 
this,  aav  tsahaiqaaa  aad  sqaipaeat  haws  baaa  iapleaeated  to  allov  aaaaaraaaata 
oa  well-eharseterised  siagla  crystal  plaaas. 

Ia  ordar  to  follow  tha  iatsraotioas  of  aothaas.  a  rathar  aareaetlve  gas* 
wa  haws  adopted  eoataot  potential  aoasaroaiaats.  Those  aaka  it  possible  to 
establish  tha  aarfaea  sewerage  vithoat  pertarbiag  the  gas  with  aaargatie 
aleetroas*  high  fields*  or  theraioaie  soareas.  Tha  systesi  whieh  has  baaa 
bailt  for  eoataot  potential  aeasar sweats  oa  single  crystal  plaaas  is  shewn 
sehewatieally  ia  Fig.  7.1.  la  it  a  no  1 scalar  baaa  of  hot  gas  iapiages  oa  a 
single  crystal  speeiaea.  Changes  ia  tha  eoataot  potential  arising  whoa  tha 
impinging  aethaae  daeeaposes  oa  tha  aarfaea  ora  aeasar  ad  by  a  vibrating 

•  This  work  was  sapported  by  tha  Natioaal  toieaea  Foaadatioa  aadar  Grants 
DML  71-2501$  aad  S2-01SS4. 


Fig.  7.1  Top  view  of  ultrahlgh  vacuum  system  for  conCact  potential 

measurements  of  activated  chemisorption  of  methane  on  single 


crystal  sample. 
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zifirtu*  electrode  of  fold*  which  it  placed  in  eloae  proximity  to  the  surface 
before  and  after  gat  boabardaeat.  The  cleanliness  of  tbe  aaaple  prior  to 
aeasureaents,  as  well  as  tbe  contsaination  rate  in  tbe  systea,  is  deterained 
by  an  Anger  systea  relying  on  a  saall  retarding  analyser  bnilt  into  tbe 
reaction  ebaaber. 

The  results  of  studies  on  tbe  (100)  and  (211)  planes  of  tungsten  are 
shown  in  Fig.  7.2.  The  activation  energies  differ  significantly  on  tbe  two 
surfaces,  tbe  (100)  being  aore  unreactive.  There  is  also  a  sixeable  isotope 
effect  for  botb.  with  tbe  barrier  for  CD^  exceeding  that  of  CH^  by  aore  than 
20%.  In  these  experiaents  only  tbe  beaa  teaperature  ia  changed;  tbe  crystal 
is  aaintained  at  roca  teaperature.  These  activation  energies  therefore  refer 
entirely  to  excitations  Imparted  to  ineoaing  gas  aoleeules  in  order  for  thea 
to  react. 

Further  experiaents*  on  differently  oriented  planes  of  tungsten*  are 
underway  in  order  to  aore  fully  expose  tbe  relation  between  surface  structure 
and  reactivity  toward  aetbane.  One  thing  has  definitely  been  established  - 
tbe  atoaic  arrangeaent  of  tbe  surface  plays  a  very  significant  role  in 
affecting  tbe  activation  energy  for  aoleeular  decomposition. 

In  a  cbeaical  reaction  on  a  solid*  tbe  initial  act  of  molecular 
deeoaposition  is  only  one  of  a  sequence  of  steps  which  finally  leads  to  a  new 
product.  At  eoae  stage  of  this  reaction  aoleeular  fragments  Bust  diffuse  over 
tbe  surface  to  fora  new  chemical  combinations.  Part  of  our  studies  has  been 
aiaed  at  exaaining  bow  diffusion  occurs  on  macroscopic  crystal  planes.  To 
carry  out  those  investigations  we  have  built  a  aoleeular  beaa  systea* 
illustrated  ia  Fig.  7.3*  ia  which  a  localised  deposit  of  gas  in  tbe  fora  of  a 
line  is  foraed  on  a  single  crystal  plane.  By  aoving  tbe  saaple  across  tbe 
fixed  eleotron  beaa  of  an  Auger  spectrometer*  tbe  concentration  &  of  adsorbed 
gas  can  be  deterained  as  a  inaction  of  botb  tiae  and  position.  In  principle 
it  should  therefore  be  possible  to  aeasure  tbe  transport  over  tbe  surface  and 
thus  to  establish  tbe  surface  diffusivity. 

A  start  in  this  effort  has  been  aade  with  preliminary  studies  of  nitrogen 
on  V(100).  Measurements  on  an  initial  deposit  of  half  a  aonolayer  are  shown 
in  Fig.  7.4.  In  diffusion  over  tbe  surface  we  expeot  a  general  redistribution 
of  tbe  gas*  leading  to  a  broadening  of  tbe  concentration  profile.  There  is 
little  indication  of  any  such  effect  ia  tbe  data.  The  peak  of  tbe 
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Temperature  (K) 


Fig.  7.2  Temperature  dependence  of  methane  decomposition  on  different 
faces  of  cungaten.  E  •  activation  energy  for  decomposition. 


Fig.  7.3.  Top:  Overview  of  mo  It  cu  Ur  bttm  system  for  turftct  diffusion 
studios.  Ssmpls  monipuUtor  is  at  top,  boom  systom  at  right. 
Bottom:  Datail  of  sampla  mounting,  as  soon  through  ball-jar  window 
Ho locular  beam  nossla  at  right,  alactron  bombardment  haator  and 
crystal  wafer  in  cantor. 


Hat  littto 


jps^wt*  Ox 


MOO 


©t»»r 


t  loo 


OX  id* 


tUM 


ry  to 


r.>^ 


7.  SURFACE  STUDIES 


Distance  (mm) 

Fig.  7.4  Concentration  profiles  for  nitrogen  on  W(100)  after  equilibra 
tion  at  1300°K. 


NORMALIZED  Si (91  «V)  AUGER  INTENSITY 
(ARB  UNITS) 


Fig.  7.5  Intensity  of  nitrogen  Auger  line  versus  substrate  intensity 
for  increasing  extent  of  nitridation. 
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Fig.  7 


Fig.  7 


6  Auger  scans  for  Si(100)  at  32  K  exposed  to  1  x  10“^  Torr  for  30  min 

(a)  No  electron  bombardment  during  oxygen  exposure 

(b)  Continuous  electron  bombardment 


.7  Auger  scans  for  Si(lOO)  at  room  temperature 

(a)  No  electron  bombardment  during  oxygen  exposure 

(b)  Continuous  electron  bombardment 
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oxidation  of  ailioon. 

Our  work  on  tka  spatial  stability  of  tksss  orsrlaysr  structures  will 
foeus  upon  surface  transport  as  the  inport ant  asekanisn.  Suck  transport  is 
quantitatively  okarsetsrised  by  tks  diffusivity  D  in  Pick's  second  lav  [7.11 , 
vkiek  for  one-dinensional  diffusion  appears  as 


9c  3  _  3c 
9t  3x  9x 


(1) 


One  of  tke  problsas  in  surface  diffusion  is  tkat  tks  diffusivity  is  generally 
not  a  constant  at  constant  temperature,  but  instead  varies  witk  concentration. 
Under  tkese  oircuastances  Eq.  (1)  kas  not  been  solved  in  closed  font  and  tke 
question  is:  How  can  D  be  deduced  fron  experiment? 

Using  tke  Auger  netkods  deecrlbed  in  Section  7.1*  tke  concentration 
profile,  tkat  is  tke  concentration  as  a  function  of  position,  can  be  obtained 
for  diffusion  intervals  of  different  lengths.  Fron  one  suck  profile  tke 
spatial  derivatives  esn  be  evaluated;  tke  tine  derivative  at  eonstant 
concentration  should  also  be  available  fron  several  suck  determinations. 
Using  computer  modeling  ve  have  begun  to  explore  bow  effective  this  technique 
is,  and  bov  muck  data  are  required  for  its  implementation.  In  Fig.  7.8  is 
shown  tke  form  of  tke  diffusivity  assumed  in  our  trials.  This  is  compared 
with  tke  values  of  D  derived  from  measurements  of  seven  concentration 
profilea.  From  this  preliminary  work  it  appears  tkat  our  technique  has 
promise.  It  still  remains  to  explore  bow  significant  experimental 
uncertainties  are  in  affecting  tke  derived  diffusivity,  snd  how  muck  data  are 
required  for  reliable  results. 


7.3  AtflBift  Exploration  fayi.tll  LftZUi* 

In  an  endeavor  to  better  understand  tke  formation  of  metal  overlayers  on 
solids  we  have  been  examining  tke  individual  atomic  steps  involved  in  suck 
proeesses  using  tke  field  ion  microscope.  Diffusion  of  individual  metal  atoms 
is  obviously  important;  we  have  shown  in  past  studies  tkat  atomic  diffusion  in 
seme  systems,  at  least,  may  be  quite  complicated  [R2] .  On  tke  (110)  plane  of 


*  Ibis  work  was  supported  by  tke  National  Science  Foundation  under  grants 
OMR  80-24061  and  82-01884. 
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Fig.  7.8  Concentration  depandance  of  dlffusivlty  derived  from  computer 
generated  concentration  profiles,  compared  with  Input  values. 
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««rt»U  ftet-snticid  cubic  metals,  atomic  diffusion  occurs  by  a  replaceaeat 
mechanism  in  which  the  sdatom  displaces  the  lattice  atom,  pushing  it  into  an 
adjacent  ohaunel  while  incorporating  itself  into  the  lattice.  According  to 
onr  studies  this  is  how  aotion  takes  plaee  on  the  (110)  plane  of  iridium.  Now 
it  ia  known  froa  work  on  large  seale  single  crystals  that  Ir(110)  is 
reeonstruoted  and  does  not  have  the  atoaie  arrangement  of  a  (110)  plane  in  the 
bulk.  In  contrast,  the  (110)  plane  of  rhodium  is  not  reconstructed,  and  on  it 
diffusion  is  normal.  Is  there  a  connection  between  oross-ehannel  diffusion 
and  surface  structure? 

Regrettably  the  atoaie  arrangement  of  the  planes  on  which  diffusion 
aeaenreaents  hare  been  aade  is  not  known ;  we  hawe  begun  an  examination  to 
better  define  the  atomic  arrangements.  From  measurements  of  the  spacing 
between  the  [llO]  rows,  we  know  that  Ir(110)  produced  by  low  temperature  field 
evaporation  has  undergone  a  (2  x  1)  reconstruction:  the  distance  between  the 
imaged  close-packed  rows  is  twice  the  normal.  This  still  does  not  define  the 
positions  of  the  individual  atoms  at  the  surface.  Ve  have  been  able  to  get 
some  idea  about  this  from  observations  of  a  sequence  of  field  ion  images 
produoed  by  evaporating  one  layer  of  the  surface,  imaging,  evaporating  another 
layer,  and  then  reimaging.  The  present  view  of  the  atomic  arrangement  of  the 
(110)  of  iridium,  as  produced  by  field  evaporation,  is  shown  in  Fig.  7.9.  The 
plane  is  seen  to  be  restructured  in  a  complicated  way  involving  a  considerable 
displacement  of  many  surf see-  and  near-surface  atoms.  The  extent  to  which 
other  fee  materials  undergo  similar  rearrangements  is  presently  being 
investigated. 

A  beginning  has  also  been  made  at  examining  the  actual  formation  of  an 
overlayer.  In  these  studies  rhenium  layers  on  W(211)  have  been  explored.  By 
ueing  a  single  rheninm  adatom  to  map  out  binding  sites  on  the  surface  it  has 
been  possible  to  establish  that  rhenium  layers  are  pseudomorphie  with  the 
substrate;  that  is,  their  structure  conforms  to  that  of  the  W(2ll)  surface. 
The  energetics  of  rheninm  atoms  bound  to  the  layer  have  also  been  examined. 
Here  most  unusual  results  have  been  obtained.  As  indicated  in  Fig.  7.10,  the 
binding  energy  of  an  adatom  ia  a  channel  adjacent  to  a  large  rhenium  cluster 
does  not  agree  even  qualitatively  with  our  standard  notions.  When  four  sites 
are  available  to  the  adatom,  the  binding  energy  at  the  two  center  positions  is 
higher  than  at  the  ends.  That,  of  course,  is  as  expected.  At  the  center 
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Fig.  7.9  Top-  and  side-view  of  reconstructed  (110)  plane  of  iridium 
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Fig.  7.10  Relative  binding  energy  of  rhenium  atom  ae  rhenium  cluster 
on  W(211) . 
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positions  the  naaber  of  neighbors  around  the  iditoa  is  largos  this  at  tk* 
ends;  however.  when  the  cluster  presents  firs  sdgs  sitss,  the  aitaatioa  ia 
qaito  difforont  -  tha  biadiag  of  tha  eaatas  aita  is  aetually  lass  than  at 
adjaoaat  positions. 

If  this  offset  proves  general  it  woald  hava  eonsidarabla  iapor tones  for 
rationalisinf  tha  propartias  of  ova r layer a,  and  work  on  atosi  bindings  ia 
tharaf ora  eontinning . 


[Rl]  CSL  Annual  Report  of  Progress  1980-1981,  7.1. 

[R2]  J.  D.  Vrigley  and  6.  Ehrlich,  Physical  Review  Letters,  vol.  44,  pp. 
661-463 ,  1980. 
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X.  Hayashi  E,  Mittra 

Graduate  jUldgall 

E.  Farr  S.  Ray 

R.  Kastner  T.  Triah 

8.1  Mill Inctsr-uave  htfiatri  Circuit a» 

la  the  pact  twelve-month  period,  we  lari  investigated  five  differaat 
aapaeta  of  millimeter-wave  iatagratad  eirenits  (MMIC) .  These  arc:  (i) 
multimode  oparatioa  of  MMIC  with  the  objaetiya  of  atiliziag  the  waveguide  and 
components  at  frequencies  that  arc  meh  higher  than  that  for  which  the 
waveguide  is  originally  designed;  (ii)  theoretical  and  experimental 
investigations  of  fin  linos  with  one  or  more  layers  of  dialaetrio  and  ferrite 
substrates;  (iii)  development  of  new  types  of  planar  waveguides  for  MMIC 
application;  (iv)  analytical  studies  of  discontinuities  in  open  waveguides 
using  the  newly  developed  spectral-domain  technique;  and  (v)  design  of 
waveguide  filters  by  introducing  periodic  discontinuities  in  the  planar  guide. 

Multlmoded  waveguides  are  espeeted  to  play  a  very  important  role  in  the 
design  of  sdbmil lime tor-wave  systems  operating,  say,  in  the  220  ffli  window. 
This  frequency  range  appears  to  be  well-suited  for  imaging  radar  application. 
Since  the  dimensions  of  the  waveguides  and  components  become  extremely  small 
at  frequencies  above  200  GIs,  it  is  worthwhile  to  determine  whether  a 
waveguide  structure  designed  for  the  E-Band  (80  (Si  range)  can  be  utilised  in 
the  220  Ob  range.  We  are  currently  performing  studies  on  different  aspects 
of  this  question. 

A  systematic  investigation  of  this  question  requires,  as  a  first  step, 

*  This  work  was  supported  by  Anqr  Research  Office  Grant  DAAG29-82-K-0084 
and  by  the  Joint  Services  Electronics  Program  (U.S.  Army,  U.S.  Navy,  U.S.  Air 
Force)  under  contract  N00014-79-C-0424. 


E.  Webb 
J.  Wilson 
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tlw  dmlofitat  of  an  officiant  ooaputar  prog  ran  for  the  datarninatiion  of  the 
propagation  constant,  and  tlia  oonpntation  of  tic  field  distribution  for  tha 
different  nodes  in  tbs  guide.  Such  a  progran  baa  bean  developed  and  tested 
for  tbe  fin-line  configuration,  vbieb  premises  to  be  a  good  candidate  for  both 
tbe  nil  line  ter-  and  sabnillineter-vave  ranges.  Tbe  progran  ia  able  to 
aeeonnodate  both  tbe  isotropio  and  anisotropic  (e.g.,  ferrite)  substrates. 
Anisotropic  substrates  are,  of  course,  employed  in  tbe  design  of  nonreeiproeal 
components,  sucb  as  isolators.  Tbe  tbeoretical  results  for  tbe 
ebaracteristies  of  tbe  fin  line  are  also  being  verified  experimentally.  Plans 
for  testing  tbe  nult inode  waveguides  are  currently  underway.  Pr el in inary 
reaults  for  tbe  propagation  neasurenents  look  encouraging  and  a  nore  thorough 
investigation  is  currently  being  planned. 

Along  with  tbe  fin  line,  we  are  also  investigating  other  types  of 
shielded,  planar  waveguide  structures  which  appear  to  have  sane  unique 
advantages.  One  of  these  is  tbe  suspended  B-guide  shown  in  Figure  8.1.  Its 
advantages  are  low  netal  lose,  total  shielding  and  em^atibility  with  other 
eonponents  fabricated  by  using  tbe  suspended  E-guide.  A  theoretical  analysis 
of  this  guide  has  been  conpleted  and  nanerical  etadiee  are  currently  being 
perfomed.  The  potential  use  of  tbe  suspended  B-guide  as  a  building  bloek  for 
MMIC  is  also  being  investigated. 

An  important  theoretical  prOblen  arising  in  tbe  developuent  of  MMIC 
components  is  tbe  understanding  of  tbe  behavior  of  discontinuities  ia  open 
waveguides  and  their  circuit  characterisation.  Conventional  approaches  to 
analysing  closed  waveguide  discontinuities  are  not  very  useful  ia  tbe 
investigation  of  open-region  problems ,  beesuse  tbe  open  waveguides  not  only 
support  tbe.  usual  bound  nodes,  but  an  infinite,  continuous  spectrua  of 
unbounded  nodes.  We  have  been  developing  a  new  approach  based  on  tbe  spectral 
doaaia  technique  which  appears  to  be  well-suited  for  geoaetries  which  are 
unbounded  ia  nature.  When  fully  developed,  tbe  netbod  is  expected  to  fill  an 
important  gap  ia  tbe  array  of  theoretical  or  aunerieal  tools  available  for  tbe 
analysis  of  open  region  discontinuity  problens. 

Other  eonponents  which  we  have  studied  sre  filters  sad  periodic  and 
leaky-wave  antennas.  We  have  investigated  two  different  kinds  of  filter 
designs  -  array  of  sandwiched  dissimilar  dielectric  sheets  sad  periodic  notch 
discontinuities  in  a  planar  dielectric  waveguide.  Ia  the  recent  past,  we  have 
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printed  metallic  strip  diseoat innit is s  frost  the  point  of  view  of 
designing  leaky-wave  sntsnnss.  Both  thsorstiesl  snd  szperinentsl 
investigstioas  have  boon  esrrisd  ont  snd  extensive  results  have  boon  obtained. 
Future  work  will  be  directed  toward  designing  waveguide  filters  using  periodic 
strip  discontinuities. 

8.2  Bleatramenetia  ggitttllll  1A4  Mil  llXIil  Identification* 

We  have  been  developing  s  new  spectral  doaaia  technique  for  solving 
electromagnetic  scattering  problems  involving  arbitrary  seatterera  which  nay 
either  be  dielectric  (lossy  or  losslsss)  or  perfectly  conducting.  This  aethod 
has  been  devised  to  serve  as  an  alternative  to  Moment  Methods  and  Asyaptotio 
Techniques  for  solving  EM  scsttsring  probleas  in  the  intermediate  frequency 
range.  The  aethod  takes  advantage  of  the  siaplieity  with  which  the  planar 
source-field  relationships  are  expressed  in  the  spectral  doaain.  The  boundary 
conditions  or  constitutive  relationships,  on  the  other  hsnd,  are  expressed 
aost  simply  in  the  spatial  domain.  Alternating  between  the  two  doaains  is 
carried  out  with  the  FFT  algorithm.  A  technique  that  has  been  applied  in  the 
past  to  planar  structures  has  been  generalised  to  aeeoaaodate  three- 
dimensional  bodies.  This  is  accoaplished  by  ssapling  the  current  distribution 
on  the  scatterer  over  a  nnaber  of  parallel  planes,  snd  using  the  staple 
spectrsl-doaain  interaction  relationships  between  the  planes.  This  new 
approach  is  applicable  to  arbitrarily  shaped  conducting,  dielectric  or  lossy 
dielectric  seatterera.  It  has  a  lsrgsr  dats-hsndling  capacity  than  that  for 
Moment  Methods  snd  is  more  general  than  asymptotic  techniques.  Thus,  it 
provides  an  efficient,  auoh-needed  approach  to  filling  the  gap  between  the 
low-  sad  high-frequency  conventional  techniques. 

During  ths  psst  year,  sdditionsl  work  has  been  done  on  the  problea  of 
radar  target  identification  by  the  aethod  of  determining  the  spatial 
distribution  of  scattering  centers.  It  has  bean  found  thst  both  specular 
points  and  sdgs  diffraction  points  css  be  located  snd  catsgorised  by  analyzing 
RC3  data  with  ths  Prony  algoritha.  A  aodified  Prosy  algorithm  has  been 
developed  to  handle  the  esse  of  split  spectrua  aessureaents,  i.e.,  when  RCS 

*  This  work  was  supported  by  Office  of  Navsl  Resesrch  Grant  N00014-81-K- 
0249  sad  by  ths  Joint  Services  Electronics  Prograa  (U.S.  A ray,  U.S.  Navy,  U.S. 
Air  Force)  under  contract  N00014-79-C-0424 . 
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data  arc  gathered  ia  2  or  3  frequency  bands  which  arc  not  contiguous.  Va  haws 
shown  that  aany  of  the  advantages  of  wide-band  aeasuraaents  oan  ba  obtained  by 
eoa&iaiag  data  froa  saaller  (and,  hanea,  aora  praetieal)  snbbands. 

Progress  has  bean  aada  towarda  answering  the  question  of  whether  or  not 
the  scattering  center  description* is  still  valid  when  corners  are  rounded  or 
costed  with  dieleetries.  Ia  order  to  deteraiae  this,  the  forward~scattering 
problea  aust  first  be  solved.  Coaputer  prograas  based  on  the  speotral 
iterative  technique,  described  in  the  last  paragraph,  are  eurreatly  being 
developed  for  solving  the  sbove-aentioaed  forward  probleas. 


61 


9.  PLASMA  PHYSICS* 


Facility  ia£  Senior  Staff 


M.  Raether 


6r adnata  Stndant 


S.  V.  Ault 


9.1  Statiatical  Properties  al  flllll  IftRlfflltatf 

The  objectiT#  of  thia  research  la  to  study  the  statiatioal  proportion  of 
plasaa  turbulence.  The  systoa  under  investigation  ia  the  positive  colon*  of  a 
dc-dischargo  in  holiun  which  is  unstable  with  respect  to  the  spontaneous 
excitation  of  ion  acoustic  waves.  Bleotrie  field  fluctuations  of  the 
turbulent  wave  field  are  picked  up  hy  *  Langmuir  probe  and  their  probability 
distribution  is  Measured  by  a  fast  stapling  aethod.  Significant  deviations 
froa  Gaussian  behavior  are  noticed  in  the  unstable  regiae.  This  together  with 
evidence  derived  froa  the  speetrun  points  to  possible  phase  and/or  aaplitude 
correlations  between  waves  of  the  hind  that  is  typical  for  wave-wave  coupling. 
This  allows  us  to  draw  conclusions  as  to  the  nature  of  the  phase  and  aaplitude 
enseable  of  the  Fourier  eoaponents. 


•This  work  was  supported  by  the  university  of  Illinois 
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Faenltr  Senior 

S.  M.  Yen 
flalMtl  Stad«»ti 

R.  OabMUii  D.  R.  Hall  S.  H.  Laa 
10.1  llglljLsA  Gaa 

Rarefied  gaa  dyaaaies  daala  with  noa-aqnilibrlna  gas  flow  problaas  la 
whiek  aicroscopie  traataaat  aoeordiag  to  kiaatio  theory  is  accessary  to 
datacalaa  the  offset  of  iateraoleealar  eolliaioaa  aad  gas  sorfaee  iataraetioas 
oa  both  tha  ■ieroseopie  aad  aaeroseopie  gas  flow  propartias.  Saab  rarafiad 
gaa  flow  problaas  oaaar  aot  oaly  la  aerodyaamies.  bat  also  la  alaotroaios. 
aeroaoay,  enriroaaentai  flald  dyaaaies.  aad  other  ralatad  f la Ids. 

Tka  ala  of  this  rasa ar ok  progs aa  Is  to  dsYSlop  oaaar ioal  aatkods  to  solra 
a  vlda  raago  of  problaas  aad  ar  eoadltloas  far  froa  aad  aaar  thsraal 
aqalllbrlaa.  A  Moats  Carlo  aatkod  has  baaa  dare loped  at  tka  Coordlaatad 
Selaaoa  Laboratory  [Rl]  to  so Its  diraetly  tha  Boltsaaaa  aqaatioa  aad  has  baaa 
asad  by  Ota  Boltsaaaa  groap  to  solra  tha  Boltsaaaa  aqaatioa  for  ssroral 
rarafiad  gas  flow  problaas  aadsr  a  vlda  raaga  of  aoasqalllbrlaa  aad  boaadary 
eoadltloas  [12-R7] .  Tha  solatloas  va  hara  obtalaed  ylaldad  datallad 
alorosooplo  aad  aaerosoople  aea-aqalllbrlaa  propartias.  aost  of  which  lure 
aarar  baaa  traatad  aad  stadlad  bafora.  Va  hara  also  stadlad  oaaar leal 
solatioas  of  othar  klastlo  aqaatloas  aad  othar  aaaarleal  aethods  to  solra 
rarafiad  gas  flow  problaas.  laeludiag  tha  dlraet  slaalatloa  tsehaiqua. 

Noallaaar  araporatioa-eoadeasatloa  problaas  era  aaeooatsrsd  la  sack 
dlrarslflad  arasa  as  appar  ataosphara  astaorology*  tha  eoollag  of  oaelaar 
reactors,  deslga  of  spaee  exparlaaats.  patroehaalesl  eaglaeeriag,  raeaaa 
technology,  aad  the  lateraetloa  of  high  power  laser  radlatloa  with  aatal 
sarfaees.  Tha  traataaat  of  those  problaas  raqalras  first  tha  eoosldaratloo  of 
rapor  klaetles  problaas  eharaeterissd  by  tha  aoaeqallibrloa  rapor  aotioa  la  a 
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inudaoa  layer  it  tki  latMflui  beudirj.  Under  i  Joint  research  of  fort*  wo 
kw*  successfully  Mini  tk«  evaporation  problem  [K8] .  Our  Boltxnann 
*o luti OB*  utikliA  the  validity  of  the  book  ooUtioa  and  Ytrehas'  kinetic 
theory  mnuk  lit] .  Ytrehas'  approach  eon  b*  naad  to  oalonlata  tba  jimp 
eoaditioaa  aoroaa  tba  bndaon  layar  and  tba  not  naaa,  aeantu  and  boat  fluxes 
tlO.1]  and*  together  with  Ike  conventional  eontinnan  natbod*  to  oalonlata  tbe 
flow  pnranetera  at  tbe  enter  edge  of  tbe  Sandses  layer  for  pr  obi  ana  in  wbieb 
evaporation  oeonra  at  tke  interpkaae  bonadary. 

Stndiea  bane  alao  been  aade  for  tke  condensation  problon  [10.2].  Tbe 
noa-equil ibr ian  vapor  near  a  eoadenaing  aarfaee  differa  froa  that  near  an 
evaporating  aarfaee  and  ita  behavior  ia  more  ooaplex.  Oar  Boltxnaaa  aolntions 
yield  distiset  aoa-equil ibr ian  behavior  in  tba  Eandaan  layer  of  eoadenaing 
vapor .  Ike  eoadenaation  rate  as  a  function  of  pressure  ratio  waa  found  to 
depend  atrongly  on  tke  anbataaoe  parameter  fi  which  la  proportional  to  tbe 
ratio  of  tbe  latent  beat  and  ita  liquid  aarfaee  temperature.  Sene  of  oar 
reunite  were  found  to  be  in  agreeneat  with  those  obtained  by  using  tke  staple 
approaches  of  Qgnehi  [Sill]  and  Ytrekua  [111].  ® 

Vo  plan  to  study  farther  tke  evaporation  aad  the  eoadenaation  problaus  by 
considering  acre  ooaplex  flow  rad  bonadary  conditions. 

10.2  Naval  a «* 

tree  aarfaee  wave  pr obi eas  aacoanterad  in  naval  hydxodynaaics  are 
ekaraeteriaed  by  complexities  in  flew  gecnetry*  flow  features  rad  boundary 
conditions.  Ike  flow  has  aa  nnkneon  free  surface  rad  it  is  propagative  and 
transient.  Bn  boundary  ooadltion  at  the  free  aarfaee  is  of  a  nixed* 
parabolic  type  and  it  eontains  highly  nonlinear  terns.  Ia  tbe  steady  state* 
tkere  also  exists  a  radiation  boundary  condition*  since  the  waves*  once 
deaerated*  propagate  downs tr era.  These  complexities  have  led  to  several 
computational  difficulties:  accurately  accomodating  the  free  surf  see 
geometry*  satisfying  tke  Irani ary  condition  uniformly  ever  tbe  free  surface* 
creating  tke  radiation  boundary  condition  rad  impl masting  tbe  boundary 
condition  at  tke  contact  line  of  ike  free  surface  with  the  solid  surface  of  a 


*  Ibis  work  was  supported  by  tke  Offiee  of  Naval  ftsscaroh  under  Contract 
MDM24-«HC-d74d. 
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partially  subaerged  body. 

la  tha  initial  phase  of  oar  study,  we  have  developed  two  tine-dependent 
nuaerieal  seheaes  with  Eulerian  grid  systeas  for  solving  steady  and  unsteady 
potential  flows  for  nonlinear  free  surfaee  probleas.  In  one  aethod 
[R12,R13,R15] ,  we  use  the  finite  eleaent  aethod  to  deal  with  the  geometrical 
coaplexity  and  the  free  surfaee  boundary  condition  of  the  nonlinear  free 
surfaee  probleas.  In  this  sehene,  the  finite  eleaent  aethod  is  used  to  aake 
the  field  ealeulation  and  the  finite  difference  aethod  is  used  for  the  tiae 
evolution.  Ve  have  used  this  aethod  to  solve  two  probleas:  a  pressure 
distribution  aoving  with  a  constant  speed  and  a  aoving  submerged  elliptical 
cylinder  or  a  syaaetrical  hydrofoil. 

Ve  have  also  developed  an  explicitly  tine-dependent  finite  difference 
scheae  [R14.R16].  Ve  have  used  this  aethod  to  solve  two  probleas:  a  pressure 
distribution  aoving  with  a  conatant  speed  and  an  elliptic  cylindrical, 
surface-piercing  strut  accelerating  froa  rest.  Our  solutions  serve  to 
deaonstrate  the  feasibility  of  using  the  aethod  to  solve  two-diaensional  as 
well  as  three-diaensional  probleas. 

The  results  of  the  initial  phase  of  our  study  were  presented  at  the 
Second  International  Conference  on  Nuaerieal  Ship  Hydrodynaaies  [116], 

In  both  nuaerieal  seheaes  aentioned  above,  the  conputational  doaain  is 
expanded  downstreaa  periodically  during  the  eoaputation  as  the  disturbance  on 
the  free  surfaee  is  propagated  close  to  the  downstreaa  boundary.  Therefore, 
the  undisturbed  condition  is  applied  on  the  cut-off  downstreaa  boundary. 

In  order  to  obtain  the  longer  tiae  and  steady  state  solutions,  we  use  a 
fixed  downstreaa  boundary  set  close  to  the  disturbance  so  that  we  can  increase 
the  eoaputation  efficiency  by  dealing  with  a  saall  conputational  doaain.  The 
seoond  phase  of  our  study  has  been  focused  on  the  application  of  our  nuaerieal 
seheaes  to  the  nonlinear  free  surfaee  probleas  with  a  fixed,  saall 
computational  deaain.  Ve  have  studied  two  approaches.  In  each  approach,  we 
have  used  Or Isa ski 'a  aethod  [R17.R18]  to  iapleaent  the  open  boundary  condition 
at  the  outflow  boundary.  The  two  approaches  differ  in  the  treataent  of  high 
frequeney  errors  that  appear  in  the  solution  and  that  inhibit  accurate 
treataent  of  the  open  boundary  condition.  In  the  first  approach,  we  used 
filtering  to  aiaiaixe  high  frequency  errors.  However,  this  was  found  to  be 
effective  only  at  an  early  tiae;  therefore,  it  is  suitable  only  for  obtaining 


10.  RAREFIED  GAS  DYNAMICS  AND  COMPUTATIONAL  GAS  DYNAMICS 


65 


trat..ient  solutions.  In  the  second  approach,  ve  used  daaping  in  our  tine 
marching  schemes.  Damping  was  found  to  be  more  effective  in  controlling  high 
frequency  errors.  Transient  and  steady  state  solutions  have  been  obtained  for 
the  pressure  distribution  and  the  accelerating  strut  problems. 

The  results  of  the  second  phase  of  our  work  were  presented  at  the  'Third 
International  Conference  on  Numerical  Ship  Hydrodynamics  [R19] . 

Ve  are  now  in  the  third  phase  of  our  research.  The  objective  of  the 
study  in  this  phase  is  to  apply  our  method  to  the  ship  wave  problem.  Ve  have 
chosen  the  ship  of  Vigley  shape  as  our  first  problem.  The  hull  shape  of  this 
ship  is  shown  in  Fig.  10.1.  The  numerical  method  we  use  has  the  following 
features: 

(1)  the  mesh  system  accommodates  both  the  free  surface  and  the  hull  surface, 

(2)  the  nonlinear  free  surface  condition  is  satisfied  accurately,  and 

(3)  the  computation  efficiency  for  the  field  calculation  is  increased  by 
using  the  successive  line  over-relaxation  method. 

The  computer  program  to  implement  our  method  for  the  Vigley  ship  has  been 
written  and  preliminary  calculations  for  the  early  times  have  been  made. 

10.3  AfrPdYUMiigf 

The  objective  of  one  of  our  studies  for  aerodynamic  problems  is  to  solve 
directly  the  Euler  equation  in  the  entire  compressible  flow  regime:  subsonic, 
transonic  and  supersonic.  There  are  several  reasons  we  wish  to  do  so.  One  is 
to  consider  the  rotationality  in  the  flow.  The  other  is  to  facilitate  error 
study  of  the  calculation  since  the  same  basic  equation  and  the  same  numerical 
method  are  used  in  the  entire  flow  regime.  Our  initial  effort  is  to  study  the 
numerical  solution  of  compressible  flow  over  an  axi-symaetric  triconic  body. 
Ve  seleoted  the  implicit  factored  seheme  which  was  developed  by  Varming  and 
Beam  [R20]  to  solve  the  Navier-Stokes  equation.  Ve  used  the  AIR3D  Navier- 
Stokes  computer  program  developed  at  NASA/ AMES  and  adapted  by  Nietubiex  [R21] 
for  axi-symmetric  flows.  Ve  developed  and  used  a  mesh  generation  system 
designed  for  the  trioonic  body.  In  addition,  we  revised  the  implementation  of 
several  boundary  conditions. 

Ve  have  obtained  the  Euler  solution  for  -  0.5  to  2.75.  In  order  to 
speed  up  the  convergence,  we  use  the  solution  for  a  lower  Mach  number  to 
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initialize  the  solution  of  a  higher  Mach  number.  The  nosh  system  bad  to  bo 
rsfinsd  whenever  necessary  for  the  oalenlation  of  flow  for  a  given  Mach 
number. 

Fig.  10.2  sbovs  the  surface  pressure  distribution.  vs  z  for  1^  •  0.5 , 
0.95.  1.2  and  2.75.  These  results  are  compared  with  calculations  by  Gustafson 
[R22].  who  used  a  panel  method  for  potential  flow  for  -  0.5.  relaxation 
method  for  transonic  flows  for  -  0.95  and  1.2  and  a  finite  difference 
method  to  solve  the  Euler  equation  for  the  shock  layer  for  Ma  a  2.75.  The 
local  Mach  number  contour  lines  is  the  physical  domain  for  Mm  a  0.95.  1.2  and 
2.75  are  shown  in  Fig.  10.3. 
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FlSfiUgatl  Fault  Modal  a  Xlfll  ZlilltU* 

Typical  physical  failmrcs  observed  is  the  field  were  s inula ted  at  the 
circuit  level  om  basic  AMDS  cad  OIOS  eelle  to  study  their  effect  at  the  logie 
level.  la  have  fouad  that  aaay  physical  failures  could  not  be  Modeled  by 
exiatiag  fault  aodals  such  as  the  stuek-at  Models.  The  results  of  these 
detailed  atudiee  were  used  to  develop  a  Multivalued  algebra  which  eould  be 
used  to  prediet  MOS  circuit  behavior  under  physical  failures  for  large 
■odulea.  Using  the  rules  of  the  algebra*  a  siaulator  has  been  i^pleaented  in 

*  This  work  was  supported  by  VHSIC  (U.S.  Navel  Electronics  SysteM  CoMMand) 
under  contract  nunber  N00039-80-C-0556 . 
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FISCAL  under  TAX/ UNIX;  tkii  is  about  two  orders  of  magnitude  faster  then 
circuit  simulation.  Tkese  studies  ksre  also  bees  used  to  derive  functional 
level  f salt  models  for  modules  smob  as  Programmable  Logie  Arrays  sad  Decoders 
which  have  soma  structural  regmlarity. 

From  tbe  eireait  simulation  study,  ve  bare  also  derived  approx las te 
formulas  for  estiastimg  tbs  obaages  la  tbe  iaversioa  voltage  sad  gate  delay 
doe  to  a  shift  ia  tbe  threshold  voltage.  Oaoe  tbe  behavior  of  a  fault-free 
eireait  is  detemiaed  by  eireait  aiaslatioa.  tbe  foraslas  eaa  be  ased  to 
estimate  tbe  obaages  ia  tbe  eireait  behavior  due  to  a  threshold  faalt. 

11.2  Test  gfatmlffll  lax  Commlea  Dialtal  Systems* 

As  iategrated  eireoita  iaerease  ia  complexity,  tbe  problem  of  testiag 
them  ia  beeomiag  more  sad  more  diffiealt.  Ve  have  fooad  that  oar  approaeb  of 
deseribiag  physical  feilares  at  the  faaotional  level  aad  deriviag  tests  at 
this  level  is  aa  exeel leat  solatioa  to  this  problem.  This  technique  baa  beea 
saeeessfally  applied  to  deriviag  teats  for  memories,  programmable  logie 
arrays,  aad  regalar  logie  straetares.  A  new  aad  effieient  memory  test 
algorithm  wbieh  will  detect  iateraetioaa  between  multiple  eells  ia  a  memory 
has  beea  derived;  previous  work  bad  restricted  tbe  number  of  interacting 
cells  to  two  or  three.  Ve  have  also  found  a  method  to  design  easily  testable 
logie  straetares  directly  from  a  high-level  language  deseriptioa  of  the 
desired  function.  Tbe  inherent  structure  is  a  tree  wbieh  provides  a  high¬ 
speed  realisation  of  tbe  function.  Aa  effieient  algorithm  for  testing  general 
tree  straetares  has  also  been  found. 

The  current  tread  in  designing  complex  systems  is  toward  tbe  use  of 
struotured  techniques,  aad  this  has  motivated  tbe  design  of  a  faalt  simulation 
technique  which  takes  advantage  of  struotured  design  in  order  to  redaee  tbe 
eoatplexity  of  simulation.  A  concurrent  faalt  simulator  which  hierarchically 
evaluates  the  input  design  has  beea  designed  and  is  being  coded  to  evaluate 
the  performance. 


*Thia  work  was  supported  by  the  Joint  Services  Eleetronies  Program  (U.S. 
Army,  U.S.  Navy,  aad  U.S.  Air  Force)  under  contract  N00014-79-C-0424. 
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11 .3  fiauagui  Error  Detoetioa  Taahaloaea* 

A  faaetloaal  orror  model  has  teas  foraalstod  for  inforastioa  traasfor  sad 
•tor ago  la  s  mieroprogrammeble  oaatrol  malt,  sad  oodiag  sad  addrosslag 
toohaiqaos  bars  beea  developed  far  do too flag  those  errors  eoaeurreatly  with 
aormal  oporstioa.  Those  teehmifues  are  also  applicable  to  doteetiag 
addrosslag  fsilaros  la  asaorjr  systems. 

Vo  have  developed  aa  error-do toot loa  s oh saw  based  oa  time  reduadaaey; 
this  scheme,  oallod  leooapmtiag  with  Shifted  Oporaads  (1ESO) ,  is  applioablo  to 
both  arithaetie  aad  -logio  operatioas.  Xt  has  booa  oztoadod  to  iaolade 
multiply  sad  divide  arrays.  This  scheme  has  also  booa  oztoadod  to  oorroot  all 
faaetloaal  errors  ia  a  typioal  adder  ssiag  time  reduadaaey  whoa  the  errors  are 
ooafiaod  to  a  saall  area  of  the  ship. 

A  study  of  prograaauiblo  logio  arrays  (FLAs)  has  foaad  that  physieal 
fsilaros  will  rosalt  ia  aaidirsetioaal  errors  at  the  oatpats.  Oesiga 
procedares  have  booa  foaad  to  dosiga  PLAs  with  eoaearreat  orror  detoetioa 
asiag  naidireotioaal  orror-do toot lag  codes  saeh  as  Berger  oodes  or  a-oat-of-n 
codes. 

11 .4  fiiiiu  al  Hrtlltl  Milk  Jillfil.  Detoetioa* 

le  have  desigaed  a  parsllel-pipoliaed  aodale  which  perforas  the  battorfly 
oporstioa.  the  coatrol.  sad  the  data  storage  of  the  Fast  Foarior  Traasfora. 
Eight  battorflios  are  perforaed  oa  oae  eight-bit  slice  chip.  Coaoarreat  error 
detoetioa  was  achieved  by  1ES0  for  the  ALU.  parity  for  the  aoaory  aad  data 
paths*  sad  by  daplicatiag  the  coatrol.  The  chip  has  boea  laid  oat  ia  aMOS 
teehaology. 

A  residae  aaaber-eaeoded  digital  filter  architeetare  has  beea  desigaed 
aad  is  boiag  laid  oat  for  aMOS  fabrieatioa.  It  is  orgaaized  as  a  multiple 
stress*  pipeliaed  straetare  si low lag  high  device  atilizstioa  aad  modularity. 
The  residae  system  si lows  fsalt  toleraaee  sad  grscefal  degrads tioa. 

Aa  architeetare  sad  sa  orgsaisstioa  have  beea  defiaed  for  a  asefal 
multiple  iastraetioa  stream  shared  pipeliae  processor  implemeatable  oa  a 
siagle  chip. 


•  This  work  was  supported  by  VHSXC  (U.S.  Naval  Electromics  System  Commaad) 
aader  eoutrset  N00039-80-C-0556 . 
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A  cell  for  •  systolic  srrsy  processor  has  been  designed  end  is  being  laid 
out  using  a  CMOS  process  in  cooperation  with  ESL,  Inc.  It  incorporates  a 
two-stage  pipeline  and  uses  RESO  and  parity  for  concurrent  error  detection. 
Its  primary  application  is  for  Finite  Impulse  Response  filtering. 

An  address  generator  chip,  which  will  rapidly  generate  memory  addresses 
for  array  elements  in  loops  under  control  of  a  host  processor,  is  being 
designed  and  laid  out  using  CMOS  rules  in  cooperation  with  ESL,  Inc.  It 
allows  up  to  eight  indices  in  the  array.  The  appropriate  atride  constants  are 
stored  in  the  chip  prior  to  entering  a  loop. 

The  design  and  nMOS  implementation  of  a  microprogram  control  unit  with 
control  capability  similar  to  the  AMD  2910  microprogram  sequencer  have  been 
undertaken.  The  concurrent  error  detection  techniques  developed  in  the  group 
for  microprogram  control  units  are  being  implemented  in  this  design. 

A  design  in  nMOS  technology  of  a  crossbar  switch  which  could  be 
potentially  used  in  a  multistage  network  of  the  delta  type  has  been  completed. 
Two  types  of  chips,  a  control  chip  and  a  data  chip,  have  been  laid  out  for 
this  purpose.  Ve  are  now  working  on  incorporating  concurrent  error  detection 
in  the  module. 

A  complex  mult ip ly-ac cumulate  chip  with  concurrent  error  detection  has 
been  laid  out  in  nMOS  technology.  This  chip  uses  RESO  for  concurrent  error 
detection. 

11.3  Faplt-lS.UXJUal  M&iliA  AUPgAthW  Is !£  Processor  Arrays* 

Classical  techniques  for  fault  tolerance  are  based  on  replications  of 
modules,  with  the  errors  from  a  failed  module  being  masked  by  the  correct 
outputs  from  the  good  modules.  Ve  have  found  very  efficient  techniques  for 
fault  tolerance  in  matrix  operations  where  the  inherent  redundancy  in  the 
algorithms  executing  on  array  processors  is  used  to  detect  and  correct  errors 
from  a  failed  processor.  The  techniques  use  the  idea  of  checksums  in  a  novel 
fashion,  and  the  redundancy  required  is  only  of  the  order  of  (1/n) ,  where  n  is 
the  size  of  the  array.  The  method  oan  be  used  either  with  mesh-connected  or 
systolic  arrays.  By  tailoring  the  fault-tolerance  technique  to  the  algorithms 


*  This  work  was  supported  by  VHSIC  (U.S.  Naval  Electronics  System  Command) 
under  contract  ND0039-80-C-0353 . 
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u«d  tad  to  tk<  itrMtv*  of  the  array  processo »  n  it*  able  to 
achieve  error  dotootioa  oaf  correction  at  a  ooat  far  last  tkaa  that  of  aaiag  a 
faaaral  replication  techaiqae . 

11.6  A  Soar » a  Matrix  Processor* 

Tba  efficiency  of  spazaa  matrix  eompatstions  oaa  bo  dramatically  improved 
by  aaiag  matrix  sparsity  pa tt arms  to  alimimata  eomditiomal  braaokas  aad  array 
iadax  ealealatioas.  Va  have  axtaadad  tkis  taekaiqaa  for  highly  eoaearraat 
parallal  proeassors.  Static  sahadaliag  of  sparsa  matrix  eoapatations  at 
compilation  tima  rasalts  ia  lowar  overhead  aad  kigkar  rasoarca  atilixatioa. 
Oaaaratiag  addrassas  aad  othar  eoatrol  iaformatioa  at  tha  sita  where  it  is 
asad  radaoas  tha  amoaat  of  iatarehip  eoomaaiaatioa  dariag  axaeatioa.  A 
compiler-architecture  system  asiag  this  approach  has  baaa  proposed  aad 
evaluated  for  its  parformaaca  oa  sparsa  LU  decompositioa  aad  sparsa  matrix 
maltiplicatioa.  High  rasoarca  atilixatioa  (over  80%)  was  obs erred  for 
hardware  coaf igaratioas  with  ap  to  128  proeassors  wkea  largo  matrices 
(400x400,  S  clamant*  pax  row)  ware  asad. 

11.7  MptUaftf  Architectures** 

We  hare  baaa  able  to  eoastract  mathematical  models  of  imstraot ioa-aai t 
pipaliaas  with  probabilistic  flow-perturbing  axcaptioa  eoaditioas.  These 
models  are  valid  over  wide  ranges  of  system  parameters.  A  traoe-processiag 
algorithm  which  yields  exact  parformaaca  maasaramaats  of  programs  which  are 
sob j act  to  data  and  eoatrol-flow  depeadeaeies  ia  pipaliaas  has  baaa  developed, 
aad  aa  exact  residaal-time  modal  has  baaa  derived  for  hardware  depeadeaeies  ia 
pipaliaas.  These  are  the  first  kaowa  models  for  iastraetioa  pipaliaas  which 
eaa  affectively  evalaate  parformaaca  aadar  affects  sack  as  cache  miss, 
rasoarca  coafliet,  error  retry,  breach  aad  data  depeadeaey,  sad  faaetion- 
dapaadaat  segmeat  tima  variatioa. 

A  stady  of  the  recovery  of  pipeliaed  systems  from  transient  f salts  has 

*  This  work  was  sapportad  by  tha  Joist  Services  Electronics  Progrsm  (U.S. 
Army,  U.S.  Navy,  U.S.  Air  Force)  aadar  contract  N00014-79-C-0424 . 

**  This  work  was  sapportad  by  the  Joint  Services  Electronics  Progrsm  (U.S. 
Any,  U.S.  Navy,  U.8.  Air  Force)  aadar  contract  N00014-79-C-0424  and  by  VSHIC 
(U.S.  Naval  Electronics  System  Command)  aadar  contract  N00039-80-C-0356 . 
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produced  an  approximation  for  the  distribution  of  faulty  periods  which 
compares  very  favorably  with  simulation. 

11.8  Structured  Memory  A £.8811  Architecture* 

Our  Structured  Memory  Access  (SMA)  architecture  represents  a  significant 
breakthrough  of  the  "won  Neumann  perforaunee  bottleneck"  by  reducing  memory 
wait  time  and  by  eliminating  a  majority  of  the  memory  references  required  to 
execute  a  program.  This  breakthrough  is  achieved  by  decoupling  the  address 
generation  process  from  the  data  computation  process.  A  simulation,  running  a 
Gaussian  elimination  program,  showed  the  performance  to  be  a  linear  function 
of  the  pipeline  lengths  of  the  Memory  Access  Processor,  the  memory,  and  the 
Computation  Processor.  The  SMA  achieves  about  a  50%  reduction  in  memory 
references  over  conventional  auehines  running  optimised  code. 

11 .9  Multiuser  Multiprocessors  with  ftlvtf  Cache  Memories** 

In  the  context  of  a  multiuser  multiprocessor  system  with  private  cache, 
we  have  considered  the  write- through  versus  the  write-back  policy  of  swin 
memory  update.  The  advantage  of  the  write-back  policy  is  that  the  bus  traffic 
is  reduced  in  comparison  to  the  write- through  policy.  It  is  usually  assumed 
that  the  coherence  problems  of  write-back  require  hardware  such  as  global 
directories  to  detect  potential  coherence  problems.  For  this  reason,  a 
write-through  cache  is  usually  used  which  provides  coherence  for  all 
transactions.  Ve  have,  however,  suggested  ways  to  avoid  these  coherence 
problems  in  user  code  and  have  examined  the  potential  savings,  in  terms  of  bus 
traffic,  by  using  a  write-back  rather  than  a  write-through  cache.  Using  a 
detailed  instruction  level  simulation,  ve  found  that  in  the  typical  case  the 
write-back  policy,  w«pared  to  write-through,  will  allow  more  than  twice  the 
number  of  processors  on  the  bus  at  a  given  traffic  level. 


*  This  work  was  supported  by  VSHIC  (U.S.  Naval  Zl^tronics  System  Cosuaand) 
under  contract  N00039-80-C-0536 . 

**  This  work  was  supported  by  the  Joint  Services  Electronics  Program  (U.S. 
Army,  U.S.  Navy,  U.S.  Air  Force)  under  contract  N00014-79-C-0424  and  by  VSHIC 
(U.S.  Naval  Electronics  System  Cosunnd)  under  contract  N00039-80-C-0556. 
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11 .10  DiliftJi  Verification  Techniques* 

A  top-down  hierarchical  design  snt oast ion  facility  is  being  developed 
which  includes  design  verification  capability.  A  graphics  input  package  has 
been  interfaced  to  the  Liveraore  SCALD  II  package*  and  a  aechanisa  for 
including  functional  descriptions  of  aodules  has  been  developed  for 
verification.  An  autoaated  reasoning  program  developed  by  Argonne  National 
Laboratory  is  being  used  to  deteraine  equivalence  between  the  functional 
description  of  a  nodule  and  the  structure-driven  coabination  of  the  functional 
descriptions  of  its  subaodules. 


*  This  work  was  supported  by  VSHIC  (U.S.  Naval  Electronics  Systea  Coaaand) 
under  contract  ND0039-S0-C-0SS6 . 
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12.1  Introduction 

Current  technology  requires  the  development  of  efficient  computational 
techniques  and  the  analysis  of  the  capabilities  of  various  models  of 
computation.  Ve  are  concerned  with  the  resources  -  such  as  time*  equipaMnt. 
memory,  interconnection  -  that  are  used  or  needed  in  the  algorithmic  solution 
of  given  problems.  This  dynamic  discipline.  concrete  eemputatiomal 
complexity,  is  important  for  both  hardware  sad  software  applieatioma.  and  it 
also  contributes  to  our  basic  understanding  of  computing.  It  is  only  natural 
that  the  great  technological  innovations  represented  by  Very-Large-Scale- 
Integrsted  (VLSI)  Circuitry  have  had  a  substantial  impact  on  this  field, 
opening  new  horizons  and  posing  challenging  problems.  Much  of  our  current 
research  (organized  below  in  five  subsections)  has  been  motivated  by  this 
revolution. 

12.2  Etaiipi  in  YLSI 

Parallel  computet  .on  provides  a  fundamental  answer  to  the  need  to 
increase  computer  processing  power,  and  VLSI  is  the  natural  habitat  for 
parallel  computation.  Ve  have  pursued  our  investigations  of  efficient  VLSI 
architectures,  suitable  for  a  variety  of  problems  (broad-purpose 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program  (U.S. 
Any,  U.S.  Navy,  U.S.  Air  Force)  under  contract  N00014-79-C-0424,  and  by  the 
National  Science  Foundation  under  grants  MCS  81-05552,  ECS  81-06939,  MCS  80- 
08845,  1ST  80-12240,  and  1ST  80-12242. 
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architectures) .  ¥e  hare  studied  is  detail  an  area-time  optiaal  nultiplier  for 
very  large  integers,  which  has  a  highly  regular  layout  in  the  fora  of  a  square 
aesh.  In  addition,  we  have  begun  the  exploration  of  a  novel  approach  to  the 
sorting  problea,  to  be  iapleaented  in  a  significant  variant  of  the  eube- 
conneoted-cycles,  and  that  offers  an  area-tiae  performance  superior  to 
previously  knows  scheaes. 

12.3  Layout  Theory  and  Desisn  Aids 

A  theory  of  layout  of  VLSI  circuits  is  a  current  technological  need, 
since  aost  plsceaent  sad  routing  tasks  are  presently  carried  out  using 
eapirical  techniques.  We  have  studied  the  channel  routing  problea  in  soae 
detail.  Ve  have,  determined  an  algor itha  for  two-layer  routing  which  bounds 
both  the  area  and  the  delay;  in  particular,  the  channel  width  and  the  aaxiaua 
wire  length  are,  in  soae  sense,  both  optimised.  We  also  developed  a  simple 
yet  effective  technique  for  calculating  a  lower  bound  on  the  number  of  tracks 
required  to  solve  a  given  two-layer  channel  routing  problem.  Moreover,  we 
have  shown  that  optimal  wiring  (i.e.,  density)  is  always  possible  if  just 
three  conducting  layers  are  available;  we  are  currently  trying  to  extend  this 
theory  to  more  complex  and  general  channel  routing  problems. 

In  addition,  we  have  investigated  design  problems  that  profitably  use  the 
techniques  of  computational  geometry,  in  particular  of  the  so-called  "geometry 
of  rectangles”.  In  this  area,  we  have  developed  a  search  technique  for  a 
collection  of  n  parallel  segments  (as  occur  in  VLSI  mask  technology) ,  which 
runs  in  time  0(log  n)  with  an  expected  storage  0(n) .  We  have  also  developed  a 
general  algorithm  for  the  computation  of  a  rectilinear  mask  expressed  as  an 
arbitrary  Boolean  function  of  a  set  of  rectilinear  masks  (papers  in 
preparation) . 

12.4  Pats  Structures 

Continuing  the  study  of  dynamic  eabeddings  between  data  structures,  we 
devised  an  optiaal  on-line  simulation  of  a  tree  by  a  d-diaensional  array  in 
tiae  0((N**( 1+1/d)) /(log  N)) .  Within  constant  multiplicative  factors,  no 
faster  simulation  is  possible.  The  simulation  is  cast  ,  terms  of  aultihead 
Turing  aaehines.  An  information-theoretic  arguaent  gives  the  worst-ease  lower 
bound  on  the  tiae  required. 
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We  have  also  bees  continuing  our  investigation  into  efficient  storage  and 
accessing  of  data  for  infomation-retrieval.  We  have  obtained  information- 
theoretic  lower  bounds  (and  shown  that  some  of  these  bounds  are  .  achievable) 
for  atorage  of  some  simple  data  structures  and  the  access  costs  for 
corresponding  operations. 

12.5  Foundations  si  Complexity  TfagfiLIY 

We  have  been  studying  the  properties  of  "context-free  graphs",  graphs  in 
which  most  of  the  important  graph  theoretic  properties  are  decidable.  Major 
results  were  obtained  relating  such  graphs  to  pushdown  automata  (papers  in 
preparation) . 

Also,  we  have  been  studying  how  alternating  automata  can  be  used  with 
languages  of  infinite  strings  and  trees.  This  subject  is  closely  related  to 
the  study  of  infinite  games.  The  use  of  such  automata  yields  a  new  and  much 
simpler  proof  of  Rabin's  theorem  concerning  the  decidability  of  second-order 
theories. 

12.6  Imaxe  Process  ins  Alaorlthms  ud.  Architectures 

We  have  studied  interleaved  pyramid  architectures  and  multiprocessor 
pyramids  for  bottom-up  image  analysis.  Current  research  focuses  on  achievable 
throughput  in  a  hierarchical  interconnection  of  processors  that  compute  non- 
oblivious  algorithms.  The  running  times  of  the  algorithms  are  modeled 
probabilistically. 
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13 .1  Nstnrsl  Lgggftjgl  Understsndlns* 

This  project  is  concerned  vith  judging  the  plausibility  of  natural 
language  sentence  Meanings*  vith  representing  Meaning  and  contest*  and  vith 
novel  processing  neehanisus  for  natural  language  understanding. 

Progress  vas  Made  on  the  understanding  and  representation  of  language 
describing  event  descriptions*  spatial  relationships,  tine,  and  causality. 
"Event  shape  diagrans,"  a  nev  representation  fon  for  dealing  vith  sequences 
and  concurrency  of  events*  vers  introduced  [13.23*13.24].  Ve  are  using  these 
and  other  seheuas  to  capture  and  represent  the  intricacies  of  fluid  flov* 
causality*  ueohanical  relationships,  and  other  properties.  To  aceosqtlish 
this*  it  is  iMportant  to  use  hierarchically  layered  Models,  representing  the 
eaMO  event  or  neehanisu  at  different  levels  of  abstraction.  In  addition, 

*  This  vork  vas  supported  in  part  by  the  Office  of  Naval  Eeseareh  under 
contract  N00014-75-C-0412  and  in£art  by  the  Air  Force  Office  of  Scientific 
Eeseareh  under  contract  F49620-82-K-0009 . 
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there  vast  be  separation  of  eoaplex  nodal a  at  a  given  level  into  siapler 
partial  aodels.  Mechanisas  for  representing  and  reasoning  abont  snob  systeas 
are  being  developed  as  an  extension  to  event  shape  diagrans. 

Event  shape  diagraas  were  nsed  to  show  soae  aeehanisas  for  deriving  the 
eorreet  aeanings  and  inferenoes  in  aetaphorieal  sentences*  sneh  as  "John  ate 
np  the  coapliaents,"  or  "John's  legs  ate  np  the  spaee  between  hiaself  and 
Mary."  [13.12]  Ve  have  investigated  a  feature  transfer  aodel  for 
nnderstanding  aetaphorieal  langnage.  This  nodal  aatehes  the  two  iteas  being 
eoapared  in  a  aetaphor*  and  transfers  aaterial  froa  the  "vehiele"  of  the 
aetaphor  to  the  target  doaain.  based  on  salieney  assessaents  of  various 
features . 

Ve  have  also  been  working  on  a  parallel*  analog  aodel  .  for  knowledge 
integration  and  decision-aaking  in  the  eontezt  of  natural  language  processing. 
The  aodel  involves  dynaaieally  eonstrueting  an  unstable  weighted  network  of 
possibilities*  while  concurrently  sifting  and  stabilising  the  network  such 
that  the  best  interpretation  is  highlighted.  The  aodel  has  been  suoeessfully 
applied  to  a  natural  language  processing  task  [13.21]. 

13.2  Scheaa  Acquisition* 

The  scheaa  learning  project  is  exploring  artificial  intelligence 
techniques  that  will  enable  a  eoaputer  systea  to  loan  general  world  knowledge 
in  the  fora  of  "seheaas"  through  its  interactions  with  an  external 
environaent.  A  scheaa  specifies  how  a  particular  real  world  situation  is 
likely  to  progress  and  why.  To  act  intelligently  in  assy  different  situations 
a  systea  requires  a  vast  noaber  of  seheaas.  Ike  problea  of  "prograaaing  in" 
large  nuabers  of  seheaas  has  greatly  Halted  the  intelligent  behavior  of  past 
AI  systeas. 

A  eoaputer  systea  ought  to  be  able  to  autcaatieally  acquire  these  seheaas 
directly  froa  its  experiences  nueh  as  people  do.  The  world  knowledge  already 
possessed  by  the  systea  is  used  to  generalise  individual  experiences  into 
seheaas.  Ve  have  devised  a  taxonoay  of  the  kinds  of  leaning  that  take  place* 
the  conditions  in  which  they  can  occur*  the  effects  eaeh  has  on  the  sys tea's 


*  This  work  was  supported  in  part  by  the  National  Science  Foundation  under 
contract  N8F  1ST  81-20234. 
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schema  memory,  and  ioat  of  the  mechanisms  by  vbieb  they  oan  occur.  The  kinds 
of  learning  are:  (1)  schema  composition,  (2)  secondary  effect  elevation,  (3) 
schema  alteration,  and  (4)  volitioaalization. 

A  computer  system  to  implement  and  test  these  ideas  is  near  completion 
[13.12,13.17.13.18], 

13.3  Coanitlve  Universale* 

In  order  to  write  natural  language  understanding  programs,  it  is 
essential  to  have  good  target  forms  for  representing  meaning.  In  this 
project,  we  have  combined  the  linguistics  and  artificial  intelligence 
approaehes,  in  the  hope  that  we  can  write  programs  that  are  better  models  for 
human  cognitive  processing,  as  well  as  better  computer  progrems  (that  is,  ones 
that  are  more  easily  extensible  and  leas  likely  to  require  radical  revision) . 

One  part  of  this  work  has  focused  on  roles  and  role  relationships  between 
the  real  world  elements  represented  by  the  noun  phrases  of  a  given  sentence , 
and  the  states  or  events  reflected  by  the  verb.  Ve  are  looking  at  role- 
related  syntactic,  semantic,  and  cognitive  phenomena  in  several  languages.  A 
questionnaire  was  prepared,  and  administered  to  speakers  of  Spanish,  Basque, 
and  Mandarin.  The  results  of  the  questionnaire  should  allow  us  to  assess  the 
variety  of  mechanisms  and  ease  role  distinctions  in  a  number  of  languages,  so 
that  we  can  better  build  proper  cognitive  distinctions  in  our  computer 
progrems. 

Ve  have  also  been  concerned  with  understanding  the  reasons  why  natural 
languages  differ  so  markedly  in  their  encoding  atrategiea  for  space-time 
meaning.  Ve  have  studied  this  phenomenon  in  English,  Jinghpaw  and  Burmese  of 
Sino-Tibetan  family  of  languages,  Blaekfoot,  Cree,  Delaware  and  Ojibwa  of  the 
Algonkian  family,  and  Tarascan. 

Finally,  we  are  investigating  the  use  of  machine-readable  dictionaries 
and  thesaurus  for  building  the  lexicon  for  a  natural  language  system. 


*  This  work  is  supported  by  the  National  Science  Foundation  under  contract 
NSF  1ST  81-17238. 
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13 .4  Hierarchical  Control  gad  Monitor  inn  with  gflB&fBlttll  Levels* 

Thin  work  denis  with  the  use  of  knowledge-base  architecture  and  planning 
control  aechanisas  to  perform  an  intelligent  aonitoring  task  in  a  coaplex 
doaain. 

The  doaain  for  onr  expert  planner  is  the  aircraft  flight  doaain.  The 
flight  doaain  knowledge  is  presently  organised  into  four  conceptnal  levels: 
the  route  level*  the  trajectory  level*  the  flight-control  level*  and  the 
aircraft  shbsystea  level.  The  levels  constitute  a  fora/ function  hierarchy 
which  facilitates  problea-solving  and  aonitoring.  As  an  exaaple  of  the 
fora/funetion  relationship*  an  aircraft  route  is  implemented  by  a  trajectory, 
which  is  in  turn  iapleaented  by  a  sequence  of  flight  control  settings. 

The  causal  fraaework  determines  the  conceptual  levels  and  the  planner  at 
each  level  only  has  to  consider  the  variables  within  the  causal  fraaework. 
The  levels  architecture  allows  us  to  focus  attention  and  also  provides  the 
ability  to  do  recovery  planning  based  on  a  change  in  the  context  at  any  level. 
This  property  is  unique  and  desirable  in  real-world  problea-solving. 

The  work  aeeoaplished  thus  far  consists  of  the  design  of  a  seaantieally 
oriented  planning  architecture  and  an  initial  design  of  the  inter-level 
planning-control  aeehanisa.  This  aechanisa  performs  aeta-plasning  in  the 
levels  context.  Probleas  and  issues  currently  under  investigation  inelude  (1) 
foeus  of  attention*  (2)  vertical  problea  decoaposition*  (3)  planning  control* 
(4)  inter-level  teleological  consistency*  (5)  doaain  knowledge  representation* 
aad  (6)  partial  planning. 

13.3  f  f  y  Mode  11ns  and  Antons  tic  piatnosis** 

Present  day  diagnosis  systeas  use  either  table-look-up  procedures  or 
production  rules  for  storing  the  information  needed.  Neither  approach  deals 
with  the  fact  that  for  a  diagnosis  systea  to  perfora  intelligently*  it  needs 
to  understand  the  basic  principles  which  are  eabodied  in  the  design  and 
structure  of  the  aechanisa.  Because  of  this  weakness  neither  approach  can 

*  This  work  is  supported  by  the  National  Aeronautics  aad  Space  Adaiai stra¬ 
ti  on  under  contract  NASA  NCCI-S2. 

*•  This  work  is  supported  in  part  by  the  National  Aeronautic  and  Space  Ad¬ 
ministration  under  contract  NASA  NCCI-52,  by  the  Departaent  of 
Trsasportatioa/Federsl  Aviation  Administration  under  contract  FA79WA-4360,  and 
in  part  by  the  United  Statee  Air  Foree  under  contract  F4M20-S2-I-0003. 
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kudli  ml  tipi*  fault*  *ff*etiY*ly. 

Hi*  purpose  of  this  work  is  two-fold.  First*  sa  will  develop  models  of 
uud*rstaudiu(  of  msehanisms,  about  th*ir  operating  priueiples*  structure— 
fuaotion  relationship  and  hierarchical  interaction.  Secondly*  w*  shall 
exploit  these  neehanisn  node Is  for  intelligently  diagnosing  aultiple  faults. 

As  a  typical  case  study  we  hare  developed  son*  initial  understanding 
aodels  for  a  generic  refrigeration  system  in  terms  of  such  concepts  as  heat- 
transfer  end  circulation.  With  the  help  of  these  concepts  models  w*  can 
analyse  any  schematic  of  a  refrigerator  and  identify  the  relevant  components 
which  arc  pertinent  to  an  operation*  such  as  circulation.  Ve  can  then  easily 
construct  the  eireulation  subsystem  as  an  abstraction. 

As  another  typical  ease  study  we  have  establiahed  a  systematic  process 
for  comprehending  the  working  of  an  amplifier.  The  essential  steps  are  (1) 
physical  instantiation*  (2)  semantic  template  verification*  and  (3)  constraint 
propagation.  With  these  steps  one  could  analyse  novel  new  structures  that 
fall  in  the  same  general  category  and  achieve  an  understanding  of  the 
funetion/struetur*  relationship  at  all  levels*  including  the  component  level. 
Ac  sub-structure  level  and  the  mechanism  level. 

Finally*  in  the  direction  of  diagnosis  some  additional  modeling  is 
needed.  Ve  develbped  flow-concept  models  for  initiation  isolation*  functional 
SMdels  for  detailed  analysis  and  constraint  models  for  verification.  These 
models  seem  to  be  necessary  for  a  theory  of  diagnosis  based  on  mechanism 
models. 

13 Air  Traffic  Control.  Problem- Solving.  gg£  Learning* 

•  The  purpose  of  air  traffic  control  is  two-fold.  First*  we  must  resolve 
all  potential  traffic  conflicts.  Secondly*  w*  must  deliver  all  aircraft 
safely  to  their  destinations*  in  a  fuel-efficient  manner,  with  as  little  delay 
as  possible. 

In  the  jargon  of  AI  w*  are  trying  to  do  problem-solving  with  conflicting 
goals.  Since  there  is  no  existing  theory  of  AFC  w*  have  the  problem  of 
generating  the  theory  as  well  as  carrying  out  its  computerisation.  Ve  have 
decided  it  would  be  a  good  idea  to  include  some  learning  capabilities  in  the 

*  This  work  is  supported  by  the  Department  of  Traasportation/Federsl  Avia¬ 
tion  Administration  under  contract  FA79VAr4360. 
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txp«rt  system  v<  in  building.  The  leaning  va  ara  attempting  is  eallad 
"advice-initiated  leaning."  Using  this  type  of  leaning*  the  system  first 
attempts  a  solution.  The  solution  is  revieved  by  an  adviser  who  presumably  is 
an  expert  in  ATC.  The  adviser  either  approves  the  solution  or  gives  adviee  in 
the  form  of  a  new  strategy.  The  computer  then  has  the  task  of  digasting  the 
advice;  in  particular*  it  must  figure  out  why  the  suggestion  is  indeed  better 
than  its  own  solution. 

13.7  fiflWtttfr  Vision* 

Hierarchical  organisations  of  small*  general  purpose  processing  elements 
to  perform  fast  bottom-up  image  analysis  wars  proposed.  Three  different 
architectures  were  given:  Pyramids*  interleaved  pyramids  and  pyramid  trees. 
Different  architectures  exhibit  different  performances  with  respect  to 
hardware  requirements  and  operating  speed  [13.16]. 

Approaches  to  dot  pattern  processing  using  Voronoi  neighborhoods  were 
proposed  [13.10].  Segmentations  of  dot  patterns  containing  homogeneously 
dense  segments  were  obtained.  Experiments  on  matching  dot  patterns  in  the 
presence  of  noise  were  performed. 

The  problem  of  generating  the  three-dimensional  occupancy  map  of  a  aeene 
from  stereo  image-pairs  was  investigated.  A  three  step  procedure  consisting 
of  obtaining  depth  map,  surface  map  and  a  volumetric  representation  was  given 
[13.15].  Algorithms  were  developed  to  update  the  octree  representation  of  a 
three-dimensional  object  as  the  objeot  undergoes  translation  and  rotation 
[13.27]. 

13.8  Iujjl  itfflltti  Processina  jM  Dynamic  Sm*  Analysis** 

The  processing  of  images  involving  motion  has  beeoaw  increasingly 
important.  Applications  include  target  tracking,  dynamic  robot  vision*  image 
bandwidth  compression,  and  medical  imaging.  Image  sequence  processing 
involves  a  large  amount  of  data.  However*  because  of  the  rapid  progress  in 


*  This  work  is  supported  by  the  National  Science  Foundation  under  grant  NSF 
IC8-81-OtfOOS. 

**  This  work  is  supported  in  part  by  the  National  Seienoe  Foundation  under 
great  NSF  ECS  Sl-12082,  and  in  part  by  the  Joint  Services  Electronics  Program 
(U.S.  Army,  U.S.  Navy,  and  U.S.  Air  Force)  under  Contract  N00014-79-C-0424. 
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coaputer,  LSI,  and  VLSI  technologies,  May  useful  processing  tusks  for  iaage 
sequences  osa  aow  bs  performed  is  a  reasonable  amount  of  tiae.  A  central 
issue  is  dynaaic  seeae  aaslysis  is  tke  deterainatioa  of  three-dimensional 
aotioa  sad  structure  froa  iaage  sequeaoes. 

Uaiqueaess  condi tioas  for  estiastiag  3-D  aotioa  psraaeters  froa  two 
fraaes  were  obtsiaed  [13.22].  Progress  was  aade  ia  deteraiaiag  aodels  for  3-D 
aotioa  estiastion  probleas  [13.23].  Ve  also  developed  efficient  algor ithas 
for  finding  corners  in  2-D  pictures  [13.20].  Applications  to  stereo  pair 
aatcbiag  aad  aotioa  estiastion  were  presented  [13.19].  Soae  techniques  for 
iaage  recognition  by  Batching  relational  structures  were  developed  in  [13.20]. 
A  simple  procedure  was  developed  [13.11]  to  generate  the  ainiana  depth  surfaee 
description  of  a  three-diaensional  scene  that  is  consistent  with  a  series  of 
iswges  froa  a  laterally  translating  caaera. 
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Faculty  Senior  Staff 

L.  V.  Lannaai  C.  G.  Robins  .  L.  C.  Smith 

S.  E.  Praaea  M.  B.  Villiaas 

fe«fa»ta  aiadlMl 

J.  Mendoza  C.  Tanopir  S.-C.  Vang 

14.1  Introduction 

During  the  1981-1982  time  period,  we  eondueted  in  information  retrieval  a 
number  of  research  and  devalopmont  projects  and  participated  in  the  University 
of  Illinois'  online  search  service.  Major  activities  included  analysis  of 
database  data;  a  teat  of  the  feasibility  of  automatically  determining  the 
overlap  between  bibliographic  databases}  development  of  a  computer-readable 
databases  directory;  and  participates  in  the  University  of  Illinois'  online 
search  service. 

1^*2  Analysis  of  Database  Data* 

Our  staff  maintains  s  body  of  information  about  commercially  available 
databases  called  the  Database  of  Dstabases.  We  have  continued  to  study 
various  statistical  eharaeteristics  of  this  population  based  on  the  material 
in  our  dstabase.  Such  statistics  sad  analyses  are  useful  both  to  researchers 
in  the  srea  and  to  the  users  and  producers  of  databases. 

Ac  particular  information  developed  in  the  studies  performed  in  the  last 
year  considered  the  age.  size,  type  (scientific,  medical,  etc.),  and  source 
(geveri*eat.  private,  etc.)  of  the  databases  covered.  At  present.  773 
Astabases  sre  covered  by  our  Database  of  Databases,  but  any  given  study  might 
include  only  a  specified  subset  of  those  depending  on  the  needs  of  the  study. 
The  most  resent  analyses  are  for  the  years  1980  and  1981. 


*  This  work  wss  supported  by  the  University  of  Illinois 


90 


14.  INFORMATION  RETRIEVAL 


The  software  used  in  this  research  consisted  of  e  pro gran  franc work 
eonteiniag  slots  into  which  logic  could  he  inserted  to  perform  specific  data 
selection*  cleanup*  end  analysis.  The  data  was  first  extracted  from  the 
database  by  our  date  access  prograa*  USER*  end  written  into  a  file  of  (tag* 
▼slue)  pairs.  The  analysis  programs  read  that. file  and  selected  the  tags  of 
interest.  An  effort  was  made  to  overlap  multiple  analyses  by  re-using  data 
seleeted  and  cleaned  for  prior  use. 

Cross  tabulations  were  made  based  on  defined  groupings  or  actual  values, 
depending  on  the  data  involved.  Thus,  cross  tabulations  were  done  involving 
both  raw  year  of  origin  versus  grouped  size  and  grouped  year  of  origin  versus 
grouped  size.  Statistical  summaries  included  means,  ranges,  and  standard 
deviations  of  raw  data  and  correlation  coefficients  between  raw  data  and 
grouped  data.  Detection  of  bad  data  included  identification  of  records 
containing  illegal  or  missing  values. 

The  work  involved  in  this  projeot  also  developed  a  framework  for  further 
studies  of  the  databases  available  to  researchers  today.  We  are  currently 
updating  the  database  and  improving  the  coverage  and  consistency  of  data 
fields  that  are  likely  to  be  important  in  future  studies. 

14.3  Federal  Emereencv  Manaaement  AfRRfiZ  PlttUlt  RffflgirWJtl  AlMlMMl 

Hid  Rf I9M61  PlWtQgy  Model* 

Ve  identified  and  analysed  word-oriented  (bibliographic*  textual, 
directory,  etc.)  databases  relevant  to  various  units  within  the  Federal 
Emergency  Management  Agency  (FEMA) . 

The  use  of  computer  databases  instead  of  amnual  searching  saves  time,  an 
iaqportant  factor  in  emergency  situations.  Computer  databases  usually  provide 
more  timely  information  and  allow  one-step  searching  of  complex  topics 
tailored  to  meet  individual  needs.  A  directory  indicating  which  databases 
contain  information  of  interest  to  FEMA  would  be  a  valuable  resource  for  all 
FEMA  units.  This  project  developed  s  methodology  for  compiling  such  a 
directory  and  also  provided  a  model  as  a  sample  directory. 

Databases  of  potential  relevance  to  FEMA  were  identified  and,  after 


*  This  work  was  supported  by  the  Federal  Emergency  Management  Agency  under 
grant  EMV-1-4058. 
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discussion  with  FEMA  personnel,  the  fifteen  judged  to  be  nost  relevant  to  the 
aost  FEMA  areas  of  interest  were  selected.  Subject  profiles  that  reflect  the 
intereats  of  all  individual  FEMA  units  were  developed  in  conjunction  with 
appropriate  FEMA  personnel.  Profiles  were  used  to  test  the  fifteen  databases 
online.  The  databases  were  ranked  by  subject  coverage,  aaount  of  inforaation, 
and  type  of  inforaation  based  on  the  online  search  results.  Te  created  a 
saaple  FEMA  Database  Directory  based  on  the  analysis  of  the  fifteen  databases 
and  including  saaples  of  citations  froa  each  database. 

PimlMY  al  Coanut  er-Ee  a  dab  1  e  Bibliographic  gltlfellM* 

The  Coaonter-Eeadable  Bibliographic  Mltllil  -  A  ViXM,t9tI  ULfl  Bill 
Sourcebook,  coapiled  and  edited  by  Professor  Martha  E.  Villiaas  and  Sandra 
Rouse,  was  published  in  1976  by  Knowledge  Industry  Publications,  Inc.  Updates 
to  the  Directory  were  issued  in  April,  1977  and  April,  1978.  A  second 
edition,  including  data  that  was  relevant  as  of  Deceaber,  1978,  was  published 
in  Autuan,  1979.  A  new  edition  appeared  ggrly  in  1982,  based  on  the  aost 
recent  updating  of  the  underlying  databaae.  The  current  edition  covers  773 
databases. 

fill*  CftlUctlgB 

1.  Database  processors  in  Europe  and  North  Anerica  were  contacted  by  nail 
for  inforaation  on  the  databases  they  provided  access  to. 

2.  The  questionnaire  designed  for  the  last  directory  was  revised  to  achieve 
a  more  consistent  foraat.  A  turnaround  questionnaire  was  coaputer- 
generated,  containing  the  inforaation  in  our  files  for  the  databases 
covered  by  the  last  directory.  Existing  data  entries  needed  only  to  be 
verified  or  changed  and  aissing  inforaation  supplied  by  the  database 
producer . 

3.  The  questionnairea  were  aailed  to  the  database  producers  and  U.  S.  and 
Canadian  producers  were  subsequently  contacted  by  telephone  for 
verification  of  the  data  and  aasistanee  with  questions  they  aight  have. 

*  This  work  was  partially  supported  by  Knowledge  Industry  Publications, 
Inc. 
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4.  bfMMtlM  oa  tli  iTiilililitj  of  mv  ditibuis  vos  obtolood  fron  til* 
ytbllilti  literature  sad  eoataets  with  dstsbss*  processor s/producers. 
Ho  producers  of  aov  databasos  vara  iatarriavad  first  by  pbona.  Tha 
eoapleted  fMStioaMir*  vas  tkam  Ballad  to  tbaa  for  varifieatioa  sad.  if 
Meossary,  follovod  ay  by  talaykoas. 

9*  latanad  fMstioaMires  vara  ebaekad  aad  data  aatrios  staadardixad. 

Data  Batw 

1.  Data  vas  aatarad  oalia*  iato  tbs  Databas*  of  Databasos  file  using  a 
preforaattiag  prograa  oa  aa  Alpha  Microsystaas  aierocoapotar.  Entries 
vara  tbaa  priated  oat.  verified,  aad  oorraoted  aad  edited  oaliae  as 
necessary  using  the  DBEDIT  prograa. 

Qraanixation 

Th*  Diraetory  is  organised  alphabetically  by  database  none  and  includes 
an  introduction  and  four  indexes:  Subject  Index,  Databas*  Naae/Acronya/- 
Synonya  Index,  Producer  Index  aad  Processor  Index. 

Inforaatioa  oa  each  databas*  follovs  th*  general  foraat: 

1.  Basic  Inforaation 

2.  Producer/Distributor/Generator  Inforaation 

3.  Availability  aad  Charges  for  Database  Tapes 

4.  Subject  Natter  and  Scope  of  Data  on  Tap* 

5.  Subject  Ana lysis/indexing  Data 

6.  Data  Eleaeats  Present  on  Tap* 

7.  Database  Services  Offered 


8.  User  Aids  Available 
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Production 

The  Directory  ves  formatted  for  printing  using  onr  software  end  printing 
■asters  produced  on  a  Xerox  9700  printer.  All  indexes  were  then  generated 
antoaatically  from  the  database  and  typeset  similarly. 

14.4  Bnlyprii.tY  lUiaail  2nlifl£  Search  Service 

The  Online  Search  Service  is  a  service  provided  for  students,  faculty, 
and  staff  of  the  University  by  the  Library  and  our  group.  Information 
retrieval  services  are  provided  by  accessing  computer-readable  databases  from 
the  following  search  systems:  Bibliographic  Retrieval  Services,  DIALOG 
Information  Systems,  New  York  Times  Information  Bank,  System  Development 
Corporation,  and  the  National  Library  of  Medicine. 

In  addition  to  serving  as  one  of  the  locations  where  searches  are 
conducted,  we  contribute  to  the  Service  by  providing  monthly  financial  and 
usage  reports  to  the  Library  administration  and  by  providing  expert  advice  in 
the  use  of  .information  retrieval  services.  The  financial  and  usage  reports 
are  generated  by  internally  produced  software  on  an  Alpha  Microsystems 
microcomputer  and  are  based  on  accounting  data  provided  by  the  online  vendors 
and  the  Library. 

14.5  MltUf  Analyses* 

We  are  conducting  ongoing  analyses  of  the  MARC  database.  This  work 
involves  processing  the  annual  files  distributed  by  the  Library  of  Congress 
and  generating  a  detailed  statistical  analysis  of  the  rate  of  occurrence  and 
sixe  distribution  of  each  field  of  the  records.  The  resulting  tables,  which 
have  appeared  in  the  Journal  of  Library  Automation,  show  for  each  field  (1) 
the  percentage  of  records  in  which  it  appears,  (2)  the  rate  of  occurrence 
within  records  where  it  appears,  and  (3)  the  length  of  the  field  in 
characters.  For  (2)  and  (3)  we  report  average,  minimum  and  maximum  values  and 
standard  deviations. 

The  distribution  of  records  in  Dewey  and  Library  of  Congress  clssses  is 
also  reports* -  along  with  average  record  lengths  for  each  class.  Summary 


*  This  work  is  supported  by  the  Council  on  Library  Resources. 
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statistics  src  also  listed  separately.  Finally,  ve  report  thi  changes  between 
reporting  periods  (shoving,  for  instance,  which  fields  are  used  sore  now  then 
at  the  lsst  report). 

Ve  have  also  begun  a  second  project  using  the  sene  computer  progress  to 
analyse  the  CONSER  database  for  the  Council  on  Library  Resources. 

This  data  sllovs  users  of  the  databases  to  optiaise  their  applications 
for  the  actual  characteristics  of  the  databases . 

14.6  Automatic  Detection  fll  Da tabs »e  Ovarian* 

This  research  prograa  addresses  the  problea  of  overlap  coverage  of 
journal  articles  aaong  aultiple  databases.  Overlap  leads  to  duplicative 
processing  of  aaterial  by  database  producers;  duplicative  processing  and 
storage  of  aaterial  by  online  and  batch  vendors  of  database  services;  and 
retrieval  of  duplicative  aaterial  by  users  of  databaaee.  The  net  result  is 
reflected  in  costs  associated  with  the  databases  and  an  unnecessary 
expenditure  of  tiae  on  the  part  of  users. 

The  true  extent  of  the  problea  is  unknown.  Several  studies  'have 
addressed  the  question  of  journal  coverage  overlaps,  soae  have  addressed  the 
question  of  article  overlaps  In  narrow  subject  areas,  and  one  has  attacked  the 
problea  of  identifying  aaxiaua  possible  overlaps  aaong  specific  dstsbases. 
None  has  studied  the  problea  of  deteraining  the  actual  overlaps  aaong  database 
services. 

Ve  sre  developing  a  aethodology  for  autoaatically  deteraining  overlap 
aaong  databases  and  testing  and  iapleaenting  this  aethodology  on  nine 
dstsbsses:  AGRICOLA,  BIOSIS  Previews.  CACon.  Coaaonweslth  Agricultural  Bureaux 
(CAB).  Engineering  Index.  Excerpts  Medics,  Inspec.  Medline,  end  Science 
Citation  Index  (SCI).  The  aethodology  involves  snalysis  of  the  databsses  to 
deteraine  the  data  eleaents  present  in  end  ooaaon  to  aultiple  dstsbsses  to 
deterains  effective  keys  for  astehing  records  in  aultiple  files.  The  aatehing 
process  will  be  done  in  st  lesst  two  pssses.  The  first  psss.  or  brosd  screen, 
will  use  s  natch  key  with  high  diseriainstion  power  to  identify  esndidate 
duplicste  records.  A  set  of  finer  screens  will  be  used  for  verifying  pairs  of 


*  This  work  is  supported  by  the  Nstional  Science  Foundstion  under  grsnt  NSF 
1ST  79-21018. 
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records  as  duplicates.  The  second  pass  s»y  enploy  Harrison  keys  (which  are 
bit  string  representations  of  titles)  and  Hanning  distance  techniques  (for 
owereoning  slight  variations  in  representation  of  titles)  together  with 
natehes  on  other  elenents  as  needed. 

Results  will  be  tested:  (1)  by  using  nanual  verification  of  a 
statistically  representative  sanple  in  order  to  detemine  nisnatches.  and  (2) 
by  using  a  known  set  of  duplicate  records  to  see  if  they  are  identified  by  the 
screens.  The  evaluation  in  terns  of  nisnatches  (wrong  natehes)  and  nissed 
natehes  will  provide  a  neasure  of  the  effectiveness  of  the  techniques. 
Further  evaluation  in  terns  of  gross  costs  will  be  carried  out  in  order  to 
deternine  the  efficiency  of  the  techniques. 

The  study  will  produce  overlap  statistics  anong  nine  databases,  data 
elenent  frequency  statistics,  and  a  nethodology  for  duplicate  detection  that 
night  be  used  by  database  producers  or  database  vendors.  The  resultant 
nethodology  could  be  used  as  a  preprocessing  step  to  reduce  storage  and 
processing  costs  or  as  a  post-retrieval  step  to  elininate  the  duplicative 
output  delivered  to  users. 

During  the  second  year  of  the  project,  tapes  have  continued  to  arrive 
fron  the  database  producers,  covering  the  1978-1979  period.  Sanple  subsets 
have  been  read  and  prelininary  planning  for  the  file  structures  and  processing 
flow  has  been  done.  The  contents  of  the  sanples  have  been  studied 
statistically  to  identify  data  elenents  useful  in  the  natehing  process  and  to 
allow  estination  of  storage  requirenents.  The  software  has  been  ported  fron 
the  DEC  Systen/10  to  the  DEC  VAX  11/780  and  further  sanples  have  been  read  and 
evaluated  preparatory  to  perforuing  large-scale  conversions. 
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15 .1  IfaltlMli-TmiBil  fiililAl  Coaannications* 

This  ia  a  aajor  research  area  which  i sc lodes  various  prohleaia  is  spread- 
speetrua  coaaunications  sad  rsadoa-aeeess  coaenmieatioas  networks.  The  topies 
that  were  investifated  during  the  past  year  inelude  the  perforaaaee  of 
spread-speetrua  aultiple-aocess  coaaunications  systeas,  the  effects  of  fading 
on  spread-speetrua  coaaunications*  coaaonications  ia  the  presence  of  jaaaiag, 
error-control  coding  for  aultiple-aocess  systeas*  conflict  resolution  ia 
raadoa-aeeess  coaaunications  systeas*  packet  radio  coaaunications  systeas*  and 
routing  algorithas  for  eoBaoaications  networks.  Progress  in  this  research  is 
described  in  the  subsections  that  follow. 


*  This  research  was  supported  by  the  A ray  Research  Office  under  contract 
DAA629-81-K-0064 ,  by  the  Joint  Services  Electronics  Prograa  (U.S.  Aray*  U.8. 
Navy*  and  U.S.  Air  Force)  under  contract  N00014-79-C-0424,  by  the  Naval 
Research  Laboratory  under  contract  N00014-80-C-0802,  and  by  the  Defense  Ad¬ 
vanced  Research  Projects  Agency  under  contract  IDA  903-82-C-0026 . 
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15.1.1  Spread-Soectran  Multiple-Access  Conmua lest ions  Systems 

fa  havs  completed  tha  development  of  bounds  and  approximation  a  for  tba 
average  probability  of  arror  in  direct-sequence  (DS)  spread-spectrum 
multiple-access  (SSMA)  comannieationa  systems,  and  tba  raanlta  have-  baaa 
pnbliabad  in  tba  TERR  Transactions  on  Commanieationa  spacial  is ana  on  apraad- 
spactrna  eonnnnicationa.  Tba  bounds,  vbicb  ara  prasantad  in  [15.30]  ara  bast 
snitad  for  SSMA  systans  with  fay  simaltaneously-active  transaittara.  Sinca 
both  nppar  and  lowar  bounds  ara  derived,  approximations  ean  also  ba  obtained 
from  tbasa  rasnlts.  For  SSMA  aystams  vitb  a  large  nambar  of  simultaneous 
transmissions,  tba  bounds  ara  difficult  to  ayaluata,  but  tba  approximation 
dayalopad  in  [15.9]  can  ba  applied.  This  approximation  is  based  on  tba 
integration  of  tba  characteristic  function  of  tba  multiple-access 
interference,  and  its  computational  requirements  increase  only  linearly  in  E, 
tba  number  of  simultaneous  transmissions.  The  methods  dayalopad  in  tbasa  two 
papers  baya  proved  to  ba  extremely  useful  in  tba  study  of  DS/SSMA  and  hybrid 
SSMA  systems. 

An  investigation  of  tba  affects  of  partial  correlation  on  tba  performance 
of  DS/SSMA  systems  vas  also  completed  during  tba  reporting  period.  Soma  of 
tba  results  of  this  study  ara  prasautad  in  [15.58].  Our  investigation  has 
primarily  focused  on  systems  for  vbiob  tba  number  of  chips  par  bit,  N,  and  tba 
period  of  tba  signature  sequences,  p,  ara  relatively  prime  numbers  (this  is 
usually  tba  case  in  practical  systems  if  p  >  N) .  In  this  ease  we  find  that 
tba  resulting  loss  of  control  of  tba  sequence  phases  leads  to  an  increase  in 
tba  multiple-access  interference  over  that  present  in  a  DS/SSMA  system  which 
baa  one  data  bit  par  period  (i.e.,  p  -  N)  and  employs  the  signature  sequences 
in  their  optimal  phases. 

Ve  have  considered  the  design  of  qusdriphase  sequences  vitb  good 
correlation  properties.  Ve  are  able  to  construct  a  set  of  N  +  1  sequences  of 
period  N  (where  N  is  a  multiple  of  4  and  N  +  1  is  a  prime  power),  for  which 
the  correlation  parameter  ia  upper  bounded  by  3(N+1)~*^2*  Ve  can  also 
enlarge  this  set  to  have  2(N+1)  sequences  with  an  upper  bound  on  of 
3(N+1)~*/^  +  5.  Experimental  results  show  that  the  bound  is  actually  quite 
weak. 

An  investigation  of  the  performance  i^rovements  achievable  by  use  of 
error-control  coding  has  also  been  initiated,  and  preliminary  results  are 
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presented  in  several  reeent  papers  (e.g.,  [15.36*15.59,15.61]).  The  general 
eonelnsions  are  tkat  error-control  coding  is  extreaely  beneficial  in  nearly 
all  coaannieation  environaents,  and  it  is  absolutely 'necessary  for  acceptable 
bit  error  rates  in  frequency- bopped  spread-spectra  ooaannieations  in  tbe 
presence  of  jaaaing,  fading,  or  aultiple-aeeess  interference.  Coding  gains  of 
20  to  30  dB  are  eoaaon  in  snob  systeas. 

15 .1 .2  Soread-Soectraa  Coaannlcations  for  Packet  Radio  Networks 

Becanse  of  tbe  severe  specnlsr-anltipstb  fading  that  arises  in  urban 
areas,  ve  bsve  investigated  tbe  effects  of  suck  fading  on  tbe  perforaanee  of 
direct-sequence  spread-spec trua  eoaaunieations.  Results  for  both  coberent 
[15.46]  end  noncoberent  [15.47]  receivers  bsve  been  published.  Although  these 
results  are  for  receivers  which  utilise  standard  aatobed  filters,  work  is  also 
in  progress  to  evaluate  tbe  perforaanee  of  direct-sequence  spread-spec trua 
coaauniestions  vis  specular  aultipatb  channels  when  post-detection  integration 
or  RAKE- type  filters  are  eaployed  in  tbe  receiver.  Of  particular  interest  is 
tbe  perforaanee  of  receivers  with  post-detection  filtering  when  eaployed  on 
channels  with  jaaaing  and  aultiple-aeeess  interference. 

In  addition  to  tbe  research  on  direct-sequence  spread  speetrua,  we  are 
also  pursuing  s  parallel  investigation  of  frequency-bopped  spread  speetrua  for 
packet  radio  networks.  It  appears  that  for  channels  with  partial-band  jaaaing 
or  aultiple-aeeess  interference,  frequency-bopped  aodulation  in  conjunction 
with  error-control  coding  provides  a  very  robust  systea.  It  seeas  at  this 
point  that  tbe  tradeoff  is  between  tbe  snti-aultipatb  capability  of  direct- 
sequence  spread  speetrua  and  tbe  snti-jaa  sad  aultiple-aeeess  capabilities  of 
coded,  frequency-bopped  spread  speetrua.  Of  course,  error-control  coding  is 
also  useful  in  s  direct-sequence  systea,  but  for  asny  types  of  interference 
(e.g.,  aultiple-aeeess  interference)  the  aaouat  of  interleaving  required  in 
direct-sequence  systeas  aay  be  prohibitive.  Furtheraore,  coding  cannot  take 
care  of  tbe  problea  of  unequal  power  levels  (i.e.,  tbe  so-called  "near-far 
problea")  in  a  packet  radio  network. 

15.1.3  Cod  ins.  Modulation.  a&4  Routine  Xflt  Coaaunicstion  Hf.tl.9lfa 

Ibis  research  area  inoludes  several  related  topics  that  arise  in 
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ooHUieitira  networks  wkisk  Mat  operate  in  tka  presence  of  j earning ,  fading, 
and  multiple-access  interference.  Ike  topics  inclnde  f r eqne  uey-hopped 
spread-spectrum  communication  wia  fading  ekanneia,  communication  in  a  koatile 
j  caning  environment,  and  ronting  in  communications  networks.  Progress  in 
research  on  these  topics  is  summarized  in  the  three  paragraphs  that  follow. 

Our  work  on  slow-frequeaey-hopped  spread-spec traa  communications  systeas 
is  concerned  with  the  probability  of  error  for  sneh  systeas  when  operating  in 
the  presence  of  fading,  multiple-access  interference,  j caning,  and  other 
radio-frequency  interference  (RFI) .  Various  types  of  fading  hare  been 
considered  including  selective  and  nonselect ive  Rayleigh  and  Rician  fading. 
Ve  have  also  investigated  the  performance  of  error-control  coding 
(particularly  burst-error-eorreeting  coding)  for  frequency-hopping  spread- 
spectrum  systeas.  One  of  the  most  premising  coding  schemes  that  we  have 
investigated  utilizes  Reed-Soloaon  codes  with  partial  interleaving  for  slow- 
frequency-hopping  (more  than  one  eede  symbol  per  hop)  or  Reed-Solomon  codes 
and  diversity  for  fast-frequency-hopping  (one  code  symbol  per  hop).  Results 
are  published  in  several  recent  papers  including  [15.9,15.41,15.44, 
15.59.15.40]. 

In  related  work,  we  have  derived  the  optimum  receiver  for  a  fast- 
frequency-hopping  scheme  proposed  for  mobile  radio  applications,  and  obtained 
bounds  on  the  error  probability  for  a  suboptimum  receiver  which  is  easier  to 
implement  than  the  optimum  receiver.  These  bounds  have  been  compared  with 
simulation  results  on  error  probsbility  for  the  originally-proposed  suboptimum 
receiver.  It  is  found  thst  the  new  receiver  structure  is  superior  to  the 
original  one. 

In  the  area  of  communication  in  the  presence  of  jamming,  our  research  has 
focused  on  information-theoretic  fundamentals.  Vs  formulated  a  two-person, 
zero  sum  game  with  mutual  information  as  the  payoff  function.  The  players  in 
this  gsme  sre  the  jsmmer  and  the  coder.  By  solving  the  game,  one  first  finds 
ssddlepoint  strstegies  which  indicate  the  general  kind  of  jsmming/anti-jamming 
strategies  the  resl  plsyers  ought  to  use  in  practice.  The  work  has  continued 
and  several  saddlepoiats  have  been  calculated.  Also,  the  performsnee  of 
seversl  explicit  coding  schemes  hsve  been  evslusted  in  the  presence  of  optimal 
jammers.  In  a  related  area  we  undertook  a  serious  study  of  the  relstive 
merits  of  interlssved  vs.  un-interlesved  coding  schemes  in  the  presence  of 
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liuril  latarfuiM*  phtiaim,  la«ii4ia|  jamming,  fading  and  other  KFI. 

A  aaia  raanlt  of  oar  research  la  communication  network  analysis  daring 
this  period  has  been  to  dbtain  beanda  on  the  delay  inearred  in  a  baf fared 
station  with  an  arbitrary  packet  arriwal  process.  This  bound  applies  to 
stations  within  a  network  aader  dynamic,  non-feedbaok  roatiag  strategies,  i 
rather  sharp  boaad  was  also  obtained  for  the  delay  inearred  in  a  station  when 
the  treffie  is  obtained  from  a  generalised  roand  robin  fanoat  of  a  Poisson 
traffic  stream. 

Recently  we  hare  also  obtained  a  new  eharaeterisation  of  optiMl  dynamic 
roating  strategies  for  interacting  queues  with  state  feedbaek  information.  In 
the  oomisg  year  we  eapeet  to  complete  a  naif  led  performance  analysis  of 
optimal  and  efficient  saboptiMl  dynamic  roatiag  strategies  for  packet 
communication  networks  both  with  and  without  state  feedbaek  intonation. 


15.1.4  Random  Aocoaa  Techniques  for  Conmanication  Networks 
Xilfe  Frequency  Hoooins 

This  research  is  primarily  concerned  with  basic  random  access  teehmiqaes 
which  are  appropriate  for  mobile  packet  radio  networks  with  frequency  hopping. 
Two  key  features  of  such  networks  which  impact  random  access  protocols  are 
first  that  channel  feedback  is  serereiy  limited  and  second  that  frequency 
modulation  in  itself  yields  a  multiple  access  capability  which  should  be 
exploited.  We  hare  introduced  new  dynamic  retransmission  control  schemes  and 
haws  shown  their  suitability  in  a  frequency  hopping  ear iroament.  This  work 
included  identifying  possible  feedbaek  inforMtion  and  determining  how  it 
could  be  used  to  achieve  stable  throughput  in  spite  of  a  fluctuating  demand. 
Other  suitable  random  access  strategies  which  require  only  acknowledgement 
feedbaek  have  also  been  studied. 

More  realistic  models  of  stations  which  include  time-varying  packet 
generation  rates  have  been  considered  in  conjunction  with  random  access 
schemes.  So  far  exact  analysis  of  H>NA  and  perfect  scheduling  protocols  have 
been  achieved  using  matrix  iteration  methods. 

Certain  aspects  of  mobile  network  organisation  have  also  been  considered. 
Ve  have  discovered  how  differential  transmission  radii  in  different  frequency 
bands  can  be  exploited  to  increase  network  throughput  and  decrease  multiaccess 
interference.  Ve  have  proposed  and  analysed  methods  whereby  traffic  oan  be 
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efficiently  distributed  uoi|  available  time  or  frequency  slots  in  s 
decentralised  way  to  achieve  balanced  loads  with  reduced  traffic  variance. 
Finally,  centralised  and  decentralised  algorithms  have  been  obtained  for 
open-loop  optimal  dynamic  routing  in  communication  networks  modeled  by  pure 
flow  equations.  The  algorithms  appear  to  be  feasible  to  iaplement  in  much 
larger  networka  than  previously  known  optimal  dynamic  routing  algorithms. 

1S.2  dismal  Detection.  iaA  Filterins* 

The  overall  purpose  of  this  research  effort  is  to  develop  effective 
statistieal  signal  processing  procedures  for  applications  in  communications 
and  control.  Specific  areas  currently  under  investigation  include  the 
development  of  robust  detection,  estimation  and  filtering  procedures  for 
operation  in  uncertain  statistical  environments,  the  study  of  techniques  for 
effieient  digitisation  of  signal  detection  systems,  and  the  analysis  of 
multistage  decision  processes.  Progress  in  these  areas  is  described  briefly 
in  the  f cl lowing  subsections. 

19.2.1  debus t  liana  1  Procoaaint  ift  ffSlfltSiB  glllillM  BRTir.Wt»tt 

In  general,  the  area  of  robustness  deals  with  the  design  of  systems  and 
procedures  which  are  relatively  insensitive  (in  terms  of  performance)  to  small 
deviatiems  from  an  assumed  model.  Robust  techniques  are  of  interest  in  a  wide 
variety  of  applications  including  rsdar,  sonar,  and  seismology,  inasmuch  as 
inaccurate  models  are  the  norm  for  such  situations;  and  robustness  has  been 
studied  widely  in  the  contests  of  communications,  control,  and  inferential 
statistics.  The  basic  assumption  for  the  snalyticsl  study  of  robustness  in 
detection  and  filtering  is  that  the  statistics  of  signals  and/or  noise  are  not 
known  esactly,  but  rather  are  known  to  be  within  some  (usually  nonparametrie) 
classes  representing  uncertainty  in  the  underlying  model.  The  primary 
technique  for  designing  a  robust  system  is  to  seek  s  system  schieving  the  best 
worst-case  performance  over  the  relevant  uncertainty  classes;  i.e.,  the 
primary  design  philosophy  is  ninimax.  Within  this  context  several  detection. 


*  This  research  was  supported  by  the  Army  Research  Office  under  contract 
DAA629-81-K-0062 ,  the  Joint  Services  Electronics  Program  (U.S.  Army,  U.S. 
Navy,  and  U.S.  Air  Force)  under  contract  N00014-79-C-0424,  the  Office  of  Naval 
Research  under  contraot  N00014-S1-K-0014,  and  the  National  Science  Foundation 
under  grant  ECS-79-164S3. 


•itlutioa  *  ad  filtering  problems  have  been  eoasidered  under  this  project*  and 
these  are  described  ia  the  following  paragraphs. 

One  area  ia  vhieh  several  ainiaax  and  robustness  results  have  beea 
obtained  is  that  of  state  estimation  sad  eoatrol  for  linear  dynamical  systems 
with  uncertain  statistieal  behavior.  Ia  particular*  minimax  results  have  beea 
developed  for  steady-state  liaear  estimation  [15.29]  sad  eoatrol  [15.20,15.53* 
15.54]  of  time- invariant  systems  and  for  finite-interval  estimation  and 
regulations  [15.53*15.67]  of  time-varying  systems.  Also,  an  algorithm  for  the 
state  estimation  of  sueh  systems  with  random  jump  inputs  has  beea  developed 
and  is  described  in  [15.4].  Along  these  same  lines  the  problem  of  robust 
filtering  and  smoothing  of  stationary  proeesses  with  uaeertain  speetra  has 
also  been  considered.  Specific  results  inelude  a  characterisation  of  miaimsx 
smoothers  for  homogeneous  random  fields  with  spectral  measures  known  only  to 
lie  in  general  uncertainty  classes  geaeri  *&  by  Choquet  capacities  [15.24],  an 
analysis  of  the  effects  on  performance  of  spectral  uncertainty  on  stationary 
linear  filtering  [15.34],  and  a  study  of  minimax  filtering  for  observed 
Poisson  processes  with  uncertain  rate  characteristics  [15.44,15.45].  Two 
related  studies  include  a  characterisation  of  the  MSB  rate-distortion  function 
for  information  sources  with  uncertain  spectra  lying  in  capacity  classes  [P10] 
and  the  development  of  jam/ antijam  strategies  for  data  transmission  through 
the  Gaussian  test  channel  with  an  intelligent  jammer  [15.5,15.37]. 

Another  aspect  of  robustness  in  which  several  results  have  been  obtained 
recently  is  that  of  robust  signal  detection.  Specific  results  ia  this  area 
inelude  the  development  of  robust  detectors  for  the  (nonlinear)  detection  of 
signals  in  additive  dependent  noise  with  uncertain  first-order  distribution 
[15.2S],  a  study  of  the  discrimination  of  point  processes  with  uaeertain  rate 
functions  [15.45],  the  design  of  matched  filters  for  digital  transmission 
through  uncertain  channels  [15.36,15.65,15.67],  and  the  characterisation  of 
robustness  in  the  quadratic  receiver  design  problem  [15.64,15.67].  Also, 
[15.19]  considers  a  design  technique  known  as  the  segment  method  as  an 
alternative  to  miaimax  for  designing  binary  decision  systems  within 
uncertainty. 

In  addition  to  the  above  specific  results,  a  general  study  of  the 
characterisation  of  solutions  to  miaimax  rbbuat  design  problems  is  reported  ia 
[15,64,15,67].  The  techniques  developed  hero  allow  the  unification  and 
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general ization  of  a any  previoualy  studied  robustness  analyses  as  well  the 
extension  to  new  design  problems.  Also  during  this  reporting  period  several 
survey  papers  on  robustness  topics  bsve  been  prepared  [15.18*15.51*15.57]. 
The  topies  of  these  surveys  include  robust  signal  processing  for  communication 
systems  [15.18]*  robust  filtering  [15.51]*  and  the  use  of  distance  measures  in 
robust  detection  and  discrimination  [15.57]. 

15.2.2  Pialtal  lM  Multistats  Signal  Detection  SllUli 

In  addition  to  the  robustness  work  surveyed  in  the  above  subsection* 
several  other  aspects  of  the  design  and  analysis  of  signal  detection  systems 
have  been  considered  recently.  Progress  in  this  area  is  described  briefly  in 
the  following  paragraphs. 

One  area  in  which  recent  progress  has  been  made  is  that  of  designing 
efficient  systems  for  detecting  signals  in  dependent  non-Gaussian  noise 
backgrounds.  For  example*  optimum  detection  of  signals  in  weakly  dependent 
noise  is  considered  in  [15.27].  It  is  shown  here  that*  for  a  moving-average 
dependence  model,  efficient  detection  can  be  achieved  by  reshaping  the 
independent-noise  optimum  detection  nonlinearity  with  an  additive  linear 
correction  term.  An  investigation  of  several  commonly  used  noise  models 
indicates  thst  the  resulting  performence  improvement  is  most  significant  for 
impulsive  types  of  noise.  In  [15.17]  similar  modified  structures  are 
considered  for  detecting  signals  in  $ -nixing  noise.  It  is  demonstrated  that 
the  ^-mixing  model  with  only  the  mixing  coefficients  given  is  not  sufficiently 
descriptive  of  dependence  structure  to  admit  a  design  which  is  uniformly 
better  than  the  independent-noise  design.  (This  is  in  contrast  to  the  moving- 
average  model  of  [15.27].)  However*  optimum  memory less  designs  for  4 -mixing 
and  other  dependence  models  osn  be  derived  using  general  results  of  [15.25] 
provided  that  the  second-order  distributions  of  the  noise-process  are  known. 
The  optimum  system  in  this  case  is  the  solution  to  s  Fredholm  operator 
equation  of  the  second  kind  and,  as  such*  can  be  characterized  in  terms  of  the 
Barrett-Lampsrd  expsnsion  of  the  noise  process. 

A  second  area  of  recent  interest  is  that  of  anslyzing  the  efficiencies  of 
multistage  detection  procedures.  Two  basio  types  of  procedures  hsve  been 
studied  in  this  context.  In  [15.(6],  a  useful  class  of  multistage  binary 
hypothesis  tests  is  proposed  and  analyzed.  These  tests  operate  by  testing  an 
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initial  aaapla  and  than  taking  additional  saaplas  if  tk#  first  saaple 
indicates  rajaetion  of  a  partionlar  ona  of  tka  hypo tka ses.  It  it  aaan  tkat 
significant  afficiancy  ean  ba  gainad  (over  f ixed-sample-size  procadnras)  by 
doing  tkis  if  ona  of  tka  kypotkasas  is  rara  (as  in  tka  easa,  for  example*  in 
saarek  radar) .  A  noaparaaatric  version  of  tkis  proeadnra  is  disenssad  in 
[15.331*  In  (15.321*  tka  (Pita an)  affieianeias  of  tmneatad  Vald  sequential 
tacts  ara  analyzad.  Thasa  taste*  altkongk  known  for  sons  tine,  kave  not  boon 
analyzed  prarionsly  in  tkis  Banner.  Tkis  analysis  elaarly  indicates  tka 
faworability  of  traneation  in  snek  situations,  particularly  wken  wary  saall 
error  probabilities  ara  desired. 

Finally,  soaa  related  work  includes  a  numerical  study  of  tka  convergence 
bekavior  of  tka  relative  affieianeias  of  some  eonaonly  used  detection  systems 
[15.231*  and  a  correct  derivation  of  tka  outpnt  signal- to-noise  ratio  for 
quadratic  raeaivars  and  fading  ekannels  [15.351  wkiek  previously  ksd  been 
incorrectly  derived  in  tka  literature. 
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16.1  fllyrftlgfciftll  Design  Ai4t  lOX  XLSl  Circuits* 

This  work  covers  two  related  topics  in  the  area  of  design  verification 
for  VLSI  circuits.  One  topic  is  concerned  with  logic  expression  extraction 
from  the  circuit  layout  or  the  interconnection  description,  and  the  other 
topic  is  concerned  with  the  development  of  fast  and  relatively  accurate 
methods  fox  timing  simulation. 

In  the  design  verification  process  of  VLSI  circuits  it  is  often 
desirable,  or  even  necessary,  to  obtain  symbolic  Boolean  expressions  at 
specified  nodes  in  the  system  from  the  extracted  circuit  description.  The 
extracted  logic  description  can  then  be  compared  with  the  original  logio 
specifications  to  verify  the  design.  This  verification,  in  fact,  should  be 
carried  out  before  any  circuit  or  timing  simulation  is  performed.  Extracted 
logic  expressions  earn  also  be  used  to  desoribe  parts  of  the  system  for  mixed¬ 
mode  simulation  [II]  or  eould  be  evaluated  using  ternary  algebra  for  switch- 
level  [R2]  or  multiple-delay  simulation  [13] . 

In  our  research  work,  an  approach  for  generating  logic  expressions  at 


This  work  was  supported  by  the  Joint  Services  Electronics  Program  (U.S. 
Army,  U.S.  Navy,  U.S.  Air  Force)  under  Contract  N00014-79-C-0424. 
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specified  nodes  in  NOS  circuits  has  been  developed.  The  approach  is  based  on 
partitioning  the  circuit  into  "one-way"  subcircuits  and  then  using  path 
finding  algorithsis  to  generate  logic  expressions  at  the  subcircuit  level. 
Composite  expressions  are  then  generated  by  interconnecting  the  subcircuits. 
The  approach  has  been  successfully  implemented  in  a  computer  program  which  can 
generate  logic  expressions  of  extracted  or  specified  circuit  descriptions.  In 
addition,  the  progrsm  is  capable  of  performing  switch-level  simulation  for 
given  input  logic  sequence  specifications. 

In  the  area  of  timing  simlation  a  new  approach  has  been  developed.  The 
approach  is  based  on  a  modified  Gauss-Seidel-Newton  method  for  solving  the 
partitioned  circuit  equations.  Then  feedback  paths  exist  in  the  circuit,  a 
predictor  formula  is  used  to  estimate  the  state  of  the  "yet  unsolved" 
subcircuits.  It  has  been  found  that  the  new  approach  gives  more  accurate 
results  than  existing  Gauss-Seidel  methods,  such  as  the  one  used  in  MDTIS-C 
[R4] ,  without  requiring  much  additional  computation.  The  approach  has  been 
implemented  and  tested  in  a  computer  progrsm  for  simulating  NOS  circuits.  The 
program  also  includes  procedures  for  analysis  sequencing,  feedback  detection, 
and  an  input/output  subcircuit  scheduler  which  schedules  only  those 
subcircuits  that  affect  the  output  directly  or  indirectly. 

In  order  to  further  improve  tne  speed  of  timing  simulation,  a  new 
approach  to  multiple-delay  logic  simulation  is  under  investigation.  Cirouit 
simulation  is  used  to  determine  the  propagation  delay  as  a  function  of  the 
rise  and  fall  times  of  the  input  signal  and  the  loading  on  the  eirenit  for  a 
given  technology.  This  information  is  stored  in  a  table.  Sealing  techniques 
can  be  used  to  compensate  for  different  sixed  transistors.  Simulation  results 
from  combinatorial  circuits  indicate  a  close  match  between  the  signal 
waveforms  obtained  from  cirouit  simulation  with  the  program  SPICE  and  those 
obtained  from  a  three  level  multiple-delay  logic  simulator  which  uses  the 
above  table  look-up  and  scaling  techniques  [16.26]. 

16.2  Swltched-Cauac itor  Filters  -  gtrifitMt*  AAd  gtiUitifiil  fitlilH* 

Research  on  switched-eapacitor  (SC)  filter  design  techniques  has  resulted 


•This  work  was  supported  by  the  National  Science  Foundation  under  Grant  ENG 
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in  new  SC  structures  which  can  be  implemented  basically  in  two  different  ways, 
either  by  using  SC  first-order  building  blocks  or  SC  biquads.  These 
structures  are  the  state-space  (SS)  [K5-R9,16.21,16 .27] ,  the  mul tipi e-input 
follow-the-leader  feedback  (MIFLF)  [R10.R11],  and  the  inverse  follow-the- 
leader  feedback  (IFLF)  [16.10,16.11].  Theory  and  procedures  for  the  optimum 
implementation  of  these  types  of  SC  filters  have  been  developed. 

The  sensitivity  of  the  transfer  function  of  the  IFLF  circuit  with  respect 
to  the  capacitance  ratios  is  surprisingly  low  as  compared  to  other  low- 
sensitivity  SC  structures,  such  as  the  bilinear  ladder  [K11.K12].  As 
compared  to  the  HlFLF,  the  IFLF  has  lower  passive-sensitivity  and  higher 
active-sensitivity.  The  MIFLF  active-sensitivity  is  also  lower  than  that  of 
the  bilinear  ladder  [R12] . 

In  the  area  of  statistical  circuit  analysis  and  design,  several  variance 
reduction  techniques  have  been  studied  with  respect  to  their  application  in 
estiauting  circuit  yields.  It  has  been  shown  that  given  an  approximation  of 
the  region  of  acceptability  with  a  given  accuracy,  importance  sampling 
techniques  deteriorate  rapidly  as  the  dimensionality  of  the  parameter  space 
increases.  However,  modest  reductions  in  the  variance  of  the  yield  estimator 
can  be  obtained  with  stratified  sampling  techniques.  Also,  a  practical  method 
has  been  proposed  for  the  implementation  of  the  control  variate  method.  This 
method  may  be  the  most  useful  in  a  general  purpose  statistical  circuit 
analysis  program.  Finally,  a  new  approach  to  Monte  Carlo  based  yield 
maximization  has  been  investigated.  In  a  variety  of  passive,  active,  and 
switched-capacitor  filter  cirouits,  it  has  been  shown  that  a  quadratic  model 
can  be  used  to  approximate  the  responses  over  a  parameter  grid  several 
standard  deviations  in  length.  Thus,  the  circuit  yield  for  different  design 
values  can  be  estimated  using  correlated  sampling  and  the  quadratic  model. 
This  method  is  very  efficient  if  line  search  techniques  are  used  to  find  the 
maximum  yield  along  a  given  search  direction.  Several  filter  examples  show 
that  the  global  behavior  of  the  quadratic  model  is  good  and  that  it  is  very 
effective  for  yield  maximization  [16.24]. 
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17.1  CgaUfll  lad.  P.t.filiiOB  StM.mlfA  XSA  Sziilll  aadix  Imoexteet  Information* 

Uncertainties  arise  because  of  unknown  system  parameters,  unknown  signal 
enr ironnents »  and  hardware  tolerances.  Related  complexities  arise  in 
situations  inyolwing  multiple  system  performance  criteria  and  multiple 
decision  makers.  The  objective  of  this  project  is  to  gain  a  basic 
understanding  of  the  behavior  and  control  of  complex  systems.  During  the  past 
year,  attention  has  been  given  to  several  topics.  Highlights  are  summarised 
in  the  following. 

An  important  issue  in  the  design  and  analysis  of  control  and  estimation 

*  This  work  was  supported  by  the  Joint  Services  Electronics  Program  (U.S. 
Army.  U.S.  Navy,  and  U.S.  Air  Force)  under  contract  N00014-79-C-0424. 
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systems  is  the  asintenanee  of  systea  perforaaoee  in  the  presence  of  the 
inevitable  differences  between  the  design  aodel  and  the  trne  physical  systea 
(aee  [17.14]  for  a  review  of  robustness  and  sensitivity  results  for  nonlinear 
systeas) .  We  have  addressed  this  issue  using  three  distinct  approaches. 

The  first  approach  continues  our  work  on  the  analysis  of  feedback 
prop**ties  using  aeasures  of  gains  of  various  loop  transfer  functions.  For 
linear  systeas#  useful  aeasures  of  gain  are  the  singular  values  of  the 
transfer  function  aatris.  One  drawback  to  singular  value  analysis  is  that  the 
singular  values  aust  be  evsluated  at  the  true  plant  (as  opposed  to  the  noainal 
design  aodel)  to  provide  a  valid  indication  of  certain  feedback  properties. 
To  circuavent  this  difficulty#  the  singular  value  sensitivity  function  was 
developed  in  [17. 56 ,17 .16]  ss  an  analysis  tool  which  coapleaents  the  singular 
value  analysis. 

The  second  spproach  to  robustness  probleas  in  estiaation  and  control  uses 
a  ainiaaz  formulation.  The  basic  idea  is  to  aodel  plant  uncertainty  by 
asswaing  that  the  true  plant  is  known  only  to  be  an  eleaent  of  a  given  set. 
The  estiaator  or  control  systea  design  is  then  chosen  to  ainiaixe  the  worst 
ease  perforaance.  Niniaaz  probleas  of  this  type  were  solved  for  linear 
estiaation  [17.32]  and  control  [17.44,17.29,17.65]  probleas  for  stochastic 
systeas  with  white  noise  processes  whose  spectral  densities  are  uncertain,  and 
for  soalar  aeasureaent  noncausal  linear  estiaation  probleas  with  noise 
processes  whose  second  order  statistics  are  uncertain  [17.54,17.15]. 

The  final  approach  to  robustness  probleas  assuaes  that  the  plant  is  one 
of  a  finite  nuaber  of  given  aodels.  A  cost  functional  is  defined  for  each 
aodel,  and  the  controller  structure  is  fixed.  The  controller  paraaeters  are 
chosen  to  ainiaixe  a  weighted  sea  of  the  objective  functionals.  Although  this 
procedure  is  not  novel,  previous  solution  algorithas  have  required  that  an 
initial  stabilizing  value  for  the  controller  paraaeters  be  available.  In 
[17.26],  a  three  level  optiaization  algorithm  was  developed  which  eliminates 
this  problea. 

In  the  area  of  fundsaental  systea  theory  basic  to  oontrol  strategy 
development,  "strange  attractors"  have  been  investigated.  Currently  there 
seeas  to  be  no  unanimity  about  the  definitions  of  the  teras  "strange 
attractors"  or  of  "chaotic  behavior"  as  spplied  to  the  solutions  of 
deterainistic  differential  equations.  Perhaps  this  lack  of  agreement  arises 
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bMitit  tin*  in  important  diffirnoit  in  thi  be  Ret  i  or  of  aolntiona  which 
superficially  appear  to  bo  ainilar .  Vo  have  dona  o  nnabor  of  oolonlotiono  to 
toot  this  poaaibility;  wo  bawo  computed  solntiono  of  o  wariety  of  difforoatiol 
oqnationa  whiok  exhibit  "chaos",  including  (o)  tko  Lorenz  ay a ton  [Rl] ,  (b) 
Veda 'a  war i ant  of  tko  forced  Duffing  equation  [R2] ,  and  (c)  aororal  difforont 
exaaplea  of  an  antonoaons  nonlinear  feedback  ay a tea,  one  anck  being  tko  systea 
proposed  and  studied  by  Sparrow  [S3].  Next  wo  intend  to  study  tko  ay a tea 
recently  investigated  by  Oraxag  and  McLangklin  [R4J;  wo  think  tkoir  example 
will  be  particularly  useful  as  tkey  hare  coaputed  tke  long-term  correlations 
for  their  solutions. 

Ve  have  found  it  useful,  in  eoaiparing  tke  solutions  of  these  different 
systeas,  to  eaploy  tke  Lorenz  asp  ping  [SI],  that  is  to  plot  tke  next 
succeeding  local  maxima  versus  tke  iasMdiately  preoeding  local  maximum  of  any 
component  of  tke  (vector)  solution.  (One  obtains,  for  tke  Lorens  systea 
itself,  a  sharply  eusped  inverted  **vee*.)  Ve  have  found  that  some  of  tke 
solutions  exhibited  tke  "vee*  and  others  did  not.  One  would  expect  that 
purely  randoa  behavior  of  a  solution  would  yield  scatter  in  tke  plot,  and  not 
a  well  defined  curve. 

In  nonlinear  systeas.  the  tracking  or  synchronisation  of  a  nonlinear 
feedback  systea  to  multiple  quasi-periodic  inputs  has  been  obtained  using  both 
averaging  and  perturbation  aetkods.  One  is  able  to  eoapute  tke  tracking 
range,  at  least  approximately,  in  this  fashion.  These  results  will  be 
coapared  with  eoaputer  generated  numerical  solutions.  Finally,  a  new  proof 
for  tke  existence  of  periodic  solutions  of  tke  generalised  Lidasrd  equation 
has  been  found. 

Major  progress  also  has  been  a ado  in  control  of  singularly  perturbed 
systeas.  A  coaposite  slow-fast  control  is  designed  in  two  stages  to  guarantee 
stability  and  near-optiaal ity  of  nonlinear  deterministic  systeas.  This  result 
has  been  generalized  to  a  class  of  stochastic  control  problems. 

17.2  laplifatltlgl  Constrained  Decomposition  Ak4  Hierarchical  SaalSSli* 

In  large-scale  systeas.  control  tasks  nay  be  decomposed  and  hierarchical 
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lirili  U7  b«  iqottd  pa  ths  basis  of  analytically  ds rived  strstsgiss.  This 
project  is  devoted  to  the  fundamental  problem  of  iavestigstiag  sa  analytical 
framework  for  imeorporstiag  reslistie  eoastrsiats  of  a  computer  aetwork. 
Qeeatioas  eoaeeraiag  iaformatioa  flow*  time  scale*,  sad  loss  of  feedback  loops 
are  examiaed  ia  eoordiaatiag  the  multiple  eoatrollers  ia  a  aetwork. 

A  deeompositioa  methodology  was  developed  ia  [17.28]  which  allows  the 
structure  sad  local  eoavergeaee  properties  of  iterative  hierarchical 
algorithms  to  be  aaalyxed.  This  methodology  was  applied  ia  [17.27]  to  aaalyxe 
deeompositioas  of  weakly  coupled  systems,  sad  ia  [17.86]  to  study  the  effects 
of  problem  simplif ieatioa  oa  the  eoavergeaee  properties  of  hierarchical 
algorithms.  The  results  oa  weakly  eoupled  hierarchical  algorithms  were  used 
to  aaalyxe  the  eoavergeaee  of  several  solutioa  algorithms  for  the  exteaded 
algebraic  Riestti  equation  [17.67]. 

Aaother  area  of  iavestigatioa  eoaceras  Staekelberg  games  sad  iaeeative 
desiga.  Staekelberg  games  coastitute  aa  importaat  class  of  multi-person 
multi-criteria  decisioa  problems  whercia  there  is  a  hierarchy  ia  decisioa 
makiag,  some  players  being  ia  a  position  to  anaouaee  their  policies  ahead  of 
time  and  enforce  them  oa  the  remaining  players.  Vhea  the  underlying  decision 
process  is  dynamic  such  game  problems  become  mathematically  chal leaping,  since 
the  standard  available  techniques  and  approaches  of  optimisation  and  optimal 
eoattol  theory  are  not  direetly  applicable.  This  part  of  the  project  is 
directed  towards  developing  indirect  approaches  sad  a  general  theory  for  the 
solutioa  of  such  deeisioa  problems.  During  this  reporting  period  we  have  made 
important  advaaees  ia  this  ares,  which  are  summarised  ia  the  sequel. 

Our  first  set  of  results  is  eoueerued  with  the  derivation  of  bounds  oa 
the  performance  of  the  leader  (deeisioa  maker  at  the  top  of  the  hierarchy)  ia 
deterministic  discrete  and  time  dynamic  games  with  two  [17.7]  or  more  [17.8] 
levels  of  hierarchy;  sad  these  bouads  have  been  shown  to  be  tight  ia  many 
instances  [17.6,17.7,17.2].  Derivation  of  such  bouads  have  also  been  exteaded 
to  eoatiauous-time  dynamic  games  for  the  ease  when  the  leader  has  access  to 
saapled  state  iaformatioa  [17.44]. 

Our  second  set  of  results  is  concerned  with  the  derivation  of  policies 
(eoatrol  laws)  for  the  leader  that  aehieve  these  bouads.  In  this  regard,  a 
act  of  very  appealing  affine  policies  has  beea  obtaiaed,  and  their  existeaee 
has  bees  established  for  a  rather  general  class  of  such  dynamic  gssMs  [17.41]. 
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la  a  aore  recent  work  [17.52]*  m  hm  the  relationship  between 
saoh  daoiaioa  problaaa  and  incant iy a  daaiga  problena  that  arise  in  economics, 
and  haye  also  introdnoad  a  new  concept*  yix.  least  sensitiye  (robnst) 
iaeentiye  polieies  in  probleas  with  paraaetrie  uncertainties. 

Tat  another  set  of  resnlts  obtained  in  the  seope  of  this  project  involves 
stoehastie  Staefcelberg  gaaes  and  stochastic  incentive  design  probleas.  In 
[17.5]  we  have  obtained  the  eoaplete  solution  of  a  class  of  cont inmons- t iae 
stochastic  Stackelberg  differential  gaaes  by  solving  a  stoehastie  control 
problea  of  a  nonstandard  type.  In  [17.48]  we  have  studied  the  derivation  of 
optiaal  inoentive  seheaes  in  two-agent  stoehastie  decision  probleas  with  a 
hierarchical  decision  structure,  in  a  general  Hilbert  space  setting*  with  the 
underlying  inforaation  pattern  allowing  the  leader  to  acquire  both  eoaaon  and 
private  inforaation.  In  this  set-up,  we  have  shown  that,  under  soae  fairly 
wild  conditions,  there  exists  an  optiaal  incentive  policy  for  the  leader, 
which  is  affine  in  the  dynaaic  inforaation  and  generally  nonlinear  in  the 
static  (eoaaon  and  private)  inforaation. 

As  a  part  of  this  project,  we  have  also  investigated  the  applicability  of 
the  feedback  Stackelberg  solution  concept  (developed  priaarily  for  diserete- 
tiae  probleas)  in  continuous-tine  Stackelberg  differential  gaaes,  and  have 
reported  soae  preliminary  results  in  [17.51].  More  extensive  results 
regarding  this  problea  will  be  reported  in  the  next  period. 

17.3  atltlflitl  ISL  SaaalUL  Svsteas  Xfi£  All  id  Atraitut  Avionics* 

Whenever  aodel  uncertainty  is  present  or  a  range  of  operating  conditions 
is  anticipated,  engineering  analysis  and  design  aust  deal  with  questions  of 
sensitivity.  Perforasnces  of  aanuf actured  components  are  aoeessarily 
specified  with  nonxero  tolerances.  The  par sac ter  values  which  characterise 
these  components  are  therefore  uncertain.  These  paraaeter  values  aay  change 
during  the  operating  lifetimes  of  the  components  due  to  aging  and  due  to 
changes  in  envirosaental  conditions,  such  as  pressure  and  teaperature.  la 
addition,  mathematical  aodel s  used  for  analysis  and  design  of  actual  systeas 
cannot  possibly  lead  to  predicted  perforaanee  which  exactly  aatches  the 
performance  of  the  actual  systeas.  It  is  useful  to  regard  soae  of  the 


*  This  work  was  supported  by  the  U.S.  Air  Force  under  grant  AFOSR-7 8-3633 . 
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pttiMttri  of  the  aodols  os  uncertain  la  order  to  asks  theae  node Is  son 
realistic.  The  eoabined  offsets  of  paraaeter  uncertainty-'  oa  overall  systea 
behavior  are  of  priaelpal  consideration  ia  any  systea  desi^a. 

Oar  objective  here  has  been  to  advance  the  earreat  state  of  knowledge 
aader  both  cooperative  sad  aoneooperative  solution  concepts.  Ia  [17.2], 
several  new  resalts  have  beea  reported  concerning  stochastic  teaa  probleas 
with  general  partially  nested  iafonatioa  patterns  and  also  the  Nash 
equilibriua  solutions  of  stochastic  sons ero- sub  dynaaie  gases  with  quasi— 
classical  infonaation  patterns.  Recently,  no  have  developed  an  eqnilibrina 
theory  for  a  aore  general  elsss  of  aulti-persoa  anlti-eriteria  stochastic 
decision  probleas  wherein  the  decision  ackers  have  different  subjective 
probability  aeasures  on  the  uncertain  quantities  [17.47],  i.e.,  for  probleas 
which  do  not  adait  a  season  underlying  probability  space.  For  such  decision 
probleas,  we  have  introduced  the  concepts  of  locally  and  globally  stable 
eqnilibrina  solutions,  and  have  atudied  the  existence  and  uniqueness  of  such 
equilibria  under  certain  conditions  on  the  loss  functionals  and  the  subjective 
probability  ae asuree.  Certain  explicit  results  and  analytic  solutions  have 
also  beea  Obtained  and  reported  in  [17.47]. 

On  another  topic,  an  iapertaat  outgrowth  of  the  singular  perturbation 
aethodology  has  been  the  analysis  of  adaptive  systeas  with  reduced  order 
aodels  [17. 5g,17. 19,17 .It, 17. 14,17.24,  and  others].  The  error  ia  aodel  order 
is  due  to  the  presence  of  fast  uaaodeled  par aae tars.  A  "doainant  richness" 
condition  has  been  introduced  under  which  adaptive  identities  and  observers 
can  be  aade  robust  with  respect  to  paraae tries.  To  extend  this  aethodology  to 
the  ease  of  feedback  adaptive  control  the  robustness  of  hyperstability 
properties  has  been  analysed. 

Finally,  as  an  application  of  chained  aggregation,  a  study  of  "near" 
unobservability  has  been  aade  [17.70].  This  concept  has  been  developed 
geoaetrically  in  terns  of  the  canonical  angles  between  an  observable  subspace 
L  and  an  A-invariant  subspace  Y.  The  role  of  eigenvalue  separation  in  this 
problea  has  been  clarified,  leading  to  a  contact  with  singular  perturbation 
theory. 
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17 .4  Ut«i  Sol*  SritMi* 

V*  highlight  two  eliiNi  of  roooat  aeoompliahments  on  |<mt.  Ou  eloos 
of  result*  portoiao  to  the  well-poaednos*  of  singularly  par tor bod  games 
[17.34].  For  praotieal  roaaoaa  It  &a  doalrablo  to  employ  reduetd  ordar  aodola 
for  the  datoraiaatloa  of  game  strategic*.  Vo  bar*  demonstrated  that  tba 
aataral  alagolar  yortarbatioa  approaab  loada  to  ill-posed  aolotioaa  bat  wa 
alao  daaonatratad  bow  to  obtaia  well-posed  radooad  ordar  aodola.  Tba  bo y  to 
tb*  aneeaaafol  daralopaaat  of  aaob  rodaead  ordar  aodola  la  tbo  proaorwatioa  of 
tbo  iaforaation  atrootore. 

m 

A  aoeoad  elaaa  of  onr  raaolta  va  wiab  to  highlight  la  tbo  aao  of 
incentives  to  achieve  global  taaa  optlaal ity  for  tb*  loador  [17.37,17.2]. 
Basically,  tb*  ldaa  la  for  tb*  loador  to  apply  a  stratagy  which  ladncaa  tb* 
followara  to  optlalaa  tb*  loador* a  objaetlwa  wblla  optimising  tbalr  own.  “ 

Anotbar  aignlf leant  aeeoapliabaont  la  tb*  two  tla*  acal*  aatbodology  r 

[17.22,17.39,17.69,17.33]  for  aodallag  of  largo  aeal*  dynamic  networks  snob  as 
power  ays tans  and  Markov  obalna.  A  fundamental  rolatlonsblp  botwaen  weak 
coupling  and  tin*  scale*  baa  boon  revealed  and  need  for  systanatic  “ 

detorminatiott  of  oobarent  areas  as  local  fast  aubsystoms  and  tbeir  interact iva 

r 

core,  as  tba  slow  aggregate  subsystem.  Ibis  type  of  modeling  naturally  load* 
to  a  schema  for  fast  decentralised  and  slow  coordinated  control,  wbicb  baa 
bean  developed  for  Markov  chains  and  queueing  networks.  * 

Finally,  a  detailed  ga one trio  analysis  of  chained  aggregation  baa  boon 
mad*  [17.23].  This  baa  lad  to  a  general  version  of  three-ooatrol-eonponent 
design,  a  biararcbical  control  stratagy  for  large  scale  systems.  This 
approach  exploits  tb*  a  priori  information  structuring  present  in  tb*  system 

e> 

and  system  output  measurements.  Connections  with  reduced  order  modeling, 
especially  the  role  of  residual  feedback,  have  been  examined. 


*  This  work  was  supported  by  the  National  Science  Foundation  under  grant 
BCS  79-19396. 
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17 .5  £gxu  System  Normal  uA  Boourltv  State  Aimmat* 

Results  Itn  been  obtained  is  tbe  two  areas  of  to 1 tags  collapse  aad 
maximum  spates  loadability.  As  a  rariable  allocation  problem,  Toltage 
collapse  has  been  found  to  be  directly  related  to  generator  excitation  and 
transformer  tap  limits.  Tbe  normal  control  algorithm  for  dispatching 
transformer  taps  during  contingencies  can  result  in  extreme  demands  for 
reactive  power  from  sources  with  finite  reserves.  Alternative  approaches  to 
voltage  control  are  currently  being  considered  to  minimise  tbe  possibility  of 
voltage  collapse. 

Direct  transient  stability  methods  are  being  integrated  with  conventional 
maximum  system  loadability  methods  to  replace  exhaustive  time  domain 
simulations.  Tbe  "first  cut"  fast  transient  stability  calculation  was 
developed  to  si low  the  dynamic  constraint  to  be  enforced  in  the  contingency 
evaluations  required  for  the  maximum  loadability  calculation.  Tbe  method 
employs  an  approximate  Potential  Energy  Boundary  Surface  (PEBS)  approach 
together  with  linear  load  flow  algorithms. 

17.6  ggltr  System  Reduced  Order  Mode lint* 

Tbe  traditional  models  used  to  analyse  large  aeale  power  systems  contain 
many  inconsistencies  resulting  from  heuristic  order  reduction  and  model 
simplification.  Several  of  these  models  have  been  successfully  derived  from  a 
singular  perturbation  approach  while  others  have  been  shown  to  be  void  of 
mathematical  justification.  In  as ;h  case,  improved  models  can  be  obtained 
through  higher  order  correctious  to  the  singularly  perturbed  model,  lorh  is 
continuing  to  develop  a  consistent  set  of  reduoed  order  models  for  all  of  the 
components  and  interconnections  in  a  large  scale  power  system. 


*  This  work  wss  supported  by  industrial  funding  through  the  University  of 
Illinois  Power  Affiliates  Program  81-82. 
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18.1  &£  PlitstlMM  in  falUfttowaiignsl.  Sim 1  Emuilu* 

This  research  project  spans  the  areas  of  multisensor  array  processing, 
linear  space  varying  processing,  and  nonlinear  processing  using  order 
statistics.  In  array  processing  we  have  developed  a  unified  treatment  of 
spotlight  node  synthetic  aperture  radar  (SAR)  and  conputor-aided  tomography 
(CAT).  Ve  have  shown  that  Fourier  domain  reconstruction,  based  on  the 
projeetion'*sliee  theorem  from  CAT,  is  also  the  underlying  principle  in 
spotlight  mode  SAR.  In  particular,  the  signal  recorded  at  eaeh  SAR 
transmission  point  is  approximately  the  Fourier  transform  of  a  central 
projection  of  the  imaged  ground  area.  A  number  of  important  issues  involving 
resolution,  sampling  rates,  waveform  curvature,  the  Doppler  effect,  and  motion 
compensation  have  been  examined  within  the  context  of  this  new  SAR 
interpretation.  In  addition,  we  have  modified  the  convolution-back-proj ection 
algorithm  from  CAT,  and  applied  it  to  the  SAR  problem.  Preliminary 
indications  are  that  this  algorithm  provides  higher  resolution  than  current 


*  This  woa»  <rss  supported  by  the  Joint  Services  Electronics  Program  (U.S. 
Army,  U.S.  Navy,  U.S.  Air  Force)  uader  contract  N00014-79-C-0424. 
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aethods  of  SAR  processing. 

York  is  also  underway  on  tke  anbjeot  of  digital  beamf oraing  for  aonar. 
In  a  praetioal  situation,  tke  maker  of  aonar  aenaora  nay  be  reatrieted  by 
econoaic  or  pbysieal  constrainta.  This  difficulty  can  be  partially  overcoae 
by  using  signal  extrapolation.  Ve  bave  applied  a  band liaited  extrapolation 
algoritka  to  this  problea  and  ckaraeterixed  ita  perforaance  botk  analytically 
and  through  coaputer  siaulation.  For  a  ainuaoidal  aignal  it  ia  posaible  to 
extrapolate  to  two  or  three  tines  the  length  of  the  physical  array.  The  beaa 
perforaance  ia  thus  enhanced  and  is  coapatible  with  that  obtained  using  aany 
sure  sensors. 

In  the  area  of  linear  space  varying  processing  we  have  derived  the 
ainiaua  permissible  sampling  rate  for  the  digital  implementation  of  an  analog 
shift  variant  systea.  The  necessary  rate  is  easily  coaputed  froa  paraaeters 
of  the  bi-frequency  asp  of  the  desired  analog  response.  Potential 
applications  include  digital  restoration  of  iawgea  degraded  by  linear 
spatially  varying  operations. 

Ve  have  aade  further  progress  ia  the  area  of  nonlinear  processing  using 
order  statistics.  A  generalized  aedian  filter  has  been  studied  whose  output 
is  a  weighted  sua  of  the  input  order  statistics.  An  explicit  foraula  has  been 
derived  for  the  optisul  weight  values  in  the  order  statistic  filter  (OSF) , 
assuming  white  background  noise.  Ve  have  filtered  several  iaages  plus  noise 
with  a  linear  averager,  a  aedian  filter,  and  the  optimal  OSF  for  a  constant 
background.  Ve  found  that  the  OSF  provides  a  coaproaise  between  the  smoothing 
properties  of  the  linear  filter  and  the  edge  sharpness  offered  by  a  aedian 
filter.  Further  generalizations  of  the  aedian  filter  are  currently  being 
investigated. 

• 

18.2  Support-Limited  Sianal  ni  Iaaxt  Extrapolation* 


One  of  the  aost  important  and  difficult  questions  in  signal  and  iaage 
processing  is:  How  can  we  reconstruct  signals  and  iaages  froa  incomplete  data 
(observations)?  Ve  have  eabarked  on  a  basie  research  projeot  to  study  a 
partieular  problea  in  this  area,  that  of  support-limited  extrapolation.  In 


*  This  work  was  supported  by  the  Joint  Services  Electronics  Program  (U.S. 
Aray ,  U.S.  Navy.  U.S.  Air  Force)  under  eontraet  N00014-79-C-0424. 
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the  one-dimensional  form,  the  problem  is  the  following:  Given  the  Fonrier 
transform  F(u)  of  a  function  f(t)  in  a  finite  interval,  determine  f(t)  under 
the  assumption  that  f(t)  is  of  compact  support.  In  the  mid  and  late  1970's 
two  efficient  algorithms  (Gerchberg-Papoulis,  Cadzow-Sabri-Steenart)  were 
proposed,  which  looked  very  promising.  However,  up  to  now  no  successful 
applications  of  these  algorithms  to  real-world  problems- have  been  swde.  We 
think  the  reasons  are  twofold. 

(1)  No  one  has  studied  the  effects  of  discretization  of  these  algorithms. 

(2)  No  one  has  modified  these  algorithms  or  come  up  with  new  algorithms  which 

take  account  of  the  noise  in  the  observed  data  in  an  effective  way. 

We  have  been  concentrating  our  effort  on  these  two  aspects.  Some  preliminary 
results  have  already  been  obtained.  In  particular,  we  have  generalized 
Cadzow's  algorithm,  and  shown  that  the  solutions  from  the  discrete  algorithm 
converge  to  the  solution  from  the  continuous  algorithm  under  certain  rather 
weak  conditions. 

18.3  Lsartto.llP  Pigiul  Filter  Structures* 

Error  spectrum  shaping  (ESS)  is  a  technique  that  uses  quantizer  feedback 
to  reduce  roundoff  noise  and  limit  cycle  amplitudes  in  digital  filters. 
Nearly  all  past  work  on  ESS  has  dealt  with  second-order  filters.  In  our 
recent  work  we  have  shown  how  to  choose  the  ESS  coefficients  optimally  for 
high  order  filters  composed  of  cascaded  second-order  sections.  We  have 
compared  the  noise  gains  of  optiswl  and  suboptimsl  ESS  structures,  and  the 
section-optimal  structure  derived  using  the  Mullis  and  Roberts  state-space 
approach;  the  ESS  structures  performed  well  in  comparison.  In  addition,  a 
heuristic  strategy  for  ordering  the  second-order  sections  has  been  developed 
based  on  the  ideas  underlying  ESS.  This  strategy  further  reduces  roundoff 
noise,  especially  for  elliptic  filters. 

We  have  also  applied  ESS,  in  a  much  different  role,  for  the  design  of 
short  wordlength  FIR  filters.  Instead  of  rounding  infinite  precision  filter 
coefficients,  we  have  used  an  ESS  quantizer  to  distribute  the  error  due  to 


This  work  was  supported  by  the  Joint  Services  Electronics  Program  (D.S. 
Army,  D.S.  Navy,  U.S.  Air  Force)  under  Grant  N00014-79-C-0424. 
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rounding  nonun if ormly  throughout  the  frequency  bend.  As  an  exaaple.  for  cesee 
in  which  the  stopband  response  aust  be  very  low,  the  ESS  technique  can 
effectively  place  into  the  pasaband  nearly  all  of  the  error  due  to  coefficient 
rounding.  We  are  presently  concentrating  on  several  factors  involved  in  the 
design  of  the  ESS  feedback  filter. 

18.4  fiftiiia  pi  IH  Digital  ElUpn  filing.  Xnltlnle  Criterion  Optimisation* 

We  are  conducting  an  investigation  of  the  applicability  of  the  concepts 
and  algorithms  of  multiple  criterion  optimization  (MOO)  to  the  problem  of 
optimally  trading  off  magnitude  and  phase  response  in  I£R  digital  filter 
design.  Particular  attention  is  being  paid  to  the  formulation  of  the 
sisiultaneous  magnitude  and  phase  design  problem.  Results  to  date  are 
extremely  encouraging.  As  compared  with  an  elliptic  design,  the  MCO  approach 
can  produce  filters  having  magnitude  responses  that  are  virtually  as  good  and 
phase  responses  that  are  much  more  linear. 

18.5  Nggfrtr  Thgqgqjtig  fiMfitatl  las.  ihs.  Realization  at  Failure 
tollltHt  Signal  Pgqqggmg  in  VLSI  Circuit  Technology** 

This  research  involves  the  study  of  number  theoretic  concepts  which  can 
improve  speed,  reduce  hardware  complexity,  and  provide  self-checking 
performance  in  specialized  digital  processors  such  as  digital  filters  or  FFT 
processors.  In  particular,  residue  number  system  (RNS)  concepts  appear  to 
provide  many  features  that  are  required  to  integrate  custom  designed  functions 
in  VLSI  circuits  [Rl, 18.15,18.16] . 

Three  significant  results  have  been  obtained  that  alleviate  same  of  the 

practical  difficulties  in  using  RNS  arithmetic  in  VLSI  processors: 

• 

(1)  the  design  of  a  high  speed  pipe-lined  RNS  error  checker  resulted  from  a 
mathematical  study  of  the  error  location  algorithm  [Rl], 


*  This  work  was  supported  by  the  National  Science  Foundation  under  Grant 
EC8  80-07075. 

**  This  work  was  supported  by  the  National  Science  Foundation  under  Grant 
ENG-79-01686. 
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(2)  a  new  error  location  algorithm  based  on  tbe  concept  of  expanded 
projection  was  derived  [R2] ,  and 

(3)  an  architecture  for  a  standard  coaputational  eleaent  was  developed  for 
the  VLSI  realization  of  digital  processors  designed  with  aodnlar 
arithmetic  [R3], 

Also*  RNS  product  codes  have  been  studied  to  determine  if  they  offer  more 
desirable  properties  for  VLSI  realization.  A  hardware  architecture  for  an 
error  checker  based  on  an  RNS  product  code  has  been  designed  [18.15].  Further 
studies  will  reveal  how  the  special  features  of  a  product  code  can  improve 
VLSI  architectures  for  fanlt  tolerant  performance. 
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An  ongoing  study  is  continuing  to  relate  cliaate  anoaalies  in  one  part  of 
the  world  at  a  prior  tine  to  crop  growing  conditions  in  other  regions  of  the 
world  at  soae  later  tiaes.  A  nuaber  of  interesting  oliaate-crop  patterns  have 
been  found  which  have  a  high  probability  of  repeating  theaselwes  every  3-7 
years. 

The  cliaate-orop  patterns  can  also  be  need  as  a  naans  of  crop 
forecasting.  Various  types  of  tiae  series  analysis  have  been  developed  for 
recognition  of  these  cliaate-orop  patterns.  Soae  of  the  crop  tiae  series 
contain  data  starting  in  1866.  A  paper  describing  the  relationship  of  crops 
in  Australia  and  Indonesia  to  cl  Lasts  anoaalies  in  the  eastern  tropical 
Paoific  Ooean  has  been  subaitted  for  publication. 
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